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Введение 


В  начале  90-х  годов  было  введено  в  обиход  специалистами  понятие 
«доказательная  медицина»,  концепция  которой  подразумевает 
добросовестное,  точное  и  осмысленное  использование  лучших  результатов 
клинических  исследований  для  выбора  лечения  конкретного  больного. 
Подобный  подход  позволяет  снизить  уровень  врачебных  ошибок, 
облегчить  процесс  принятия  решения  для  практических  врачей  и 
администрации  лечебных  учреждений,  а  также  уменьшить  расходы  на 
здравоохранение  и  использовать  сэкономленные  средства  для  претворения 
в  жизнь  социально  ориентированных  медицинских  проектов  [1]. 

Согласуясь  с  принципами  доказательной  медицины,  врач  обязан 
критически  подходить  к  решению  возникающих  перед  ним  проблем  и,  как 
правило,  самостоятельно  предпринимать  те  или  иные  действия. 
Современную  медицину  невозможно  представить  себе  без  использования 
точных  и  надежных  методов  анализа. 

Обоснованность  принимаемых  решений  в  настоящее  время  в 
мировой  практике  оценивается  с  помощью  статистических  методов. 

Однако  проблема  состоит  в  том,  что  в  большинстве  случаев  объекты 
медицинских  исследований  являются  сложными  медицинскими  системами 
с  множеством  элементов  в  своем  составе  и  связей  между  ними, 
определяющих  внутреннюю  структуру  системы.  Такие  системы  (объекты, 
процессы)  формируются  и  функционируют  под  воздействием 
многочисленных  изменчивых  и  взаимодействующих  между  собой 
факторов  среды  [2  -  4].  Следовательно,  для  применения  многомерного 
статистического  анализа  имеются  все  основания.  На  сегодняшний  день 
разработана  теория  и  практика  адекватного  кодирования  медицинских 
данных.  Современные  технические  и  программные  средства 
информатизации  обеспечивают  возможность  создания  баз  данных  в 
объемах,  достаточных  для  получения  достоверных  результатов.  Наличие 
пакетов  прикладных  программ  по  математико-статистической  обработке  и 
анализу  данных  позволяют  применять  многомерные  математические 
методы  и  процедуры  специалистами-медиками.  Создаваемые  на  этой 
основе  многомерные  математико-статистические  модели  сложных 
медицинских  систем  позволяют  перейти  от  интуитивных  представлений 
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врача  о  медицинских  процессах  и  системах  к  математически 
обоснованному  количественному  анализу. 

Современные  компьютерные  программы  призваны  оказать 
эффективную  помощь  исследователю,  однако,  при  работе  с  компьютером 
возникают  проблемы  интерпретации  полученных  результатов. 

Наука  означает  знание.  Понимание  и  объяснение  являются 
непрерывными  процессами,  получающими  развитие  по  мере  того,  как 
выявляются  все  новые  и  новые  взаимосвязи  между  различными  научными 
законами. 

В  монографии  просто  и  наглядно  изложены  результаты  научных 
исследований,  проведенных  автором  совместно  с  врачами- 
исследователями. 

На  реальных  медицинских  примерах  статистическими  методами 
(множественный  регрессионный  анализ,  дискриминантный  анализ,  анализ 
временных  рядов)  была  исследована  структура  медицинских  данных. 
Исследование  и  анализ  данных  позволили  решить  задачу  идентификации 
состояния  больных  и  прогнозирования.  В  качестве  инструментария 
исследований  использованы  пакеты  прикладных  программ 
8ТАТІ8ТІСА  10,  8Р88  21.  Особое  внимание  уделено  интерпретации 
полученных  результатов. 

В  первой  и  второй  главах  представлены  результаты  применения 
дискриминантного  анализа  при  определении  хронической  сердечной 
недостаточности  (ХСН)  и  дисфункции  миокарда  у  подростков.  Для 
каждого  из  заболеваний  выделено  небольшое  количество  клинико¬ 
лабораторных  показателей,  по  которым  можно  с  высокой  достоверностью 
определить  степень  тяжести  состояния  больного.  Показана  возможность  и 
целесообразность  использования  дискриминантного  анализа  для 
идентификации  состояния  больных  с  указанными  заболеваниями. 

В  третьей  главе  рассмотрено  применение  бинарной  логистической 
регрессии  для  определения  сердечно-сосудистого  риска  у  пациентов  с 
почечной  недостаточностью.  Были  изучены  факторы  риска  сердечно¬ 
сосудистых  заболеваний,  среди  которых  выделены  наиболее  значимые, 
оказывающих  наиболее  значительное  влияние  на  прогрессирование 
хронической  ишемической  болезни  сердца  в  виде  ишемической 
дилятационной  кардиомиопатии. 
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Четвертая  глава  посвящена  вопросам  прогнозирования 
заболеваемости  расстройствами  психики  и  поведения  непсихотического 
характера  среди  сельского  населения  Украины.  Посредством  анализа 
временного  ряда  построены  модели  ряда,  по  которым  можно  делать 
прогноз  интересующего  показателя  на  будущий  временной  интервал  в 
различных  регионах  Украины. 

Книга  рассчитана  на  широкий  круг  читателей  -  студентов, 
аспирантов  и  специалистов,  использующих  анализ  данных  в  медицинских 
исследованиях. 

Результаты  анализа  данных  и  построенные  вероятностно¬ 
статистические  модели  могут  быть  использованы  при  принятии 
медицинских  решений  по  стратегии  и  тактике  лечения  больных. 
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ГЛАВА  1.  ОПРЕДЕЛЕНИЕ  РАННЕЙ  СТАДИИ 

ХРОНИЧЕСКОЙ  СЕРДЕЧНОЙ  НЕДОСТАТОЧНОСТИ  У 
ПОДРОСТКОВ  С  ИСПОЛЬЗОВАНИЕМ  8Р88 

Изучение  связей  основных  клинических  параметров  и  их 
комбинаций  является  важнейшей  задачей  теоретической  и  практической 
медицины.  В  настоящее  время  медицина  все  чаще  обращается  в  своих 
исследованиях  к  достижениям  математических  наук.  Не  вызывает 
сомнения  тот  факт,  что  статистическое  моделирование  основных 
жизненных  функций  в  норме  и  патологии  имеет  большое  теоретическое  и 
практическое  значение.  Оно  позволяет  объяснить  суть  патологического 
процесса,  правильно  оценить  прогноз  заболевания  и  разработать 
эффективные  лечебно  -  профилактические  мероприятия. 

Рост  сердечно-сосудистых  заболеваний  в  мировом  масштабе  и  их 
«омоложение»  в  последнее  время  является  общеизвестным  фактом. 
Большинство  проблем  общества  на  современном  этапе  развития  тем  или 
иным  образом  связано  с  демографическим  кризисом  и  инвалидизацией 
лиц  трудоспособного  возраста,  ростом  числа  больных  с  различными 
нарушениями  в  работе  сердечно-сосудистой  системы  и  снижением 
продолжительности  жизни. 

В  настоящее  время  профилактика  формирования  и  прогрессирования 
ХСН  является  приоритетным  направлением  кардиологии.  В 
педиатрической  практике  большой  удельный  вес  занимают  заболевания 
сердца,  которые  включают  врожденные  аномалии  (врожденные  пороки 
сердца,  множественные  малые  структурные  аномалии  сердца)  и 
приобретенные  патологии  (кардиты,  вторичные  кардиомиопатии  и  т.д.), 
на  фоне  которых  возможно  формирование  ХСН  уже  в  раннем  возрасте  [5  - 
6]. 

Как  известно,  развитие  ХСН  является  патофизиологически 
необратимым  процессом,  поскольку  механизмы  ее  формирования, 
выполняя  адаптационную  функцию  на  начальных  этапах,  постепенно 
становятся  дезадаптационными  факторами,  приводящими  к 
прогрессированию  процесса  [7].  Основные  звенья  данного 
патологического  процесса  включают  активацию  нейрогуморальных  систем 
регуляции,  системы  цитокинов  и  оксидативного  стресса,  ремоделирование 
миокарда,  которые  тесно  связаны  между  собой.  В  то  же  время,  на  ранних 
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этапах  при  отсутствии  выраженных  клинических  проявлений  (одышка, 
отечный  синдром)  диагностика  ХСН  затруднена. 

Наряду  с  общеизвестными  медицинскими  методами  диагностики 
ХСН  и  определения  тенденции  развития  систолической  дисфункции 
миокарда  все  чаще  исследователи  обращаются  к  математическому 
описанию  данных  процессов. 

В  детерминированных  моделях  входные  значения  моделируемого 
процесса  представляются  в  виде  аналитических  закономерностей,  для 
которых  может  быть  получен  единственный  результат.  Эти  модели  могут 
отображать  как  детерминированную,  так  и  вероятностную  системы,  но  при 
этом  описание  последней  бывает  неполным,  так  как  модель  в  таких 
случаях  является  упрощенной  [8]. 

Стохастические  модели,  в  которых  информативные  признаки, 
характеризующие  состояние  моделируемого  объекта,  представлены 
случайными  величинами  и  связаны  стохастическими  (случайными, 
нерегулярными)  зависимостями,  позволяют  решить  эту  проблему. 

Так  как  параметры,  характеризующие  систолическую  дисфункцию 
миокарда,  могут  принимать  одно  из  множества  значений,  причем 
появление  значения  до  измерения  нельзя  точно  предсказать,  то  для 
разработки  диагностической  модели  систолической  дисфункции  миокарда 
у  подростков  с  учетом  изменений  морфофункциональных  характеристик 
сердца,  показателей  систем  нейрогуморальной  регуляции, 
иммуновоспалительной  активации  и  свободнорадикальных  процессов 
необходимо  применять  стохастические  методы  моделирования. 

Дискриминантный  анализ  дает  возможность  быстро  и  качественно 
классифицировать  исследуемые  объекты  и  синтезировать  адекватную 
линейную  математическую  модель.  Данный  метод  математического 
анализа  имеет  такие  преимущества:  учитывает  вариабельность 
используемых  в  моделировании  параметров,  рассматривает  совокупность 
всех  клинических  и  параклинических  показателей,  взятых  со  своими 
коэффициентами,  которые  отражают  удельный  вес  влияния  каждого 
показателя  на  постановку  диагноза  [9  -  11]. 

Метод  диагностики  хронической  сердечной  недостаточности  у 
подростков  на  ранней  стадии  развития  включает  в  себя  ряд  этапов. 
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Этап  1.  Сбор  информации  о  состоянии  организма  пациента, 
включающий  определение  морфофункциональных  характеристик  сердца, 
путем  проведения  эхокардиографии  в  покое  и  после  физической  нагрузки, 
показателей  систем  нейрогуморальной  регуляции,  аутоиммунным 
активации  и  свободнорадикальных  процессов. 

Были  проанализированы  данные  159  подростков  с  заболеваниями 
сердечно-сосудистой  системы.  Все  пациенты  по  клиническим  критериям 
были  разделены  на  3  группы:  1-я  группа  -  27  подростков  с  заболеваниями 
сердечно-сосудистой  системы  без  признаков  ХСН  (17%);  2-я  группа  -  65 
подростков  с  заболеваниями  сердца  и  начальной  стадией  ХСН  (40,9%);  3-я 
группа  -  контрольная,  которую  составили  67  практически  здоровых 
сверстников  (42,1%). 

Так  как  наименьший  процент  составили  подростки  с  заболеваниями 
сердечно-сосудистой  системы  без  сердечной  недостаточности,  то  объем 
выборки  для  создания  математической  модели  определялся  с  учетом  этого, 
наименьшего  показателя. 

Все  пациенты  рассматривались  как  совокупность  объектов  с 
варьирующими  количественными  и  качественными  характеристиками,  на 
основании  которых  определялась  группа,  к  которой  принадлежит  объект. 
Это  позволило  предсказать  группы  (наличие  /  отсутствие  сердечной 
недостаточности)  для  новых  объектов  из  той  же  совокупности. 

Метод  дискриминантного  анализа,  как  и  вся  группа  методов, 
основанных  на  наблюдениях,  предполагает  использование  определенной 
терминологии.  Классифицируются  объекты,  а  различные  группы,  к 
которым  необходимо  эти  объекты  относить,  чаще  всего  называются 
классами.  Алгоритм  отнесения  объектов  к  разным  классам  называется 
решающим  правилом. 

Нами  рассматривались  две  выборки,  по  одной  из  которых  строили 
решающее  правило  (обучающая),  а  по  второй  проверяли  качество  работы 
этого  правила  (контрольная). 

Этап  2.  Кодирование  качественных  и  количественных  показателей. 

Для  определения  значимых  для  формирования  диагноза  свойств, 
было  проанализировано  73  показателя.  Все  показатели  были  закодированы 
и  поставлены  в  соответствие  73-мерному  вектору,  который  учитывает 
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отсутствие,  наличие,  направленность  и  величину  каждого  показателя. 
Кодировочная  таблица  представлена  в  табл.  А.1  (Приложение  А). 

Ни  одна  переменная  не  была  линейной  комбинацией  других 
переменных.  Т.е.  мы  не  пользовались  суммой  переменных  или  их  средним 
арифметическим  совместно  с  самими  переменными. 

Математическая  обработка  результатов  проводилась  с 
использованием  и  пакета  прикладных  программ  8Р88  21  для  \Ѵі псіохѵз  . 

Для  ввода  всех  закодированных  данных  в  компьютер,  после  запуска 
8Р88  для  \УіпсІо\У8,  с  помощью  двойного  нажатия  левой  кнопки  мыши  на 
значке  8Р88  открывается  редактор  данных  8Р88  (рис.  1.1). 


Рисунок  1.1-  Редактор  данных:  просмотр  данных 

В  редакторе  данных  двойным  нажатием  на  ячейках  с  надписью 
«пер»  перешли  в  режим  просмотра  переменных,  который  обеспечивает 
редактор  данных  (рис.  1.2).  Здесь  последовательно,  строка  за  строкой 
определялись  необходимые  переменные.  Для  этого  двойным  нажатием  в 
таких  полях,  например,  как  «Имя»,  «Тип»,  «Ширина»  и  т.д.,  задавалась  вся 
необходимая  информация  для  вводимых  переменных. 

При  кодировании  качественной  переменной  8П,  по  которой  мы 
группировали  исследуемых  пациентов,  учли  что  этот  показатель  может 
принимать  три  значения:  1-  «Без  сердечной  недостаточности»,  2-  «1-я 
стадия  сердечной  недостаточности»,  3  -  «Контрольная  группа». 
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После  открытия  диалогового  окна  «Метки  значений»  (рис.  1.3)  метки 
значений  были  определили  следующим  образом: 

1)  вначале  ввели  в  поле  «Значение»  число  «1»  и  нажали  клавишу 
<ТаЪ>; 


Рисунок  1.2  —  Диалоговое  окно  редактирования  переменных 

2)  ввели  в  поле  «Метка»  текст  «Без  сердечной  недостаточности»; 

3)  нажали  на  кнопке  «Добавить»,  после  чего  метка  значения  была 
добавлена  в  список; 

4)  повторили  эти  действия  для  значений  «2»  —  «1-я  стадия 

сердечной  недостаточности»  и  «3»  —  «Контрольная  группа». 

Максимально  допустимая  длина  метки  значения  составляет  60  знаков. 

Таким  образом,  закодировали  все  качественные  показатели. 


Этап  3.  Расчет  необходимого  объема  выборки  пациентов. 

С  учетом  предельной  ошибки  выборки  Аа  =5  —  6% .  При  надежности 
у  =  0,95  и  рассчитанной  предельной  погрешности  необходимый  объем 
выборки  составляет: 

Г2-со-(1-со)  ,  ,0  . 

п..  = - ^ - -  =  148  ч-  214  человек, 

д2 

со 

где  1  -  значение  коэффициента  доверия  Стьюдента  при  заданном 
значении  у;  со  -  больные  с  заболеваниями  сердца  без  сердечной 
недостаточности,  %. 
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Этап  4.  Описательный  анализ  данных. 

На  этом  этапе  определяли  имеет  ли  каждая  переменная  нормальное 
распределение  при  фиксированных  остальных.  Если  экспериментальные 
данные  не  вполне  удовлетворяют  этим  предположениям,  то  статистические 
выводы  не  будут  точным  отражением  реальности. 


Рисунок  1.3  —  Диалоговое  окно  «Метки  значений» 

Выбрали  в  меню  «Анализ»  пункт  «Классификация» 
«Дискриминантный  анализ»  (рис.  1.4). 


- 

Файл  Правка  Вид  Данные 

Преобразо  в  ать 

Анализ  Графика  Сервис  Дополнения  Окно  Справка 

^  Н  ^ 

ИТ  пита 

в-  И* 

#4  ♦( 

Отчеты  ► 

1  :  зп 

1.0 

Описательные  статистики  ► 

1 

пот 

сйадпог 

1  0' 

Таблицы  ► 

ѵоггазі  екд  сНзз  роѵ 

іеіа 

ѴуНіГ  Р  2  | 

1 

6  1 

Анализ  РРМ  ► 

15 

101 

1  .54 

2.17 

2 

20 

2 

Сравнение  средних  ► 

16 

54 

1  .90 

2.75 

3 

57 

5 

Общая  линейная  модель  ► 

11 

67 

1  .30 

2.30 

4 

59 

6 

Обобщенные  линейные  модели  ► 

16 

62 

1  .95 

2.70 

5 

68 

5 

Смешанные  модели  ► 

16 

59 

2,00 

2,60 

6 

89  1 

Корреляции  ► 

17 

73 

1  .82 

2.40 

7 

170  1 

Регрессия  ► 

16 

79 

1  .57 

2.52 

8 

188  1 

Логлинейный  ► 

11 

104 

1  .48 

2.50 

9 

207 

6 

Нейронные  сети  ► 

_ _ 

_ лд 

1—50 — 

2.50 

10 

216 

5 

Классификация  ► 

Двух  этапный  кластерный  анализ 

2.23 

11 

217 

5 

Снижение  размерности  ► 

ЕЗН  Кластеризация  К-средними  ... 
гёІІЛ  Иерархическа я  кластеризация... 

2.64 

12 

220 

1 

Шкалирование  ► 

2.40 

13 

238 

5 

Непараметрические  критерии  ► 

Дере  в  ь  я  классификации . . . 

2,60 

14 

249 

5 

Прогнозирование  ► 

^  Дискриминантный  анализ... 

2.65 

15 

273 

2 

Выживаемость  ► 

..  „ 

2,50 

16 

287  1 

Множественные  ответы  ► 

2^1  _  М 

ЕЯэІ  Анализ  пропущенных  значений... 

17 

288 

5 

16 

71 

1  .96 

2.90 

18 

309 

8 

Множественная  импутация  ► 

18 

65 

1  .69 

2.80 

19 

315 

7 

Сложные  выборки  ► 

17 

79 

1  .63 

2,38 

20 

340 

9 

Контроль  качества  ► 

|^*Э  РОС  кривые... 

16 

59 

1  .61 

2.30 

21 

348 

8 

17 

62 

1  .79 

1  .82 

22 

359 

10 

4  1  2 

17 

73 

1  /4 

2.67 

Рисунок  1.4-  Выбор  метода  анализа  данных 
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После  открытия  диалогового  окна  «Дискриминантный  анализ» 
(рис.  1.5)  единственную  зависимую  переменную  вп  поместили  в  поле, 
предназначенное  для  группирующих  переменных  (рис.  1.6). 


Рисунок  1.5  -  Диалоговое  окно  «Дискриминантный  анализ» 


Нажатием  на  кнопке  «Задать  диапазон»  ввели  минимальное  и 
максимальное  значение  этой  переменной:  1  и  3  соответственно,  так  как 
данная  переменная  имеет  три  уровня  (рис.  1.7).  Затем  нажимали  на  кнопку 
«Продолжить». 


Рисунок  1.6  -  Диалоговое  окно  «Дискриминантный  анализ». 

Выбор  переменных 
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искриминантныи  анализ:  задать  диап 


Минимум:  1 1 


Максимум:  [3 


Продолжить 


Отмена 


Справка 


Рисунок  1.7  —  Диалоговое  окно  «Дискриминантный  анализ: 

Задать  диапазон» 

Всем  остальным  переменным  присвоили  статус  независимых 
переменных,  для  этого  перенесли  их  в  поле  «Независимые».  Далее 
выбрали  метод  «Шаговый  отбор»  (рис.  1.8),  который  рекомендуется  при 
наличии  большого  количества  независимых  переменных. 


Рисунок  1.8  -  Диалоговое  окно  «Дискриминантный  анализ» 

После  нажатия  на  кнопку  «Статистики...»  активировали  опции: 
«Средние»,  «Однофакторный  дисперсионный  анализ»,  «Фишера»  и  другие 
(рис.  1.9). 
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При  выборе  «Среднее»  мы  получили  средние  значения  и 
стандартные  отклонения  для  каждой  переменной  каждой  группы  и  всей 
выборки  в  целом. 


Дискриминантный  анализ:  Статистики 


Описательные - 

0  Средние 

0  Однофакторный  дисперсионный  анализ 
0  М  Бокса 


Коэффициенты  функции 
0  Фишера 


0  [Нестандартизованные 


Продолжить 


Матрицы - 

0  Внутригрупповая  корреляция 
0  Внуіригрупповая  ковариация 
0  Групповые  ковариации 
0  Обш^я  ковариация 


Отмена 


Справка 


Рисунок  1.9  -  Диалоговое  окно  «Дискриминантный  анализ:  Статистики» 


Среди  ограничений,  касающихся  статистических  свойств 
дискриминантных  переменных,  есть  равенство  ковариационных  матриц 
для  различных  классов.  Это  предположение  проверили  с  помощью  М 
критерия  Бокса  [12]: 


8 

М=(п-8)\о%\А\-^ (пк  - 1)  1о§ | К 

к= 1 


где  А  -  объединенная  внутригрупповая  матрица;  К  -  ковариационная 
матрица  для  к-й  группы;  §  -  число  групп;  п  -  общее  число  наблюдений  по 
всем  группам;  пк  -  число  наблюдений  в  к-й  группе.  Алгоритм  получения 
данных  матриц  приведен  ниже.  Результат  проверки  равенства 
ковариационных  матриц  для  различных  классов  с  помощью  М  критерия 
Бокса  приведен  на  рис.  1.10 

Если  значимость  превышает  0,05,  данные  характеризуются 
многомерной  нормальностью.  Если  значимость  равна  или  меньше  0,05, 


15 


многомерная  нормальность  считается  нарушенной,  однако  это  не  означает, 
что  дискриминантный  анализ  не  применим  к  данным;  более  того, 
включение  в  уравнение  переменной,  чье  распределение  отличается  от 
нормального,  иногда  может  улучшить  результаты  анализа.  Отсутствие 
многомерной  нормальности  указывает  лишь  на  то,  что  не  все 
используемые  переменные  имеют  нормальное  распределение  [13]. 


Рисунок  1.10  -  Результаты  теста  Бокса 


Этап  5.  Определение  коэффициентов  дискриминантных  функций. 
Ядром  дискриминантного  анализа  является  построение  так 
называемой  дискриминантной  функции: 

Щт  =  Ъ\  'Х\кт  +Ъ2-Х2Ьп+...  +  Ър  • Хркпг  +Ъ0  ,  (1.1) 
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где  ОР/ап  -  значение  дискриминантной  функции  для  т- го  пациента  в 

к-й  группе;  хікт  -  значение  дискриминантной  переменной  і  для  т-го 

пациента  в  к-й  группе;  р  -  число  дискриминантных  переменных,  Ц  ~ 
коэффициенты,  которые  предстоит  оценить  с  помощью  дискриминантного 
анализа. 

Целью  являются  определение  таких  коэффициентов,  чтобы  по 
значениям  дискриминантной  функции  можно  было  с  максимальной 
чёткостью  провести  разделение  на  группы. 

В  дискриминантном  анализе  коэффициенты  (или  веса)  Ьо,  Ъь  Ь2 ,  ..., 
Ър  определяют  таким  образом,  чтобы  группы  максимально  отличались 
значениями  дискриминантной  функции.  Это  происходит  тогда,  когда 
отношение  межгрупповой  суммы  квадратов  к  внутригрупповой  сумме 
квадратов  для  дискриминантных  показателей  максимально. 

Мы  выбрали  пошаговый  (віершве)  дискриминантный  анализ,  при 
котором  переменные  вводятся  последовательно,  исходя  из  их  способности 
различать  (дискриминировать)  группы.  При  пошаговом  анализе  «с 
включением»  на  каждом  шаге  просматриваются  все  переменные,  и 
находится  та  из  них,  которая  вносит  наибольший  вклад  в  различие  между 
совокупностями.  Эта  переменная  должна  быть  включена  в  модель  на 
данном  шаге,  и  происходит  переход  к  следующему  шагу. 

Щелчком  на  кнопке  «Метод»  в  окне  «Дискриминантный  анализ» 
(рис.  1.8)  мы  вошли  в  диалоговое  окно  «Дискриминантный  анализ: 
Шаговый  отбор»,  показанное  на  рис  1.11.  Кнопка  «Метод»  доступна  лишь 
в  том  случае,  когда  установлен  переключатель  «Шаговый  отбор». 

В  этом  окне  при  помощи  группы  переключателей  «Метод»  можно 
выбирать  один  из  пошаговых  методов  составления  дискриминантного 
уравнения,  в  области  «Критерии»  задать  критерии  для  включения  в 
дискриминантное  уравнение  и  исключения  из  него  предикторов,  а  также 
включить  в  окно  вывода  желаемые  величины  при  помощи  флажков  в 
группе  «Вывести  на  дисплей».  Метод  Лямбда  Уилкса  основан  на 
минимизации  коэффициента  Уилкса  после  включения  в  уравнение 
регрессии  каждого  нового  предиктора. 
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Рисунок  1.11  —  «Дискриминантный  анализ:  Шаговый  отбор» 

По  умолчанию  критериями  для  включения  и  исключения 
предикторов  являются  пороговые  значения  Р-критерия.  Как  правило 
Р=3,84  соответствует  величина  уровня  значимости  р,  приблизительно 
равная  0,05,  а  значению  Р=2,71  -  величина  р,  приблизительно  равная 
ОД  [14]. 

Рассмотрим  основные  этапы  получения  коэффициентов 
«канонической  дискриминантной  функции. 

Для  получения  коэффициентов  Д  канонической  дискриминантной 
функции  нужен  статистический  критерий  различения  групп.  Проверка 
качества  дискриминации  (различия)  основана  на  сравнении  средних 
дискриминантной  функции  для  исследуемых  групп.  Эти  средние  играют 
настолько  важную  роль  в  дискриминантном  анализе,  что  получили  свое 
название  -  центроиды  {сепігоісіз).  Центроидов  столько,  сколько  групп,  т.е. 
один  центроид  для  каждой  группы.  Очевидно,  что  классификация 
переменных  будет  осуществляться  тем  лучше,  чем  меньше  рассеяние  точек 
относительно  центроида  внутри  группы  и  чем  больше  расстояние  между 
центроидами  групп.  Разумеется,  что  большая  внутригрупповая  вариация 
нежелательна,  так  как  в  этом  случае  любое  заданное  расстояние  между 
двумя  средними  тем  менее  значимо  в  статистическом  смысле,  чем  больше 
вариация  распределений,  соответствующих  этим  средним.  Один  из 
методов  поиска  наилучшей  дискриминации  данных  заключается  в 
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нахождении  такой  канонической  дискриминантной  функции  БР,  которая 
бы  максимизировала  отношение  межгрупповой  вариации  к 
внутригрупповой: 


Я  =  5/Е, 


(1.2) 


где  5  -  межгрупповая  и  Е  -  внутригрупповая  матрицы  рассеяния 
наблюдаемых  переменных  от  средних. 

Рассмотрим  максимизацию  отношения  (1.2)  для  произвольного  числа 
классов.  Введем  следующие  обозначения: 

хі  -  среднее  значение  переменной  і  по  всем  классам  (общее 
среднее); 

хікт  -  значение  переменной  і  для  т-го  наблюдения  в  к-й  группе; 

Хік  -  среднее  значение  переменной  і  в  к-й  группе. 

Рассмотрим  задачу  максимизации  отношения  (1.2)  когда  имеются  § 
групп.  Оценим  сначала  информацию,  характеризующую  степень  различия 
между  объектами  по  всему  пространству  точек,  определяемому 
переменными  групп.  Для  этого  вычислим  матрицу  рассеяния  О,  которая 
равна  сумме  квадратов  отклонений  и  попарных  произведений  наблюдений 

от  общих  средних  хі,і  =  1,...,р  По  каждой  переменной.  Элементы 
матрицы  О  определяются  выражением: 


Запишем  это  выражение  в  матричной  форме.  Обозначим  р-мерную 
случайную  векторную  переменную  к-й  группы  следующим  образом: 


Тогда  объединенная  /^-мерная  случайная  векторная  переменная  всех 
групп  будет  иметь  вид: 
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X  = 


ХхХ2..Х 


Общее  среднее  этой  р-мерной  случайной  векторной  переменной 
будет  равно  вектору  средних  отдельных  признаков: 


X  = 


X  ^  X  2  •  •  •Х  ^ 


Матрица  рассеяния  от  среднего  при  этом  запишется  в  виде 

0  =  '%(Хк-х)(Хк-х)т  _ 

к= 1 


Если  использовать  векторную  переменную  объединенных 
переменных  X,  то  матрица  О  определится  по  формуле: 

О  =  (X -х)(Х -х)т . 

Матрица  О  содержит  полную  информацию  о  распределении  точек  по 
пространству  переменных.  Диагональные  элементы  представляют  собой 
сумму  квадратов  отклонений  от  общего  среднего  и  показывают,  как  ведут 
себя  наблюдения  по  отдельно  взятой  переменной.  Вне  диагональные 
элементы  равны  сумме  произведений  отклонений  по  одной  переменной  на 
отклонения  по  другой. 

Если  разделить  матрицу  О  на  (п  -  1),  то  получим  общую 
ковариационную  матрицу  Н  (рис.  А.1,  Приложение  А).  В  пакете  8Р88 
общая  ковариационная  матрица  получается  путем  установки  опции 
«Общая  ковариация»  в  диалоговом  окне  «Дискриминантный  анализ: 
Статистики»  (рис.  1.9).  Для  проверки  условия  линейной  независимости 
переменных  полезно  рассмотреть  вместо  О  нормированную 
корреляционную  матрицу. 

Для  измерения  степени  разброса  объектов  внутри  групп  рассмотрим 
матрицу  Е,  которая  отличается  от  О  только  тем,  что  ее  элементы 
определяются  векторами  средних  для  отдельных  групп,  а  не  вектором 
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средних  для  общих  данных.  Элементы  внутригруппового  рассеяния 
определятся  выражением: 


8  пк 


^  ^  '  (-^ікт  Д  ) 


к—  1  /72=1 


(1.3) 


Выражения  в  скобках  являются  отклонениями  значений  переменных 
от  общего  среднего.  Если  і  =  /,  то  сомножители  равны,  и  получается 
среднеквадратичное  отклонение.  Таким  образом,  диагональные  элементы 
представляют  собой  сумму  квадратов  отклонений  от  общего  среднего.  Они 
показывают  как  ведут  себя  наблюдения  по  отдельной  переменной.  При 
і  ф  ]  получаем  сумму  произведений  отклонения  по  одной  переменной  на 
отклонение  по  другой.  В  этом  состоит  один  из  способов  измерения 
корреляций  (ковариаций)  между  двумя  переменными,  так  как  он 
показывает,  насколько  хороню  большое  отклонение  по  одной  переменной 
согласуется  с  большим  отклонением  по  другой.  Рассматривая  целиком  всю 
матрицу,  мы  имеем  полную  информацию  о  распределении  точек  по 
пространству,  определяемому  переменными. 

Запишем  выражение  (1.3)  в  матричной  форме.  Данным  групп  будут 
соответствовать  векторы  средних: 


ч 


Х\кХ2к’”Х  рк 


Тогда  матрица  внутригрупповых  вариаций  запишется  в  виде: 

Е  =  ^(Хк-7к)(Хк-7к)Т 

к= 1 


Если  разделить  каждый  элемент  матрицы  Е  на  (п  -  §),  то  получится 
внутригрупповая  ковариационная  матрица  А,  она  является  взвешенным 
средним  ковариационных  матриц  отдельных  классов.  Данная  матрица  А 
(рис.  А.2,  Приложение  А)  была  получена  путем  установки  опции 
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«Внутригрупповая  ковариация»  в  диалоговом  окне  «Дискриминантный 
анализ:  Статистики»  (рис.  1.9). 

Степень  зависимости  двух  переменных  можно  выяснить,  исследуя  их 
корреляцию.  Для  этого  мы  в  окне  «Дискриминантный  анализ:  Статистики» 
(рис.  1.9)  отметили  флажком  «Внутригрупповая  корреляция». 

Коэффициент  корреляции  принимает  значения  от  -1  до  +1. 
Преобразование  матрицы  Е  в  матрицу  коэффициентов  корреляции 
(объединенную  внутригрупповую  матрицу)  осуществляется  делением 
каждого  элемента  на  квадратный  корень  произведения  двух 

соответствующих  диагональных  элементов,  а  именно 


на 


к=  1  т= 1 


О  -/С  8  пк 

ЕЕ  (А ікт  )  2ЕЕ  ^  ]кт  ^  А  ) 

к=  1  т= 1 


.  На  рис.  А.З  (Приложение  А) 


представлена  объединенная  внутригрупповая  матрица  (),  состоящая  из 
соответствующих  средних  значений  корреляционных  матриц  для  уровней 
зависимой  переменной. 

Для  получения  групповой  ковариационной  матрицы  сначала 
определяются  элементы  8  у : 


пк 

§і]  (  ^ікіп  )(  ^/кт  Х.)к ) 

т= 1 


Если  элементы  матрицы  8 у  разделить  на  (щ  -  1),  то  получится 

групповая  ковариационная  матрица  К.  На  рис.  А.4  (Приложение  А) 
представлена  полученная  в  8Р88  групповая  ковариационная  матрица  К, 
представляющая  собой  отдельные  ковариационные  матрицы  для  каждого 
уровня  зависимой  переменной.  Для  получения  данной  матрицы  в  окне 
«Дискриминантный  анализ:  Статистики»  (рис.  1.9)  предварительно 
отметили  флажком  «Групповые  ковариации». 

Когда  центроиды  различных  групп  совпадают,  то  элементы  матриц 
О  ж  Е  будут  равны.  Если  же  центроиды  групп  различные,  то  разница 
5=0-Е  будет  определять  межгрупповую  сумму  квадратов  отклонений  и 
попарных  произведений.  Если  расположение  групп  в  пространстве 
различается  (т.е.  их  центроиды  не  совпадают),  то  степень  разброса 
наблюдений  внутри  групп  будет  меньше  межгруппового  разброса. 
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Отметим,  что  элементы  матрицы  51  можно  вычислить  и  по  данным 
средних: 


л- 


У 


н  ч  _ 

к= 1  т=1 


Хі)(х]к-Хі),  /,7  =  1, 


Р- 


Матрицы  Е  и  5  содержат  всю  основную  информацию  о  зависимости 
внутри  групп  и  между  группами.  Для  лучшего  разделения  наблюдений  на 
группы  нужно  подобрать  коэффициенты  дискриминантной  функции  из 
условия  максимизации  отношения  межгрупповой  матрицы  рассеяния  к 
внутригрупповой  матрице  рассеяния  при  условии  ортогональности 
дискриминантных  плоскостей.  Тогда  нахождение  коэффициентов 
дискриминантных  функций  сводится  к  решению  задачи  о  собственных 
значениях  и  векторах.  Это  утверждение  можно  сформулировать  так:  если 
спроектировать  §  групп  /7-мерных  выборок  на  {§  -1)  пространство, 

порожденное  собственными  векторами  ( У\к  )’  ^  —  1» •••»<?  —  ^  то 

отношение  (1.2)  будет  максимальным,  т.  е.  рассеивание  между  группами 
будет  максимальным  при  заданном  внутригрупповом  рассеивании.  Если 
бы  мы  захотели  спроектировать  §  выборок  на  прямую,  при  условии 
максимизации  наибольшего  рассеивания  между  группами,  то  следовало  бы 

использовать  собственный  вектор  (ѵц»--»ѵрі),  соответствующий 

максимальному  собственному  числу  .  При  этом  дискриминантные 
функции  можно  получать:  по  нестандартизованным  и  стандартизованным 
коэффициентам. 

Пусть  2  і>Я  2---2р  и  ѵ\кт-іѵрк  соответственно  собственные 

значения  и  векторы.  Тогда  условие  (1.2)  в  терминах  собственных  чисел  и 
векторов  запишется  в  виде: 


Л-  = 


к 


что  влечет  следующее  равенство,  представленное  в  матричной  форме: 
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(1.4) 


(5-ЛіЕ)ѵі=0,ѵ1тЕѵі=3^ 


где  ду  -  символ  Кронекера  -  функция  двух  целых  переменных, 
которая  равна  1,  если  они  равны,  и  0  в  противном  случае  [15].  Таким 
образом,  решение  уравнения  I  5  -  АЕ 1=  0  позволяет  нам  определить 
компоненты  собственных  векторов,  соответствующих  дискриминантным 
функциям.  Решение  системы  уравнений  (1.4)  можно  получить  путем 
использования  разложения  Холецкого  матрицы  Е'1  и  решения  задачи  о 
собственных  значениях  [16  -  17]. 

Каждое  решение,  которое  имеет  свое  собственное  значение  Я,  и 
собственный  вектор  ѵ„  соответствует  одной  дискриминантной  функции. 
Компоненты  собственного  вектора  ѵг  можно  использовать  в  качестве 
коэффициентов  дискриминантной  функции.  Однако  при  таком  подходе 
начало  координат  не  будет  совпадать  с  главным  центроидом.  Для  того 
чтобы  начало  координат  совпало  с  главным  центроидом,  нужно 
нормировать  компоненты  собственного  вектора: 


Р 


Кк  =-^ь<кхі- 


(1.5) 


/= 1 


Нормированные  коэффициенты  (1-5)  получены  по 


нестандартизованным  исходным  данным,  поэтому  они  называются 
нестандартизованными.  Нормированные  коэффициенты  приводят  к  таким 
дискриминантным  значениям,  единицей  измерения  которых  является 
стандартное  квадратичное  отклонение.  При  таком  подходе  каждая  ось  в 
преобразованном  пространстве  сжимается  или  растягивается  таким 
образом,  что  соответствующее  дискриминантное  значение  для  данного 
объекта  представляет  собой  число  стандартных  отклонений  точки  от 
главного  центроида. 

Расчет  методом  дискриминантных  функций  значений 
диагностических  коэффициентов  позволил  выделить  десять  значимых  для 
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определения  хроническом  сердечной  недостаточности  показателен  со 
своими  нестандартизованными  коэффициентами  (табл.  1.1). 


Таблица  1.1 

Канонические  коэффициенты  дискриминантных  функций  для  значимых 
показателей  при  определении  состояния  сердечно-сосудистой  системы 


Показатели 

Единицы  измерения 

Канонические 
коэффициенты,  Ьі 

1-я 

функция 

2-я 

функция 

*і 

пг/мл 

-  0,019 

0,080 

*2 

пг/мл 

-  0,016 

-  0,059 

пг/мл 

0,001 

-  0,010 

*4 

%  клеток  с  маркером 
апоптоза 

0,053 

0,130 

нмоль/мл 

-  0,003 

-  0,038 

нг/мл/ч 

0,156 

0,296 

х7 

ед/мл 

-0,165 

0,270 

нмоль/сут 

0,001 

-0,011 

*9 

% 

0,120 

-  0,014 

х10 

мл/м2 

0,020 

0,009 

Константа 

— 

-7,188 

-  1,122 

Таким  образом,  развитие  ХСН  у  подростков  может  быть  описано 
следующими  дискриминантными  функциями  [18  -  19]: 


БГ1 
+0,156*  *6 


=  -  0,019* *і  -  0,016*  *2  +  0,001**3  +  0,053**4-  0,003* *5 
-0,165**7  +0,001**8  +0,120**9  +  0,020*  *ш  -7,188, 


+ 
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где  хг  -  интерлейкин- 1  (пг/мл);  х2  -  интерлейкин-6  (пг/мл);  т3  - 
фактор  некроза  опухоли  ФНО-а  (пг/мл);  х4  -  СО-95  (%  клеток,  несущих 
маркер  апоптоза);  *5  -  циклический  3,5-аденозинмонофосфат  (цАМФ, 
нмоль/мл);  х6  -  ренин  (нг/мл/час);  х7  -  карбонилированные  белки  (Ед/мл); 
х8—  норадреналин  (нмоль/сут);  х9-  фракция  выброса  левого  желудочка 
после  физической  нагрузки  (%);  тІ0  —  индекс  конечного  диастолического 
объема  левого  желудочка  (мл/м2). 

Стандартизованные  коэффициенты  можно  получить 
преобразованием  нестандартизованных  коэффициентов  к 

стандартизованной  форме  [20] : 


где  е9  -  сумма  внутригрупповых  квадратов  і- й  переменной, 
определяемой  по  формуле  3.  Стандартизованные  коэффициенты 
применяют  для  уменьшения  размерности  исходного  признакового 
пространства  переменных.  Если  абсолютная  величина  коэффициента  для 
данной  переменной  для  всех  дискриминантных  функций  мала,  то  эту 
переменную  можно  исключить,  тем  самым  сократив  число  переменных. 
Полученные  стандартизованные  коэффициенты  канонической 
дискриминантной  функции  приведены  на  рис.  Е12 

Этап  6.  Определение  коэффициентов  канонической  корреляции. 

Оценка  степени  удачного  распределения  на  группы,  полезность 
дискриминантных  функций  и  количество  функций,  имеющих  реальное 
содержание  при  определении  различий  между  группами,  оценены  с 
помощью  коэффициентов  канонической  корреляции  (табл.  Е2). 
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Рисунок  1.12  -  Стандартизованные  коэффициенты  канонической 

дискриминантной  функции 

Полученные  собственные  значения  приведены  в  порядке  убывания 
их  величин.  Величина  собственного  значения  связана  с 
дискриминирующими  возможностями  функции:  чем  больше  собственное 
значение,  тем  лучше  различие.  Таким  образом,  первая  функция  обладает 
наибольшими  дискриминирующими  возможностями,  вторая 
обеспечивает  максимальное  различие  после  первой.  Фактические  числа, 
представляющие  собственные  значения,  трудно  интерпретировать 
непосредственно. 

Чтобы  облегчить  сравнение,  собственным  значением  приписывают 
процент  дисперсии.  Для  этого  сначала  суммируют  все  собственные 
значения  (13,010+1,264=14,274),  затем  делят  каждое  собственное  значение 
на  общую  сумму  (13,010/14,274=0,911  и  1,264/14,274=0,089).  Так  первая 
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функция  содержит  91,1%  общих  дискриминантных  возможностей,  а 
вторая  -  8,9%. 


Таблица  1.2 


Характеристика  дискриминантных  возможностей  функций 


Функция 

Собственное 

значение  Л 

%  дисперсии 

Общий  % 

Каноническая 

корреляция 

гі 

1 

13,010 

91,1 

91,1 

0,964 

2 

1,264 

8,9 

100,0 

0,747 

Мы  видим,  что  вторая  функция  малозначима,  но  не  знаем,  можно  ли 
ею  пренебречь,  так  как  нет  правила,  которое  помогло  бы  определить,  как 
велико  должно  быть  относительное  процентное  содержание,  чтобы 
функция  представляла  для  нас  интерес.  Относительное  процентное 
содержание  показывает  только  то,  насколько  функция  слабее  по  сравнению 
с  другими  и  ничего  не  добавляет  к  определению  различий  между 
группами. 

Для  того  чтобы  определить,  представляет  ли  интерес,  например, 
вторая  функция  необходимо  определить  коэффициент  канонической 
корреляции. 

В  дискриминантном  анализе  различают  две  задачи:  интерпретации  и 
классификации.  Задача  интерпретации  -  определение  количества, 
значимости  дискриминантных  функций  и  их  значений  для  объяснения 
различий  между  классами.  Задача  классификации  -  определение  класса,  к 
которому  принадлежит  новый  объект.  В  некоторых  случаях  нет 
необходимости  решать  задачу  интерпретации.  В  этой  ситуации 
используют  простые  классифицирующие  функции,  называемые 
дискриминантными,  а  при  решении  задачи  интерпретации  -  канонические 
дискриминантные  функции.  Каноническая  корреляция  показывает  степень 
зависимости  между  дискриминантными  функциями  и  классами.  Нулевое 
значение  говорит  об  отсутствии  связи,  а  большие  числа  (всегда 
положительные)  означают  большую  степень  зависимости  (максимальное 
значение  равно  1). 
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Следовательно,  оценку  реальной  полезности  дискриминантной 
функции  получили,  определив  коэффициент  канонической  корреляции  г, 
[21]: 


П 


(1.6) 


где  і  —  номер  соответствующей  дискриминантной  функции; 

1і  -  собственное  значение. 

Например,  коэффициент  канонической  корреляции  для  первой 
дискриминантной  функции  определяется  как: 


г\ 


13,010 
1  +  13,010 


=  0,964 


Значение  коэффициента  канонической  корреляции,  первой  функции 
(0,964)  позволяет  сделать  вывод  о  существовании  достаточно  высокой 
связи  между  ХСН  и  значением  первой  дискриминантной  функции,  что  и 
предполагалось  процентом  дисперсии  этой  функции.  Также  выявлена 
высокая  связь  между  случаями  развития  ХСН  (0,747)  и  значение  второй 
дискриминантной  функции.  Качественная  оценка  силы  связи  гХу  величин  X 
и  У  была  сделана  на  основании  шкалы  Чеддока  (табл.  1.3)  [22]. 


Таблица  1.3 


Таблица  анализа  силы  связи  между  переменными 


Значение  г, 

Интерпретация 

от  0  до  0,3 

Очень  слабая 

от  0,3  до  0,5 

Слабая 

от  0,  5  до  0,7 

Средняя 

от  0,7  до  0,  9 

Высокая 

от  0,9  до  1 

Очень  высокая 

Этап  7.  Определение  значимости  дискриминантных  функций 
методом  ^-статистики  Уилкса. 

Оценка  значимости  дискриминантных  функций  была  проверена 
методом  ^-статистики  Уилкса  (табл.  1.4). 
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^-статистика  Уилкса  -  это  мера  различий  между  классами  по 
нескольким  дискриминантным  переменным,  ^-статистика  Уилкса 
принимает  значения,  обратные  привычным,  и  оценивает  остаточную 
дискриминантную  способность  Чем  ближе  значение  X  к  0,  тем  лучше 
различие  классов,  а  чем  ближе  к  1,  тем  различие  хуже  (классы  совпадают). 

Определяется  ^-статистика  Уилкса  согласно  формуле: 


«  1 

я*  =  ГТ  — 

«Ый+ѵ 


(1.7) 


где  к  -  число  вычисленных  функций; 
Хі  -  собственное  значение. 


Таблица  1.4 

Определение  значимости  дискриминантных  функций  методом 

^-статистики  Уилкса 


Тест 

функций 

Лямбда 

;* 

Уилкса  лі 

X? 

Степени 

свободы,  СІ/і 

Значимость 

р, 

1 

0,032 

523,746 

20 

0,0001 

2 

0,442 

123,819 

9 

0,0001 

Например,  ^-статистика  Уилкса  для  первой  функции  определяется 

как: 


*  1  1 

А,  = - =  0,032 

1  +  13,010  1  +  1,264 

Для  второй  дискриминантной  функции  ^-статистика  Уилкса 
рассчитывается  как: 


- =  0,442 

1  +  1,264 
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Поскольку  X  является  «обратной»  мерой,  этот  результат  означает,  что 
все  используемые  переменные  эффективно  участвуют  в  различении 
классов. 

Из  табл.  1.4  видно,  что  вычисление  первой  функции  уменьшает 

л  * 

количество  оставшейся  информации,  и  величина  л-  становится  равной 
0,442  (остается  небольшой). 

2 

Этап  8.  Расчет  критериального  значения  Х\  ,  которое  сравнивается  с 
критическим,  с  заданным  числом  степеней  свободы  сі/і  =  (р-(к-1))(§-(к-1)-1. 

Число  степеней  свободы  зависит  от  количества  дискриминантных 
переменных  р  и  количества  классов  §: 


'  Р  +  8Л 

I  2  ) 


(1.8) 


Так,  например,  число  степеней  свободы  для  первой  и  второй 
дискриминантных  функции  определяется  как: 

4Гі'=  (10-(1-1))(3-(1-1)-1)=20, 

#2  =  (10-(2-1))(3-(2-1)-1)=9. 


Критериальное  значение  Хі  ,  например,  для  второй  дискриминантной 
функции  определяется  как: 


^10  +  3Л 

к  2  у 


ІиО, 442  «123,819. 


В  нашем  случае,  уровень  значимости  Р;<0,0001,  то  есть  различие 
между  группами  значимое. 

Этап  10.  Расчет  квадрата  расстояния  Махалонобиса  до  среднего 
значения  группы  значений  дискриминантной  функции. 

Рассчитанные  значения  дискриминантных  функций  рассматриваются 
как  точки  в  некотором  пространстве.  Для  каждой  группы  можно  рассчитать 
центр  группирования  (центроид)  [23]: 
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с 


пк 

X  ОРЬп 
т= 1 


В  табл.  5  приведены  полученные  в  8Р88  значения  центроидов  для 
каждой  группы. 


Таблица  1.5 

Значения  центроидов  для  каждой  группы 


Группы 

Функция 

1 

2 

Без  сердечной  недостаточности 

2,099 

2,374 

1  стадия  сердечной  недостаточности 

3,409 

-,815 

Контрольная  группа 

-4Д53 

-,166 

Для  каждого  объекта  рассчитывают  расстояние  от  него  до  каждой 
точки  группирования.  Для  этого  используют  квадрат  расстояния 
Махалонобиса  (расстояние  от  объекта,  который  необходимо 
классифицировать,  до  центра  класса)  [24]. 

Для  того  чтобы  в  8Р88  получить  информацию  о  квадрате  расстояния 
Махалонобиса  в  диалоговом  окне  «Дискриминантный  анализ».  Выбор 
метода  (рис.  1.6)  выбрали  кнопку  «Классифицировать»  и  активировали 
опцию  «Поточечные  результаты»  (рис.  1.13). 

В  полученной  в  8Р88  после  включения  опции  «Поточечные 
результаты»  (рис.  1.13)  таблице  А.2  (Приложение  А)  в  столбце  «Квадрат 
расстояния  Махалонобиса  до  центра»  приводится  квадрат  расстояния 
Махалонобиса  до  среднего  значения  группы  значений  дискриминантной 
функции.  В  столбце  «Значения  дискриминантных  функций»  приводятся 
величины,  получаемые  при  подстановке  значений  переменных  пациента  в 
уравнение  дискриминантных  функций. 
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Рисунок  1.13  -  Диалоговое  окно  «Дискриминантный  анализ: 

Классификация» 

Этап  11.  Определение  группы,  к  которой  принадлежит  пациент. 

Относя  объект  к  ближайшему  классу  в  соответствии  со  значением 
И2,  мы  неявно  приписываем  его  к  тому  классу,  для  которого  он  имеет 
наибольшую  вероятность  принадлежности. 

Группа,  к  которой  фактически  принадлежит  пациент,  отображается  в 
колонке  «Фактическая  группа»  (табл.  А.2,  Приложение  А).  В  следующих 
колонках  содержится  информация  о  прогнозе  принадлежности  к  группе, 
сделанном  на  основании  значений  дискриминантных  функций. 

Р(0>с1  I  Сг=§)  -  это  вероятность  того,  что  пациент,  принадлежащий  к 
прогнозируемой  группе,  действительно  имеет  значения  параметров, 
соответствующие  дискриминантной  функции  или  некоторые  другие 
крайние  значения. 

Р(0=§І0=с1)  -  мера  принадлежности  к  одной  из  двух  групп.  Это 
вероятность  того,  что  некоторый  случай  принадлежит  к  прогнозированной 
группе,  которая  рассчитывается  на  основе  подстановки  в 
дискриминантную  функцию  значений  набора  переменных, 
соответствующих  данному  случаю.  Вероятность  того,  что  данный  случай 
принадлежит  к  другой  группе  получается  вычитанием  меры 
принадлежности  из  1 . 
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Этап  12.  Построение  территориальной  карты  на  основе  значений 
дискриминантных  функций. 

Запрос  на  вывод  территориальной  карты  осуществляется  с  помощью 
опции  «Территориальная  карта»  в  диалоговом  окне  «Дискриминантный 
анализ:  Классификация»  (рис.  1.13). 

На  основании  значений  дискриминантный  функций  была  построена 
территориальная  карта  распределения  больных  в  зависимости  от  диагноза 
(рис.  1.14),  где  1  -  группа,  у  которой  отсутствует  сердечная 
недостаточность;  2  -  группа,  у  которой  начальная  стадия  сердечной 
недостаточности;  3  -  группа  здоровых  подростков;  *  -  центр  группы. 
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Рисунок  1.14  -  Территориальная  карта  распределения  больных  с 
хронической  сердечной  недостаточностью  и  без  неё 

Этап  13.  Построение  графика  распределения  пациентов  в 
зависимости  от  диагноза. 

В  диалоговом  окне  «Дискриминантный  анализ:  Классификация» 
(рис.  1.13)  сделали  запрос  на  вывод  графика  распределения,  установив 
флажок  в  поле  «Объединенные  группы». 
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На  рис.  1.15  представлен  график  распределения  пациентов 
зависимости  от  диагноза,  где  по  оси  X  располагаются  значения 
дискриминантной  функции  1  (2)Г7),  а  по  оси  У  -  значения 

дискриминантной  функции  2  (ОУ2). 


Сердечная 

недостаточность 

О  Без  сердечной  недостаточности 
О  Первая  стадия  сердечной  недостаточности 
О  Контрольная  группа 
#  Центроид  группы 


Рисунок  1.15  -  График  распределения  пациентов  с  различными 
состояниями  сердечно-сосудистой  системы 


Этап  14.  Определение  1 -статистики  ошибок  и  доли  правильной 
диагностики.  Стандартизованной  мерой  эффективности  для  любого 
количества  классов  является  т -статистика  ошибок  [25  -  27]: 


8 

пс~Т-  рл 

- И - ,  (1.9) 

О 

п~ТРкпк 

к= 1 


где  пс  -  число  правильно  классифицированных  объектов,  а  рк  - 
априорная  вероятность  принадлежности  к  классу. 


35 


8 

Выражение  ркпк  представляет  собой  число  объектов,  которые 

к= 1 

будут  правильно  предсказаны  при  случайной  классификации  их  по 
группам  пропорционально  априорным  вероятностям.  Если  все  группы 
считаются  равноправными,  то  априорные  вероятности  полагаются 
равными  единице,  деленной  на  число  классов.  Максимальное  значение  1  - 
статистики  равно  1  и  оно  достигается  в  случае  безошибочного 
предсказания.  Нулевое  значение  указывает  на  неэффективность 
процедуры,  1  -статистика  может  принимать  и  отрицательные  значения,  что 
свидетельствует  о  плохом  различении  или  вырожденном  случае. 

Поскольку  пс  должно  быть  целым  числом,  числитель  может  стать 
отрицательным  чисто  случайно,  когда  нет  различий  между  группами. 

Вывод  в  8Р88  таблицы,  представленной  на  рис  1.16,  осуществляется 
путем  установки  опции  «Итоговая  таблица»  в  диалоговом  окне 
«Дискриминантный  анализ:  Классификация»  (рис.  1.14). 


Рисунок  1.16  -  Результаты  классификации  ХСН  в  8Р88 
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К  группе  1  (подростки  с  заболеваниями  сердца  без  ХСН)  относится 
27  наблюдений.  24  наблюдений  (88,  9  %)  были  определены  корректно,  а  3 
(11,1  %)  по  ошибке  отнесены  к  группе  2.  В  группе  2  (подростки  с 
начальной  стадией  ХСН)  насчитывается  65  наблюдений.  Все  они  были 
определены  корректно.  В  группе  3  (контроль)  насчитывается  67 
наблюдения.  Все  они  также  были  определены  корректно  (рис.  1.16). 

Анализ  результатов,  полученных  с  помощью  8Р88,  показал,  что 
прогноз  для  исследуемых  групп  сделан  верно  (корректно 
классифицированы  в  общей  сложности  98,1%  всех  случав  заболеваний 
сердца). 

Для  определения  т-статистики  ошибок  и  определения  априорных 
вероятностей  рк  для  каждой  группы  пациентов  установили  в  поле 
диалогового  окна  «Дискриминантный  анализ:  Классификация»  (рис.  1.13) 
флажок  «Вычислить  по  размерам  групп»  (рис.  1.17). 


Рисунок  1.17  -  Априорные  вероятности  для  всех  групп 
%ркпк  =(0,170- 27)  +  (0,409-65)  +  (0,421 -67)  =  59,382. 

к=  1 


_  156-59,382 
Т_  159-59,382 


0,9698. 
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Т.е.  классификация  с  помощью  дискриминантных  функций  делает  на 
96,9%  ошибок  меньше,  чем  ожидалось  при  случайной  классификации. 

Нередко  как  показатель  качества  алгоритма  диагностики 
(прогностической  «силы»)  используют  долю  правильной  диагностики: 

ц=рік+р2к+р3т|,  (1.10) 


где  к  -  доля  правильной  диагностики  в  1-й  группе; 

X  -  доля  правильной  диагностики  во  2-й  группе; 
т|—  доля  правильной  диагностики  в  3-й  группе; 

рк(к=1,  2,  3;  рі+р2+рз=1)  -  доли  групп  в  объединенной  совокупности. 

р1=27/159=0,170;  р2=65/159=0,409;  р3=67/159=0,421. 

8 

X  Ркпк=  0,170  +  0,409  +  0,421  =  1. 
к= 1 


24 

к  =  —  =  0,889; 
27 


X 


65 

65 


=  1; 


//  =  0,170-0, 889  +  0, 409-1  +  0, 421-1  =  0, 9811. 

Данный  результат  соответствует  полученному  в  8Р88  (рис.  1.16). 

Для  корректного  применения  дискриминантных  функции  при 
решении  поставленной  задачи  использовались  две  выборки  пациентов,  по 
159  человек  каждая:  одна  для  вычисления  дискриминантной  функции 
(анализируемая),  вторая  -  проверочная  (для  проверки  результатов  расчета 
на  основании  первой  выборки).  Т.е  была  проведена  кросс-проверка 
(перекрестная  проверка). 
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В  результате  проведения  клинической  апробации  были  получены 
следующие  результаты  (рис.  1.18). 

К  группе  1  относится  27  наблюдений.  25  наблюдений  (92,6  %)  было 
определено  корректно,  2  ошибочно  отнесены  к  группе  2.  В  группе  2 
насчитывается  65  наблюдений.  65  из  них  были  определены  корректно.  В 
группе  3  насчитывается  67  наблюдений.  Корректно  определено  67. 


Рисунок  1.18  -  Результаты  классификации 

Анализ  полученных  результатов  показал,  что  прогноз  для 
исследуемых  групп  сделан  верно  (корректно  классифицированы  в  общей 
сложности  98,  7  %  всех  случав  заболеваний  сердца). 

Результаты  оценки  показателя  качества  алгоритма  диагностики 
(прогностической  «силы»)  приведены  ниже. 

Рі=27/159=0,170;  р2=65/159=0,409;  р3=67/1 59=0,421. 

X  ркпк  0,170  +  0,409  +  0,421  =  1. 

к= 1 


25 

к- =  —  =  0,926; 
27 
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//  =  0,170-0, 926  +  0, 409-1  +  0, 421-1  =  0, 9874. 

Этап  15.  Формирование  диагностического  заключения. 

Таким  образом,  используя  аппарат  дискриминантного  анализа,  на 
основе  измерения  морфофункциональных  характеристик  сердца  пациента, 
показателей  систем  нейрогуморальной  регуляции, 

иммунновоспалительной  активации  и  свободнорадикальных  процессов 
классифицировали  пациентов,  т.е.  определили  группу,  к  которой 
принадлежит  пациент:  (с  заболеваниями  сердечно-сосудистой  системы  без 
признаков  ХСН;  с  заболеваниями  сердца  и  начальной  стадией  ХСН; 
контрольная  группа,  которую  составили  практически  здоровые  подростки. 
Это  позволило  предсказать  группы  (наличие  сердечной  недостаточности) 
для  новых  объектов  из  той  же  совокупности. 

В  данной  главе  показана  возможность  и  целесообразность 
использования  дискриминантного  анализа  для  идентификации  состояния 
больных  ХСН. 
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ГЛАВА  2.  ОПРЕДЕЛЕНИЕ  ДИСФУНКЦИИ  МИОКАРДА  У 
ПОДРОСТКОВ  С  ИСПОЛЬЗОВАНИЕМ  ПАКЕТА  ПРОГРАММ 

8ТАТТ8ТІСА 


Одним  из  ранних  проявлений  доклинической  стадии  ХСН  является 
систолическая  дисфункция  миокарда.  Это  состояние  у  детей  и  подростков 
часто  протекает  бессимптомно,  латентно  и  выявить  его  можно  только  с 
помощью  проведения  функциональных  проб.  Одной  из  таких  проб 
является  проба  с  дозированной  физической  нагрузкой,  позволяющая 
определить  оптимальный  режим  физических  нагрузок  для  подростков  с 
различным  уровнем  здоровья.  В  этом  плане  актуальна  разработка 
математической  модели  прогнозирования  функционального  состояния 
миокарда  у  пациентов,  которая  позволит  определить  больных  с  адекватной 
реакцией  сердечно-сосудистой  системы  на  физическую  нагрузку  и  с 
неадекватной,  т.е.  систолической  дисфункцией. 

Этап  1.  Сбор  информации. 

Для  определения  систолической  дисфункции  миокарда  у  подростков 
была  исследована  совокупность  объектов  (137  пациентов  10-18  лет), 
разделенная  на  три  группы:  1-я  группа  -  пациенты  с  адекватной  реакцией 
на  физическую  нагрузку  (52  подростка);  2-я  группа  -  пациенты  с 
неадекватной  реакцией  на  физическую  нагрузку  (23  подростка);  3-я  группа 
-  контрольная  (62  здоровых  сверстника). 

Этап  2.  Определение  оптимального  объема  выборки. 

Так  как  из  предварительного  исследования  известно,  что  среди  всех 
пациентов  больные  с  позитивной  динамикой  составляют  37,9%,  при 
уровне  надежности  у  =  0,95  и  предельной  ошибки  выборки  5-6% 
необходимый  объем  выборки  составляет: 


п 


со 


I2  ■  СО  -  (1  —  со) 


119  ч- 168  человек, 


где  і  -  стандартизированный  коэффициент  Стьюдента. 
Этап  3.  Кодирование  информации. 
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Для  определения  признаков,  которые  существенно  влияют  на 
процесс  дифференциальной  диагностики  систолической  дисфункции 
миокарда,  было  исследовано  83  признака,  включающие  в  себя  паспортные 
данные,  данные  анамнеза  заболевания  и  жизни,  клинические  симптомы, 
показатели  эхокардиографии  (после  нагрузки)  и  т.д..  Все  признаки  были 
закодированы  и  поставлены  в  соответствие  83 -мерному  вектору, 
учитывающему  отсутствие,  присутствие,  направленность  и  величину 
каждого  показателя.  Кодировочная  таблица  представлена  в  табл.  Б.1  [29]. 

Для  ввода  данных  двойным  щелчком  левой  клавиши  мыши  на 
значке  8ТАТІ8ТІСА  10  запустили  программу  8ТАТІ8ТІСА.  В 
открывшемся  редакторе  данных  (рис.  2.1)  ввели  все  подготовленные 
данные.  Для  этого  дважды  нажимали  на  ячейке  с  надписью  «Ѵаг»,  после 
чего  и  появилось  окно  «ѴагіаЫе»,  в  котором  ввели  всю  информацию  о 
переменных. 


Рисунок  2.1  -  Редактор  данных:  просмотр  данных 

Для  кодирования  качественных  переменных,  например,  показателя 
гагпіса,  по  которому  мы  и  будем  группировать  исследуемых  пациентов, 
нажали  на  кнопку  «Техі  ЬаЬеІв»  диалогового  окна  редактирования 
переменных  (рис.  2.2).,  после  чего  открывается  диалоговое  окно  «Техі 
ЬаЬеІв  Еёііог». 
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Рисунок  2.2  -  Диалоговое  окно  «ѴагіаЫе» 

Метки  значений  определялись  следующим  образом  (рис.  2.3): 

1)  вначале  ввели  в  поле  «Тех*  ЬаЬеІ»  название  1  группы 
«положительная  нагрузка»,  а  затем  в  поле  «ІЧитегіс»  число  «1»; 

2)  повторили  эти  действия  для  значений  «2»  —  «отрицательная 
нагрузка»  и  «3»  —  «норма». 

Таким  образом,  были  закодированы  все  качественные  показатели. 

Результат  ввода  всех  значений  в  диалоговом  окне  показан  на  рис.  2.4. 

Этап  4.  Определение  дискриминантных  переменных,  включаемых  в 
модель. 

Задача  определения  количества  дискриминантных  переменных 
решается  путем  процедуры  последовательного  отбора. 

Нам  необходимо  принять  решение  по  отбрасыванию  переменных, 
которые  заведомо  не  могут  быть  хорошими  предикторами  [30]. 

После  ввода  всех  необходимых  данных  для  проведения 
дискриминантного  анализа  выбрали  в  меню  «Зіайвйсв»  пункт  «Микіѵагіаіе 
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Ехріогаіогу  ТесЬпщиез  — >  Вівсгітіпапі;  Апаіувіз»  (рис.  2.5),  после  чего 
открывается  диалоговое  окно  «Оізсгітіпапі;  РипсГіоп  Апаіузів»  (рис.  2.6). 


Рисунок  2.3  -  Диалоговое  окно  «Техі  ЬаЬеЬ  Есіііог» 


0  6ТАТ15ТІСА  -  [Рай:  $ргеа<І5Нее(2  (89ѵ  Ьу  137с)] 


ЕІІе  Е8Л  Ѵіеѵѵ  Іп5егі  Еогтаі  $1аІІ5ІІС5  Оа(а  М'піпд  СгарЬі  ТооІ5  Даіа  Епіегргіае  ѴѴіпсІоѵѵ  Неір 

М  Асісі  Іо  ѴѴогкЬоок  7  АсЗсІ  Іо  Керогі  ^  АсІсЗ  Ео  М5  ѴѴогсІ ' 


АгіаІ 


в  /  і  А-  » 


Й  Й"  ^  &  гі  “=?  Ѵаг5  ’  Са5К  ’ 


1 

пот 

2 

Оіаалог 

3 

0ІЭ0П02І 

4 

5П 

5 

Іѵ 

6 

ека  піѵ 

7 

(і/  піѵ 

8 

ГЭ2ПІСЭ 

ѴАК0000 

10 

ЛОІ 

11 

Ѵ02Г35І 

12 

№ѵѵѴаг 

13 

ека  сИзз 

14 

т  Іеіа 

V 

1 

33 

Миокард;  Миокард; 

Без  серд 

66 

72 

7 

положите 

1 

м 

15 

50 

50 

1.77 

2 

53 

Миокард;  Миокард; 

Без  серд 

57 

66 

9 

положите 

1 

м 

17 

69 

69 

1.45 

3 

57 

аритмия  Дизпгаст 

Без  серд 

48 

57 

9 

положите 

1 

м 

11 

67 

67 

1.30 

4 

89 

Миокард;  Миокарді 

Без  серд 

63 

65 

2 

отрицатег 

2 

м 

17.1 

73 

73 

1.82 

5 

182 

Миокард;  Миокарді 

Без  серд 

61 

61 

0 

отрицатег 

2 

м 

14 

87 

87 

1.96 

6 

207 

Синдром  Синдром 

Без  серд 

54 

73 

19 

положите 

1 

ж 

15 

49 

49 

1.50 

7 

216 

аритмия  Дизпгаст 

Без  серд 

65 

73 

9 

положите 

1 

ж 

15 

76 

76 

1.50 

8 

220 

Миокарді  Миокарді 

Без  серд 

67 

80 

13 

отрицатег 

1 

ж 

15 

52 

52 

1,75 

9 

238 

аритмия  Дизпгаст 

Без  серд 

57 

66 

9 

положите 

1 

м 

17 

111 

111 

1,71 

10 

287 

Миокарді  Миокарді 

Без  серд 

59 

74 

16 

положите 

1 

м 

13 

66 

66 

1.85 

11 

309 

стеноз  ао  Пороки  с 

Без  серд 

57 

60 

3 

положите 

1 

м 

18 

65 

65 

1.69 

12 

315 

Коарктаці  Пороки  с 

Без  серд 

31 

47 

16 

положите 

1 

ж 

17 

79 

79 

1.63 

13 

339 

стеноз  ле  Пороки  с 

Без  серд 

66 

67 

1 

положите 

2 

м 

15 

63 

63 

1.55 

14 

340 

стеноз  ле  Пороки  с 

Без  серд 

68 

77 

10 

положите 

1 

м 

16 

59 

59 

1.61 

15 

348 

стеноз  ао  Пороки  с 

Без  серд 

69 

78 

9 

положите 

1 

м 

17 

62 

62 

1.79 

16 

358 

стеноз  ле  Пороки  с 

Без  серд 

65 

71 

6 

положите 

1 

м 

16 

55 

55 

1.46 

Рог  Неір.  ргеаз  Р1 


5ргеа8а1іееі2  С67.Ѵ12 


58  ДеІОД.  УУадЫДРР  ;  д^ЩЖ  I  НЕС 


Рисунок  2.4  -  Внесенные  данные  в  редактор  данных  ЗМііДіса 
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Рисунок  2.5  -  Выбор  метода  анализа  данных 

В  диалоговом  окне  «ВІ8сгітіпап1:  Рипсііоп  Апа1у8І8»  (  рис.  2.6)  нажали 
на  кнопку  «ѴагіаЫев».  В  открывшемся  окне  «8е1ес1:  опе  §гоиріп§  ѵаг.  апсі 
іпсіерепсіепі  ѵагіаЫе  1І8І;»  выбрали  группирующую  переменную  «гагпіса»  и 
необходимые  независимые  переменные  (рис.  2.7)  и  нажали  ОК. 


Рисунок  2.6  -  Диалоговое  окно  «Візсгітіпапі;  Ршісйоп  Апа1у8І8» 
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Рисунок  2.7  -  Диалоговое  окно  «8е1есі  опе  §гоиріп§  ѵаг.  апё 
іпсіерепсіепг  ѵагіаЫе  1І8І» 


В  открывшемся  окне  (рис.  2.8)  «ѴагіаЫез  сопіаіп  Іехі  ѵаіиез/іехі 
1аЬеІ8»  выбрали  из  предложенного  перечня  пункт  «Сопііпие  \ѵітН  сиггепі 
зеіесііоп». 


ѴагіаЫез  сопіаіп  Іехі  ѵаіиез/іех!  ІаЬеІз 


ЗеІесІеЬ  ѴагіаЫез  ѵѵііН  Іех!  ѵаІиезДех!  ІаЬеІз 
Тех!  ѵагіаЫез:  попе 

Тех!  ІаЬеІз:  82  81  18  27  26  47  60  64-65  68  80  13 


Уои  Наѵе  зеіесіесі  ѴагіаЫез  ІогіНе  сиггепі  Апаіузіз  іНа!  аге  Ы  !уре  Тех!'  ог  сопіаіп 
Тех!  І_аЬеІ5  Іог  а  гоіе  іп  ІЬе  Апаіузіз  Іог  ѵѵЬісН  сопііпиоиз  питегіс  ѵаіиез  аге  ехресіесі. 

РІеазе  зеіес!  Ігот  ІЬе  ІоІІоѵѵіпд  орііопз.  №іе  іЬз!  зеіесііпд  ІЬе  "Сопііпие  ѵчтііі  сиггепі 
зеіесііоп"  (г!  аѵаіІаЫе}  тау  ІеасІ  !о  ипехресіесі  гезиКз  іі  питегіс  ѵаіиез  ѵѵеге  аззідпесі 
!о  Тех!  ІаЬеІз.” 


Яетоѵе  ЯаддеЬ  ѵагіаЫез  Ігот  зеіесііоп 


Сопііпие  ѵѵііЬ  сиггепі  зеіесііоп 


5сап  ЗргеасізЬееІ  Іог  Іехі  ІаЬеІз 


(РІеазе  зее  ІЬе  Неір  Іоріс 
аззосіаІеЬ  ѵяіЬ  іЬіз  Ьіаіод  Іог 
тоге  сіеіаііз  аЬои!  ѴагіаЫе 
Іурез  (е.д..  Тех!.  БоиЫе)  апЬ 
Тех!  ІаЬеІз). 


Рисунок  2.8  -  Диалоговое  окно  «ѴагіаЫе8  сопіаіп  Іехі  ѵаіиез/іехі 

1аЬеІ8» 
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В  диалоговом  окне  «Вівсгітіпапі  Рітсйоп  Апаіувів»  (рис.  2.6)  нажали 
на  кнопку  «Сойев  іог  §гоиріп§  ѵагіаЫе».  В  появившемся  окне  «8е1есі  сойев 
{'о г  §гоиріп§  ѵагіаЫе»  выбрали  все  3  сравниваемые  группы  с  помощью 
нажатия  кнопки  АН  (рис.  2.9)  и  нажали  ОК. 


Рисунок  2.9  -  Диалоговое  окно  «8е1ес1  сойев  іог  §гоиріп§  ѵагіаЫе» 

В  диалоговом  окне  «Вівсгітіпапі  Рітсйоп  Апаіувів»  (рис.  2.6) 
установили  «галочку»  на  опции  Асіѵапсе  орііот  (зіертзе  апаіузіз)  и 
нажали  на  кнопку  ОК,  открылось  окно  диалога  «Мойеі  Вейпійоп» 
(рис.  2.10). 


Рисунок  2.10  -  Диалоговое  окно  «Мойеі  Вейпійоп» 

Перейдя  на  вкладку  «Аёѵапсесі»  (рис.  2.11),  в  поле  «Мейюй» 
указали  метод  дискриминантного  анализа:-  ГогхѵаМ  8іер\ѵІ8е  (пошаговый 
вперед). 
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Рисунок  2.11  -  Диалоговое  окно  «Моёеі  БейпШоп».  Вкладка 

«Асіѵапсесі» 

При  данном  выборе  программа  на  последовательных  шагах  включает 
переменные  в  модель. 

Опция  Тоіеюпсе  (толерантность)  задает  нижнюю  границу 
толерантности.  Толерантность  является  мерой  избыточности  переменные. 
Чем  меньше  значение  толерантности,  тем  избыточнее  переменная  в 
модели,  так  как  переменная  несет  малую  дополнительную  информацию 
[31].  Переменные  с  толерантностью  меньше  заданного  значения  в  модель 
не  включаются. 

Тест  толерантности  может  обеспечить  точность  вычислений. 
Толерантность  еще  не  отобранной  переменной  равна  единице  минус 
квадрат  множественной  корреляции  К 2  между  этой  переменной  и  всеми 
уже  отобранными  переменными,  когда  корреляции  определяются  по 
внутригрупповой  корреляционной  матрице. 

Если  проверяемая  переменная  является  линейной  комбинацией  (или 
приблизительно  равна  линейной  комбинации)  одной  или  нескольких 
отобранных  переменных,  то  ее  толерантность  равна  нулю  (или  близка  к 
нулю).  Переменная  с  малой  толерантностью  (скажем,  меньше  0,001)  может 
привести  к  ошибке  при  вычислении  матрицы,  обратной  Е,  ввиду  быстрого 
накопления  ошибок  округления.  Помимо  вычислительных  проблем, 
нежелательно  использовать  переменную,  которая  является  линейной 
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комбинацией  отобранных  переменных,  потому  что  она  не  дает  никакой 
новой  информации. 

Статистика  ^-включения  представляет  собой  частную  Е-статистику, 
оценивающую  улучшение  различения  от  использования  рассматриваемой 
переменной  по  сравнению  с  различением,  достигнутым  с  помощью  других 
уже  отобранных  переменных.  Если  величина  статистики 
Е-включения  мала,  то  такая  переменная  не  будет  отобрана,  потому  что  она 
не  дает  достаточно  большого  вклада  в  различение.  Переменная  должна 
пройти  проверку  толерантности  и  Е-включения,  прежде  чем  она  будет 
рассмотрена  в  соответствии  с  критерием  отбора. 

Статистика  Е-удаления  также  является  частной  Е-статистикой. 
Однако  она  оценивает  значимость  ухудшения  различения  после  удаления 
переменной  из  списка  уже  отобранных  переменных.  Эта  процедура 
проводится  в  начале  каждого  шага,  чтобы  проверить,  имеется  ли  какая- 
нибудь  переменная,  уже  не  вносящая  достаточно  большого  вклада  в 
различение,  поскольку  отобранные  позже  переменные  дублируют  ее  вклад. 
На  заключительном  шаге  статистика  Е-удаления  может  быть  использована 
для  ранжирования  дискриминантных  возможностей  отобранных 
переменных.  Переменная  с  наибольшим  значением  статистики  Е-удаления 
даст  наибольший  вклад  в  различение,  достигнутое  благодаря  другим 
переменным.  Переменная,  имеющая  вторую  по  величине  статистику  Е- 
удаления,  является  второй  по  значимости  и  т.  д. 

После  выбора  в  диалоговом  окне  «Мобеі  РеРіпШоп»  метода  РопѵаЫ 
зіертзе  и  нажатия  по  ОК,  открылось  окно  результатов  (рис.  2.12). 

Информационная  часть  окна  «РНслгтпаЩ  РипсНоп  Апаіузіз  КезиЙв» 
сообщает,  что: 

-  МитЬег  о/  ѵагіаЫек  іп  тосіеі  (число  переменных  в  модели)  равно  12; 

-  ШІЬ  ЬатЪсІа  (значение  лямбда  Уилкса)  равно  0,1070642; 

-  арргох.Р  Приближенное  значение  Е-статистики  с  числом  степеней 
свободы  2  и  123  равно  1,725747,  а  с  числом  степеней  свободы  24  и  246  - 
21,07576. 

Статистика  лямбда  Уилкса  ( X )  вычисляется  как  отношение 
детерминанта  (определителя)  матрицы  внутригрупповых  ковариаций  к 
детерминанту  общей  ковариационной  матрицы.  Значения  X  принадлежат 
интервалу  [0,1].  Значения  X,  лежащие  около  0,  свидетельствуют  о  хорошей 
дискриминации  [32-33]. 
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Рисунок  2.12  -  Диалоговое  окно  «Вівсгітіпапі  Рипсііоп  Апа1у8І8 

Ке8и1і8» 

Далее  в  диалоговом  окне  «Оізсгітіпапі  Рипсііоп  Апа1у8І8  Ке8и1і8» 
(рис.  2.12)  перешли  на  вкладку  «Асіѵапсесі»  и  нажали  на  кнопку  «Зшптагу: 
ѴагіаЫез  іп  Ійе  тоёеі»  (итоги:  переменные,  включенные  в  модель) 
(рис.  2.13).  Появится  итоговая  таблица  анализа  данных  (рис.  2.14). 

После  нажатия  на  кнопку  «8иттагу:  ѴагіаЫез  іп  іЬе  тоёеі» 
(итоги:  переменные,  включенные  в  модель)  появляется  итоговая  таблица 
анализа  данных  (рис.  2.13). 

В  первом  столбце  таблицы  приведены  значения  УѴіІкз  ЬатЪсІа 
(лямбда  Уилкса),  являющиеся  результатом  исключения  соответствующей 
переменной  из  модели.  Чем  больше  значение  X,  тем  более  желательно 
присутствие  этой  переменной  в  процедуре  дискриминации. 

Значение  Рагііаі  ЬатЪсІа  (частная  лямбда)  есть  отношение  лямбда 
Уилкса  после  добавления  соответствующей  переменной  к  лямбде  Уилкса 
до  добавления  этой  переменной.  Частная  лямбда  характеризует  единичный 
вклад  соответствующей  переменной  в  разделительную  силу  модели.  Чем 
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меньше  статистика  лямбда  Уилкса,  тем  больше  вклад  в  общую 
дискриминацию . 


Рисунок  2.13  -  Диалоговое  окно  «Вівсгітіпапі  Рипсііоп  Апаіувів 
Яешіід».  Вкладка  «Асіѵапсесі» 


N=137 

Оізсгітіпапі  Рипсііоп  Апаіузіз  Зиттагу  (5ргеасізЬееі2) 

Зіер  12,  N  сЯѵагз  іп  тосіеі:  12;  Сгоиріпд:  гагпіса  (3  дгрз) 

ѴѴІІкз’  І атЬс)а:  .10706  арргох  Р  (24.246)=21,076  р<0.0000 
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0.000000 
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0  310839 

Р* 
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0.782612 

17.08298 

0.000000 

0.702534 

0.297466 

гааз  апд 
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0.099818 
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0.131742 
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0.159918 

ікзо 
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0.838368 
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0.646519 
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0.013782 

0.746399 

0.253601 

кт  погасі 

0.113107 

0.946578 
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0.034168 

0.920736 

0.079264 
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0.036613 
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Рисунок  2.14  -  Итоги:  переменные,  включенные  в  модель 


Г-гетоѵе  (Р-удаление)  -  это  значения  Г -критерия,  связанные  с 
соответствующей  частной  лямбдой  Уилкса.  Значения  р-іеѵеі  -  это  уровни 
значимости  критериев  Г-гетоѵе. 


51 


Наибольший  вклад  в  различение  дисфункции  миокарда  дает 
поверхность  тела  пациента  (рис.2.14). 

Этап  5.  Построение  диагностирующей  модели  для  дисфункции 
миокарда. 

В  диалоговом  окне  «Оізсгітіпапі  РипсНоп  Апаіувів  Кевиіів»  (рис. 
2.13)  нажали  кнопку  «Регіогт  сапопісаі  апаіувів»  (выполнение 
канонического  анализа),  программа  вычислила  независимые 
(ортогональные)  дискриминантные  функции. 

В  открывшемся  диалоговом  окне  «Сапопісаі  Апаіувів»  (рис.  2.15) 
нажали  кнопку  «СоеШсіепів  Іог  сапопісаі  ѵагіаЫез»  (коэффициенты 
канонических  переменных). 

В  результате  появились  две  таблицы  с  коэффициентами 
дискриминантных  (канонических)  функций.  В  первой  таблице 
представлены  исходные  (нестандартизированные)  коэффициенты 
дискриминантных  функций  (рис.  2.16,  а).  Эти  коэффициенты  могут  быть 
использованы  для  вычисления  значений  канонических  переменных  для 
каждого  наблюдения  каждой  дискриминантной  функции. 


Рисунок  2.15  -  Диалоговое  окно  «Сапопісаі  Апаіузів» 

Во  второй  таблице  (рис.  2.16,  б)  приведены  стандартизованные 
коэффициенты  дискриминантных  функций.  Эти  коэффициенты, 
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основанные  на  стандартизованных  переменных,  принадлежат  к  одной  и 
той  же  шкале  измерений  (абсолютной),  поэтому  их  можно  сравнивать, 
чтобы  определить  величины  и  направления  вкладов  переменных  в  каждую 
каноническую  функцию. 


ѴагіаЫе 

Раш  СоеЛісіепІз  (ЗргеасІзНееЬ?) 
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Рисунок  2.16  -  Коэффициенты  дискриминантных  функций: 
а)  нестандартизированные  и  б)  стандартизированные 

Стандартизованные  коэффициенты  показывают  вклад  переменных  в 
значение  дискриминантной  функции.  Это  является  одним  из  подходов  к 
определению  значимости  переменной.  Однако  этот  подход  имеет 
серьезное  ограничение.  Если  две  переменных  сильно  коррелированны,  то 
их  вклад  в  дискриминантное  значение  должен  разделяться,  даже  при 
значительном  совместном  вкладе.  Соответственно  их  стандартизованные 
коэффициенты  могут  быть  меньше  по  сравнению  с  теми  случаями,  когда 
используется  одна  из  этих  переменных.  Или,  другими  словами,  вклад 
одного  коэффициента  частично  погашается  отрицательным  вкладом 
другого.  Это  происходит  из-за  того,  что  в  стандартизованных 
коэффициентах  одновременно  принимается  во  внимание  влияние  всех 
переменных. 

Таким  образом,  получены  следующие  дискриминантные  функции: 
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т  1  =  -0,061*  хі  -  0,118*  х2+  0,029*  хЗ  -  0,057*  х4  +  0,091*  х5  -  0,004*х6  + 
+2,367*  х7  -  0,019*  х8  -  0,022*  х9+  0,034*  хЮ-1,644*  хі  1  +  0,005*  х12  - 
-1,214; 


Ш2  =  -  0,022*  хі  +  0,152*  х2  -  0,013*  хЗ  +  0,052*  х4  +  0,138*  х5  +  0,003* 
*  хб  +  3,794  *  х7  +  0,718*  х8  +  0,0008*  х9  -  0,004*  хЮ  +  0,281  *  хіі  + 

+  0,007*  х12-  16,378, 

где  хі  -  фактор  некроза  опухолей  (ФНО-а,  пг  /  мл); 
х2  -  СБ-95  (%  клеток,  несущих  маркер  апоптоза); 
хЗ  -  циклический  аденозинмонофосфат  (цАМФ,  нмоль  /  мл); 
х4  -  ИКСО  (мл/м2) 

х5  -  показатель  сокращения  левого  желудочка  (%); 

хб  -  скорость  транс-аортального  кровотока  (см  /  с); 

х7  -  поверхность  тела  (м2); 

х8  -  ренин  (нг  /  мл  /  ч); 

х9  -  ангиотензин-П  (пмоль  /  л); 

хЮ  -  простациклин  (пг  /  мл); 

хіі  -  правый  желудочек  (см); 

х12  -  норадреналин  (нмоль  /  сут). 

В  диалоговом  окне  «Сапопісаі  Апа1у§І8»  выбрали  вкладку 
«Асіѵапсесі»  (рис.  2.17),  нажали  кнопку  «Расіог  вігасйіге»  (факторная 
структура).  В  появившейся  таблице  (рис.  2.18)  приведены  объединенные 
внутригрупповые  корреляции  переменных  с  соответствующими 
дискриминантными  функциями. 

Эти  корреляции  называют  еще  структурными  коэффициентами. 
Обычно  структурные  коэффициенты  используют  для  содержательной 
интерпретации  функций,  в  отличие  от  коэффициентов  дискриминантной 
функции,  которые  обозначают  вклад  каждой  переменной  в  функции. 
Структурный  коэффициент  можно  вычислить  и  для  переменной  в  пределах 
отдельно  взятой  группы.  Тогда  получаем  внутригрупповой  структурный 
коэффициент  2ц  для  і- й  переменной  /- й  функции,  который  вычисляется  по 
формуле: 
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где  г,у  -  внутригрупповые  структурные  коэффициенты  корреляции 
между  переменными  і  и  у;  с#  -  стандартизованные  коэффициенты 
канонической  функции  для  переменной  у  и  функции  /  [20]. 


СапопісаІ  Апаіухіз:  $ргеасІ5Нееі2 


(Зиіск  Асіѵапсесі  СапопісаІ  $соге$  | 

тшЗ  Зиттащ 

ттш  Зигптагѵ:  СЫ  $диаге  1е$($  о(  $иссе$$іѵе  гооІ$ 

СапсеІ 

Соейісіепіз  (ог  сапопісаІ  ѵагіаЫе$ 

[Я  ОрЬоп$  ▼ 

Ш  Расіог  зігисіиге 

Д|  Вѵ  Згоир 

[Щ  Меап$  о(  сапопісаІ  ѵагіаЫе$ 

Рисунок  2.17  -  Диалоговое  окно  «СапопісаІ  Апаіувів». 

Вкладка  «Асіѵапсесі» 

С  другой  стороны,  структурные  коэффициенты  являются  просто 
двуместными  корреляциями,  поэтому  на  них  не  влияют  взаимные 
зависимости  прочих  переменных. 

В  работе  [34]  описывается,  как  структурные  коэффициенты  могут 
использоваться  для  графического  представления  различия  между 
групповыми  центроидами  в  случае  двух  канонических  дискриминантных 
функций.  На  графике  с  осями,  которые  относятся  к  этим  двум  функциям, 
представлены  групповые  центроиды  и  главный  центроид,  изображены 
векторы,  исходящие  из  главного  центроида  и  направленные  в  каждую 
дискриминантную  переменную.  Направляющие:  углы  этих  векторов 
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вычисляются,  исходя  из  структурных  коэффициентов.  Длина  вектора 
определяется  межгрупповыми  и  внутригрупповыми  вариациями 
соответствующей  переменной. 


Расіог  Зігисіиге  Маігіх  (ЗргеасізІіееі2) 
Соггеіаііопз  ѴагіаЫез  -  Сапопісаі  Коо!з 
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ѴагіаЫе 

Кооі  1 

Рооі  2 

5(195 

-0.3511531 

0.176843 

РГ02  СІ5 

0.318809 

0.308907 

роѵ  Іеіа 

0.057711 

0.532739 

Р* 

-0.169708 

0.348437 

гааз  апд 

-0.300870 

0.088785 

іт  Іпоа 

-0.284637 

-0.061710 

ргозіаіз 

0.243468 

-0.091718 

ікзо 

-0.242858 

0.080039 

Ізаті 

0.003884 

-0.088459 

кт  погас! 

0.012844 

0.253040 

гааз  геп 

0.006142 

0.270468 

сІ1_а  к_ѵ 

-0.085334 

0.168433 

Рисунок  2.18  -  Объединенные  внутригрупповые  корреляции  переменных  с 
соответствующими  дискриминантными  функциями 


Полученная  диаграмма  дает  наглядное  представление  о  различиях 
групп  с  помощью  дискриминантных  переменных,  а  также  о 
потенциальных  возможностях  этих  переменных. 

В  окне  «Сапопісаі  Апаіузіз»  нажали  на  кнопку  «8шшпагу:  СЫ  вциаге 
Іезів  оР  виссеввіѵе  гооів»  (итоги:  хи-квадрат  критерий  последовательности 
корней)  (рис.  2.15).  Появится  таблица  результатов  с  пошаговым  критерием 
для  канонических  корней  —  дискриминантных  функций  (рис.  2.19). 


Кооіз 

Кетоѵесі 

СЫ-Здиаге  Тезіз  ѵѵіі 

т  Зиссеззіѵе  Кооіз  КетоѵесІ  (Зрге 

Еідеп- 

ѵаіие 

СапопісІ 

К 

ѴѴіІкз’ 

ІатЬсІа 

СИі-Здг. 

(И 

р-ѵаіие 

0 

3  9765121  0.893899 

0,107064 

287,1109 

24 

0.000000 

1 

0  876854я  0.683516 

0.532806 

80.9032 

11 

0.000000 

Рисунок  2.19  -  Таблица  результатов  с  пошаговым  критерием  для 

канонических  корней 

В  диалоговом  окне  «Сапопісаі  Апа1у8І8»  перешли  на  вкладку 
«Сапопісаі  8СОге8»  (канонические  значения)  и  нажали  кнопку  «Сапопісаі 
8СОге8  Рог  еасЪ  сазе»  (канонические  значения  для  каждого  наблюдения) 
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(рис.  2.20).  Появится  таблица  (рис.  2.21)  со  значениями  дискриминантных 
функций  для  каждого  наблюдения. 


Рисунок  2.20  -  Диалоговое  окно  «Сапопісаі  Апаіувів». 
Вкладка  «Сапопісаі  всогев» 


Упзіапсіагсіігесі  Сапопісаі  Зсогез  (ЗргеасІзИееі2) 


Сазе 

Сгоир 

Кооі  1 

Кооі  2 

1 

попожитепьная  физнагрѵзкаі  -1.77625 

0.30163 

2 

попожитепьная  физнагрузка 

-4.57447 

-0.79041 

3 

положительная  физнагрузка 

-4. 31024 

-1.24178 

4 

отрицательная  физнагрузка 

-0.90177 

2.46267 

5 

отрицательная  физнагрузка 

-1.67667 

1.17803 

6 

положительная  физнагрузка 

-0.94063 

-1.28100 

7 

положительная  физнагрузка 

-3.71312 

-0.28700 

8 

отрицательная  физнагрузка 

0.16817 

0.57280 

9 

положительная  физнагрузка 

-3.43594 

1.19067 

10 

положительная  физнагрузка 

-0.59223 

0.52584 

11 

положительная  физнагрузка 

-2.83045 

-0.73086 

12 

положительная  физнагрузка 

-2.52023 

-0.68473 

13 

положительная  физнагрузка 

-0.52004 

0.04671 

14 

положительная  физнагрузка 

-2.72984 

0,12341 

15 

положительная  физнагрузка 

-2.97920 

-1.61743 

16 

положительная  физнагрузка 

-1.63367 

-0.43257 

17 

положительная  физнагрузка 

-4.44589 

-0,78120 

136 

норма 

2.53059 

-0.64247 

137 

норма 

1 ,67824 

0,27540 

Рисунок  2.21  -  Таблица  со  значениями  дискриминантных  функций 

для  каждого  наблюдения 

В  окне  «Сапопісаі  Апаіувів»  нажали  кнопку  «Зсайегріоі  ей  сапопісаі 
всогев»  (диаграмма  рассеяния  для  канонических  значений).  Появилась 


57 


диаграмма  рассеяния  (рис.  2.22)  канонических  значений  для  пар  значений 
дискриминантных  функций.  На  диаграмме  видно,  что  наблюдения, 
принадлежащие  одинаковым  группам,  локализованы  в  определенных 
областях  плоскости. 


ГСоо1 1  Ѵ5.  ГСооІ  2 


Кооі  1 


о  положительная  физнагрузка 
□  отрицательная  физнагрузка 
о  норма 


Рисунок  2.22  -  Диаграмма  рассеяния  исследуемых  групп  пациентов 

Вернувшись  в  окно  результатов  «Ойсгітіпапі  Рипсйоп  Апаіузіз 
Кезиііз»,  нажали  кнопку  «Сапсеі»,  и  активизировали  вкладку 
«Сіаввійсайоп»  (классификация).  Открылось  окно  результатов 
классификации  (рис.  2.23). 

В  рамке  А  ргіогі  сіаззі/ісаііоп  ргоЬаЫІіііез  (априорные  вероятности 
классификации)  приведены  различные  опции  задания  априорных 
вероятностей  того,  что  наблюдение  при  классификации  попадет  в  одну  из 
групп:  Ргорогііопаі  іо  §гоир  зііез  (пропорциональные  размерам  групп); 
8ате  /ог  аіі  §гоирз  (одинаковые  для  всех  групп);  ІІзег  сіе/іпесі  (заданные 
пользователем).  Априорные  вероятности  могут  существенно  влиять  на 
точность  классификации.  Если  есть  предварительные  сведения  (оценки)  о 
возможном  количественном  соотношении  наблюдений  в  группах,  то 
желательно  выбрать  опцию  ІІзег  сіе/іпесі,  если  таких  сведений  нет  и  число 
наблюдений  в  группах  примерно  одинаково,  то  надо  выбрать  8 сипе  /о г  аіі 
§гоирз,  в  противном  случае  -  Ргорогііопаі  іо  §гоир  зеез. 
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Рисунок  2.23  -  Диалоговое  окно  «Вівсгітіпапі;  Рипсйоп  Апаіувів 
КезиЙз».  Вкладка  «Сіаввійсайоп» 

В  рамке  8соге  іо  заѵе  /ог  еаск  сазе  (сохранить  для  каждого 
наблюдения)  приведены  опции,  при  выборе  которых  можно  сохранить  тот 
или  иной  результат  классификации:  8аѵе  сіаззі/ісаііоп  /ог  сазе  (сохранить 
результаты  классификации);  8аѵе  сіізіапсе /ог  сазе  (сохранить  расстояния); 
8аѵе  розіегіог  ргоЪаЫШу /ог  сазе  (сохранить  апостериорные  вероятности). 

В  строке  Мах.  питЪег  о/  сазез  іп  а  зіп§1е  гезиііз  зргеасізкееі 
(максимальное  число  наблюдений  в  таблице  результатов).  Если 
наблюдений  больше  указанного  числа,  то  результаты  будут  выделены 
несколькими  таблицами. 

В  диалоговом  окне  «Оізсгітіпапі:  Рипсйоп  Апаіувів  КезиЙз»  на 
вкладке  «Сіаззійсайоп»  нажали  кнопку  «СІаввіРісайоп  йдпсйопв»  (функции 
классификации).  Функции  классификации  -  это  линейные  функции, 
которые  вычисляются  для  каждой  группы  и  могут  быть  использованы  для 
классификации  наблюдений  (пациентов).  Исследуемого  пациента 
приписывают  к  той  группе,  для  которой  классификационная  функция 
имеет  наибольшее  значение. 
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В  таблице,  изображенной  на  рис.  2.24,  приведены  коэффициенты  и 
свободные  члены  при  переменных  линейных  функций.  Например, 
классификационные  уравнения  для  всех  3-х  групп  имеют  вид: 

Группаі  =  -  хі*  0,292+  х2  *  3,588+  хЗ  *  0,089+  х4  *  1,202+  х5* 
*2,546+  хб*  0,060+  х7*  61,200+  х8*  6,133+  х9*  0,067+  х10*  0,253  +  хіі* 
*5,770  +  х12*  0,049  -  159,393. 

Группа2  =  -  хі*  0,452+  х2  *  3,789+  хЗ  *  0,104+  х4  *  1,244+  х5* 
*3,052+  хб*  0,061+  х7*  74,958+  х8*  7,971+  х9*  0,032+  х10*  0,298  +  хіі* 
*3,807  +  х12*  0,075  -  204,465. 


ГруппаЗ  =  -  хі*  0,563+  х2  *  3,118+  хЗ  *  0,213+  х4  *  0,969  +  х5* 
*2,975+  хб*  0,043+  х7*  72,414+  х8*  6,245-  х9*  0,029+  хЮ*  0,401  -  хіі* 
*  1,228  +  х12*  0,071  -  168,451 


где  группаі  -  пациенты  с  адекватной  реакцией  на  физическую 
нагрузку;  группа2  -  пациенты  с  неадекватной  реакцией  на  физическую 
нагрузку;  группаЗ  -  контрольная  группа  . 


ѴагіаЫе 

Сіаззійсаііоп  Рипсііопз;  дгоиріпд:  гагпіса  (5| 

положительна 
я  физнагрузка 
р=, 37956 

отрицательная 
физнагрузка 
р=  16788 

норма 

р=:45255 

зсІ95 

3.588 

3.789 

3.118 

РГ02  СІ5 

2.546* 

3.052 

2.975 

роѵ  Іеіа 

61.200 

74.958 

72.414 

Рг 

5.770 

3.807 

-1.228 

гааз  апд 

0.067 

0.032 

-0.029 

іт  Іпоа 

-0.292 

-0.452 

-0.563 

ргозіаіз 

0.253 

0.298 

0.401 

ікзо 

1.202 

1.244 

0.969 

Ізат^ 

0.089 

0.104 

0.213 

кт  погасі 

0.049 

0.075 

0.071 

гааз  геп 

6.133 

7.971 

6.245 

<К  ак  ѵ 

0.060 

0.061 

0.043 

Сопзіапі 

-159.393 

-204.465 

-168.451 

Рисунок  2.24  -  Классификационные  функции 
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В  случае  двух  групп  разделяющая  точка  С  (константа 
дискриминации)  между  двумя  группами  равна  полусумме  величин  двух 
центроидов  этих  групп. 

В  случае  трех  групп  разделяющие  точки  С1,  С2  ,  С  3. 

Правило  классификации  имеет  вид:  если  1"  >  С  1  -  пациент 

принадлежит  первой  группе;  С1  <  1"  <  С2  -  второй  группе;  1"  <  С  3  -  3-й 
группе. 

В  диалоговом  окне  «Оізсгітіпапі  РііпсПоп  Аиаіувів  Кезийз» 
(рис.  2.23)  нажали  кнопку  «Сіаввійсайоп  о!"  сазез»  (классификация 
наблюдений).  Программа  вывела  таблицу  классификации  для  каждого 
наблюдения  (рис.  2.25).  Классификации  упорядочены  по  первому,  второму 
и  третьему  выбору.  Столбец  1  содержит  первый  классификационный 
выбор,  т.е.  группу,  для  которой  соответствующее  наблюдение  имеет 
наивысшую  апостериорную  вероятность  и  наибольшее  значение 
классификационной  функции.  Наблюдения,  которые  не  удалось  правильно 
классифицировать,  отмечены  *. 


СІаззіЯсаііоп  оІСазез  (ЗргеасІ5Ьееі2) 

Іпсоггесі  сіаззііісаііопз  аге  тагкесі  ѵѵіІРі  * 

ОЬзегѵесІ 

і 

2 

3 

Сазе 

СІаззіІ. 

р=. 37956 

р=. 16788 

р=.45255 

1 

положительная  физнагрѵзка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

2 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

3 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

4 

отрицательная  физнагрузка 

отрицательная  физнагрузка 

положительная  физнагрузка 

норма 

*  5 

отрицательная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

6 

положительная  физнагрузка 

положительная  физнагрузка 

норма 

отрицательная  физнагрузка 

7 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

*  8 

отрицательная  физнагрузка 

норма 

отрицательная  физнагрузка 

положительная  физнагрузка 

9 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

10 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

11 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

12 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

13 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

норма 

14 

положительная  физнагрузка 

положительная  физнагрузка 

отрицательная  физнагрузка 

_ норма 

132 

норма 

норма 

положительная  физнагрузка 

отрицательная  физнагрузка 

133 

норма 

норма 

отрицательная  физнагрузка 

положительная  физнагрузка 

134 

норма 

норма 

отрицательная  физнагрузка 

положительная  физнагрузка 

135 

норма 

норма 

отрицательная  физнагрузка 

положительная  физнагрузка 

136 

норма 

норма 

отрицательная  физнагрузка 

положительная  физнагрузка 

137 

норма 

норма 

отрицательная  физнагрузка 

положительная  физнагрузка 

Рисунок  2.25  -  Таблица  классификации  для  каждого  наблюдения 


На  вкладке  «Сіаввійсайои»  в  диалоговом  окне  «Вівсгітіпапі  РипсНоп 
Аиаіувів  Кевиіів»  (рис.  2.23)  нажали  кнопку  «Зсщагесі  МаНаІапоЬіз  ёівіаисез» 
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(квадрат  расстояния  Махаланобиса).  Была  выведена  таблица  квадратов 
расстояний  Махалонобиса  каждого  наблюдения  от  центроида  группы  (рис. 
2.26). 


Зциагесі  МаЬаІапоЬіз  Оізіапсез  Ігот  СЗгоир  СепІгоіВз  (Зргеа<1зИееі2) 
Іпсоггесі  сІаззіЛсаІіопз  аге  тагкесі  ѵѵііИ  * 

Сазе 

ОЬзегѵесі 

СІаззіЬ 

положительная 

физнагрузка 

р=,37956 

отрицательная 
физнагрузка 
р=, 16788 

норма 
р=, 45255 

1 

положительная  физнагрѵзка 

8.3923 

11.91088 

22  67144 

2 

положительная  физнагрузка 

17.3010 

35.68066 

56.24844 

3 

положительная  физнагрузка 

14.5754 

34.43840 

51.49346 

4 

отрицательная  физнагрузка 

40.7414 

30.14211 

46.32682 

*  5 

отрицательная  физнагрузка 

12.8190 

11.44853 

25.76708 

6 

положительная  физнагрузка 

15  1661 

24.18356 

23  11127 

7 

положительная  физнагрузка 

11.2858 

24.22021 

42.54937 

*  8 

отрицательная  физнагрузка 

10.5497 

6.28034 

7.95111 

9 

положительная  физнагрузка 

12.5423 

16  87510 

40.62140 

10 

положительная  физнагрузка 

9.5842 

8.05292 

13.55404 

11 

положительная  физнагрузка 

7.6272 

19.97784 

31.53586 

12 

положительная  физнагрузка 

40.4120 

51.50520 

61.62639 

13 

положительная  физнагрузка 

6.7096 

7.43588 

10,31746 

14 

положительная  физнагрузка 

6.8642 

14.43755 

29.44491 

132 

норма 

29.5131 

32.87676 

8.64295 

133 

норма 

34.2684 

25.79112 

7,22740 

134 

норма 

16.1526 

14.02534 

4.20093 

135 

норма 

20.5923 

16.28628 

1.60706 

136 

норма 

29.0732 

23.38361 

6.80447 

137 

норма 

19.3685 

11.69558 

3.93736 

Рисунок  2.26  -  Таблица  квадратов  расстояний  Махалонобиса  каждого 
наблюдения  от  центроида  группы 


Эти  расстояния  аналогичны  квадратам  евклидовых  расстояний,  но 
учитывают  корреляции  между  переменными  в  модели.  Пациента 
приписывают  группе,  к  которой  он  ближе  всего.  Пациенты,  которых  не 
удалось  правильно  классифицировать,  также  отмечены  *. 

Этап  6.  Определение  т-статистики  ошибок  и  доли  правильной 
диагностики. 

На  вкладке  «Сіаввійсайоп»  в  диалоговом  окне  «Ріасгітіпапі  Рітсйоп 
Апаіувів  КезиЙв»  (рис.  2.23)  нажали  кнопку  «Сіаввійсайоп  тайіх»  (матрица 
классификации).  Матрица  (рис.  2.27)  содержит  информацию  о  количестве 
и  проценте  корректно  классифицированных  наблюдений  (пациентов)  в 
каждой  группе.  Строки  матрицы  —  исходные  классы,  столбцы  — 
предсказанные  классы. 
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Огоир 

Сіаззійсаііоп  Маігіх  (5ргеасІзІіееі2) 

Коѵѵз:  ОЬзегѵесІ  сіаззііісаііопз 

Соіитпз:  Ргесіісіесі  сіаззійсаііопз 

Регсепі 

Соггесі 

положительная 
физнагрузка 
р=. 37956 

отрицательная 
физнагрузка 
р=  16788 

норма 

р=.45255 

положительная  физнагрузка 

98  07692 

51 

1 

0 

отрицательная  физнагрузка 

78.260871 

2 

18 

3 

норма 

98.38710 

1 

0 

61 

Тоіаі 

94.89051 

54 

19 

64 

Рисунок  2.27  -  Результаты  классификации 

К  группе  «Положительная  физнагрузка»  относится  52  наблюдения. 
51  наблюдение  (98,07692  %)  было  определено  корректно,  1  ошибочно 
отнесено  к  группе  «Отрицательная  физнагрузка».  В  группе 
«Отрицательная  физнагрузка»  насчитывается  23  наблюдения.  18  из  них 
были  определены  корректно,  2  наблюдения  ошибочно  отнесено  к  группе 
«Положительная  физнагрузка»  и  3  к  группе  «Норма».  В  группе  «Норма» 
насчитывается  62  наблюдения.  Корректно  определены  61,  ошибочно  -  1  к 
группе  «Положительная  физнагрузка». 

Анализ  полученных  результатов  (рис. 2.27)  показал,  что  прогноз  для 
исследуемых  групп  сделан  верно  (корректно  классифицированы  в  общей 
сложности  94,89  %  всех  случав  заболеваний  сердца).  Данное  значение 
получено  отношением  корректно  определенных  наблюдений  к  общему  их 
количеству  (130/137=94,  89).  Как  видно  из  данного  расчета,  не 

учитывались  априорные  вероятности  каждой  группы  пациентов. 

В  пакете  8ТАТІ8ТІСА  есть  возможность  задавать  различные  опции 
априорных  вероятностей.  В  диалоговом  окне  «БщсгишпаЩ  Рипсііоп 
Апа1у§І8  Ке8иЙ8»  (рис.  2.22)  на  вкладке  «С1а88Ійсайоп»  нажали  кнопку 
«Ро8Іегіог  ргоЬаЫ1ійе8»  (апостериорные  вероятности).  В  открывшейся 
таблице  (рис.  2.28)  каждому  пациенту  поставлена  в  соответствие 
вероятность  принадлежности  к  группам.  Эта  вероятность  определяется 
посредством  расстояний  Махалонобиса  и  априорных  вероятностей.  Чем 
дальше  наблюдение  расположено  от  центра  группы,  тем  менее  вероятно, 
что  оно  принадлежит  этой  группе.  Пациента  приписывают  к  той  группе, 
для  которой  имеется  наибольшая  апостериорная  вероятность 
классификации.  Априорные  вероятности  могут  быть  заданы 
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пользователем,  могут  быть  равны  для  всех  групп,  могут  быть 
пропорциональны  размерам  групп. 


Сазе 

Розіегіог  РгоЬаЬІІіІіез  (5ргеасІ5Ьееі2) 

Іпсоггесі  сіаззійсаііопз  аге  тагкесі  ѵѵііН  * 

ОЬзегѵегі 

СІаззіІ. 

положительная 

физнагрузка 

р=,37956 

отрицательная 
физнагрузка 
р=, 16788 

норма 
р=  45255 

1 

положительная  физнагоѵзка 

0.928423 

0.070699 

0  000878 

2 

положительная  физнагрузка 

0.999955 

0.000045 

0.000000 

3 

положительная  физнагрузка 

0.999978 

0.000022 

0.000000 

4 

отрицательная  физнагрузка 

0.011154 

0.988031 

0.000815 

*  5 

отрицательная  физнагрузка 

0.532065 

0,466956 

0  000979 

6 

положительная  физнагрузка 

0.973411 

0.004741 

0.021848 

7 

положительная  физнагрузка 

0.999313 

0,000687 

0.000000 

*  8 

отрицательная  физнагрузка 

0.109756 

0.410417 

0.479827 

9 

положительная  физнагрузка 

0.951760 

0,048239 

0.000001 

10 

положительная  физнагрузка 

0.472827 

0.449718 

0.077455 

11 

положительная  физнагрузка 

0.999073 

0.000919 

0.000008 

12 

положительная  физнагрузка 

0.998248 

0.001722 

0.000029 

13 

положительная  физнагрузка 

0.664919 

0,204546 

0,130535 

14 

положительная  физнагрузка 

0.990058 

0.009928 

0.000015 

132 

норма 

0.000025 

0.000002 

0.999973 

133 

норма 

0.000001 

0,000035 

0.999964 

134 

норма 

0.002120 

0.002716 

0.995165 

135 

норма 

0  000063 

0.000241 

0.999696 

136 

норма 

0.000012 

0.000093 

0.999895 

137 

норма 

0.000371 

0.007606 

0.992023 

Рисунок  2.28  -  Таблица  с  апостериорными  вероятностями 


Например,  в  последней  строке  (пациент  №137)  (рис.  2.28) 
максимальное  значение  вероятности  соответствует  группе  3  (норма),  т.е.  с 
вероятностью  0,992023  данный  пациент  относится  к  группе  №3.  Согласно 
таблице  квадратов  расстояний  Махалонобиса  (рис.  2.26)  расстояние  от 
данного  пациента  до  центроидов  групп  минимально  для  группы  3 
(3,93736).  Это  дополнительное  подтверждение  того,  что  исследуемый 
пациент  относится  к  группе  «норма». 

В  случае,  когда  флажок  как  на  рис.  2.23  установлен  в  поле  «8ате  й)г 
аіі  §гоир8»  (одинаковые  для  всех  групп),  каждая  группа  пациентов  имеет 
априорную  вероятность,  равную  0,333.  Следовательно,  число  объектов, 
которые  будут  правильно  предсказаны  при  случайной  классификации  их 
по  группам  пропорционально  априорным  вероятностям,  равно: 


8 

1ІРкпк  =  (0,333-52)  +  (0,333-23)  +  (0,333 -62)  =  45,62. 
к=  1 
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Для  130  правильных  предсказаний  из  137  возможных  т-статистика, 
согласно  формуле  (1.9),  составит: 

130-45,62 

137-45,62 

Это  означает,  что  классификация  с  помощью  дискриминантных 
функций  делает  на  92,3%  ошибок  меньше,  чем  ожидалось  при  случайной 
классификации. 

В  случае  если  флажок  установлен  в  поле  «Ргорогйопаі  Іо  §гоир  зігез» 
(пропорциональные  размерам  групп)  (рис.  2.  2),  каждая  группа  пациентов 
имеет  свою  априорную  вероятность.  Следовательно,  число  объектов, 
которые  будут  правильно  предсказаны  при  случайной  классификации  их 
по  группам  пропорционально  априорным  вероятностям,  равно: 

8 

1ІРкпк  =  (0,379-52)  +  (0,168-23)  +  (0,453 -62)  =  51,66. 

к= 1 

Для  130  правильных  предсказаний  из  137  возможных  т  -статистика 
согласно  формуле  (1.9)  составит: 


_  130-51,66 
1  "137-51,66 


0,918. 


В  случае  если  флажок  установлен  в  поле  «ІІзег  бейпеё»  (рис.  2.22), 
то  априорные  вероятности  задаются  пользователем. 

В  качестве  показателя  качества  алгоритма  диагностики 
(прогностической  «силы»)  использовали  долю  правильной  диагностики  ц, 
определяемую  по  формуле  (1.10). 

рі=52/137=0,379;  р2=23/137=0,168;  р3=62/1 37=0,453. 


8 

X  ркпк  =  0,379  +  0,168  +  0,453  =  1. 
к=\ 
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К- =  —  =  0,981; 
52 


1  О 

Я  =  —  =  0,783; 
23 


іі  =  —  =  0,984. 

62 

/7  =  0, 379  •  0, 981  +  0, 168  •  0, 783  +  0, 453  •  0, 984  =  0, 9489 . 


Для  корректного  применения  дискриминантных  функции  при 
решении  поставленной  задачи  использовались  две  выборки  пациентов,  по 
137  человек  каждая:  одна  для  вычисления  дискриминантной  функции 
(анализируемая),  вторая  -  проверочная  (для  проверки  результатов  расчета 
на  основании  первой  выборки).  Т.е  была  проведена  кросс-проверка 
(перекрестная  проверка). 

В  результате  проведения  клинической  апробации  были  получены 
следующие  результаты  (рис.  2.29).  К  группе  «Положительная 
физнагрузка»  относится  52  наблюдения.  51  наблюдение  (98,07692  %)  было 
определено  корректно,  1  ошибочно  отнесено  к  группе  «Отрицательная 
физнагрузка».  В  группе  «Отрицательная  физнагрузка»  насчитывается  23 
наблюдения.  20  из  них  были  определены  корректно,  2  наблюдения 
ошибочно  отнесено  к  группе  «Положительная  физнагрузка»  и  1  к  группе 
«Норма».  В  группе  «Норма»  насчитывается  62  наблюдения.  Корректно 
определены  61,  ошибочно  -  1  к  группе  «Положительная  физнагрузка». 


Сгоир 

СІаззФсабоп  Маігіх  (ЗргеасізИеегё) 

Коѵѵз:  ОЬзегѵесІ  сІаззФсаІіопз 

Соіитпз:  Ргесіісіесі  сІаззФсабопз 

Регсепі 

Соггесі 

положительная 
физнагрузка 
р=. 37956 

отрицательная 
физнагрузка 
р=. 16788 

норма 

р=,45255 

положительная  физнагрузка 

98.07692 

51 

1 

0 

отрицательная  физнагрузка 

86.95652  2 

20 

і] 

норма 

98.38710  1 

0 

61 

Тоіаі 

96.35036] 54 

21 

62 

Рисунок  2.29  -  Результаты  классификации 
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Анализ  полученных  результатов  показал,  что  прогноз  для 
исследуемых  групп  сделан  верно  (корректно  классифицированы  в  общей 
сложности  96,  35  %  всех  случав  заболеваний  сердца). 

Результаты  оценки  показателя  качества  алгоритма  диагностики 
(прогностической  «силы»)  приведены  ниже. 

рі=52/137=0,379;  р2=23/137=0,168;  р3=62/1 37=0,453. 


8 

X  Ркпк  =0,379  +  0,168  +  0,453  =  1. 

к= 1 


к- =  —  =  0,981; 

52 

Я  =  —  =  0,869; 

23 

77  =  —  =  0,984. 

62 

//  =  0, 379  •  0, 98 1  +  0, 168  •  0, 870  +  0, 453  •  0, 984  =  0, 9635 

Этап  7.  Формирование  диагностического  заключения. 

Таким  образом,  использование  разработанного  метода  позволяет 
выявить  систолическую  дисфункцию  миокарда  у  подростков  и 
своевременными  превентивными  мерами  предотвратить  снижение 
физической  активности  пациента. 

Предметом  дальнейшего  совершенствования  метода  является 
использование  в  диагностирующей  модели  других  клинических 
показателей  в  соответствии  с  новыми  концепциями  диагностики 
систолической  дисфункции  миокарда. 
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ГЛАВА  3.  ИСПОЛЬЗОВАНИЕ  БИНАРНОЙ  ЛОГИСТИЧЕСКОЙ 
РЕГРЕССИИ  ДЛЯ  ОПРЕДЕЛЕНИЯ  СЕРДЕЧНО-СОСУДИСТОГО 
РИСКА  У  ПАЦИЕНТОВ  С  ПОЧЕЧНОЙ  НЕДОСТАТОЧНОСТЬЮ 

В  течение  последних  десятилетий  в  мире  наблюдается  повсеместный 
неуклонный  рост  числа  больных,  страдающих  хронической  почечной 
недостаточностью  (ХПН).  В  частности,  в  Украине  насчитывается  около 
130,8  больных  на  100  тыс.  населения,  а  также  наблюдается  высокая 
смертность  пациентов  с  данной  патологией  (8,9%)  [35].  Это  делает 
проблему  лечения  ХПН  одной  из  центральных  в  современной  нефрологии. 
Значительный  прогресс  медицины,  научно-технические  достижения  в 
области  гемодиализа,  широкое  внедрение  в  клинику  перитонеального 
диализа  и  возможность  трансплантации  почки  создали  реальные 
предпосылки  для  успешного  решения  этой  проблемы  [36]. 

Однако,  если  диагностика  и  лечение  ХПН  достаточно  хорошо 
разработаны,  то  при  ведении  пациентов  с  сочетанной  патологией  часто 
возникают  трудности.  Особое  место  среди  сочетанных  состояний  занимает 
ишемическая  болезнь  сердца  (ИБС),  поскольку  пациенты  с  нарушением 
функции  почек  значительно  чаще  умирают  от  острой  ишемии  миокарда, 
чем  от  собственно  терминальной  почечной  недостаточности.  Не  менее 
20%  всех  случаев  смерти  при  ХПН  обусловлены  инфарктом  миокарда 
(ИМ)  [37].  Своевременная  профилактика  и  лечение  сердечно-сосудистых 
заболеваний  при  ХПН  имеет  важное  значение  для  улучшения  результатов 
диализа.  В  этой  связи  приоритет  приобретает  исследование  факторов 
риска,  влияющих  на  возникновение,  клинические  проявления  и  прогноз 
сердечно-сосудистых  заболеваний  у  пациентов  с  патологией  почек. 

Еще  более  актуальным  представляется  решение  этих  вопросов  в 
свете  современной  демографической  проблемы,  а  именно  в  аспекте  данных 
о  непрерывном  росте  в  последние  15-20  лет  числа  больных,  страдающих 
хронической  болезнью  почек  (ХБП),  и  результатов  популяционных 
исследований,  демонстрирующих  вклад  ХБП  в  кардиоваскулярную 
смертность  населения. 

В  настоящее  время  можно  считать  доказанным,  что  снижение 
функции  почек  является  независимым  фактором  риска  развития  сердечно¬ 
сосудистой  патологии.  Почечная  недостаточность  может  быть  либо 
маркером  пока  неизвестных  факторов  кардиоваскулярного  риска,  либо 
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сама  по  себе  вызывать  и  ускорять  прогрессирование  сердечно-сосудистой 
патологии,  следовательно,  выявление  факторов,  влияющих  на 
прогрессирование  ИБС  у  этой  категории  больных  является  актуальной 
задачей. 

Для  решения  задачи  определения  факторов  риска  развития  сердечно¬ 
сосудистой  патологии  используют  дискриминантный  анализ, 
представленный  в  главах  1  и  2  данной  монографии.  Применение 
дискриминантного  анализа  подтверждается  широким  распространением 
пакетов  прикладных  программ,  включающих  его. 

Несмотря  на  известные  достижения  в  области  классификации  и 
анализа  информативности  данных  средствами  дискриминантного  анализа 
имеется  ряд  сдерживающих  факторов.  Отдельные  измерения  могут  быть 
тесно  связанными  (коррелированными).  Все  это  приводит  к 
вырожденное™  или  плохой  обусловленности  оценок  ковариационных 
матриц  (плохо  обусловленные  задачи).  Плохая  обусловленность  является 
сдерживающим  фактором  в  использовании  методов  дискриминантного 
анализа. 

В  рассматриваемом  нами  случае,  при  изучении  факторов  ,  влияющих 
на  прогрессирование  ИБС,  предпосылка  дискриминантного  анализа 
отсутствия  значимых  различий  между  ковариационными  матрицами  для 
исследуемых  групп  не  выполняется. 

Подход  на  основе  логистических  моделей  обладает  тем 
преимуществом,  по  сравнению  с  дискриминантным  анализом,  что 
обладает  гораздо  менее  жесткими  ограничениями,  и  поэтому  имеет  более 
широкую  область  применения. 

Вероятностная  природа  причинно-следственных  связей  в  развитии 
ИБС  определяет  концепцию  многофакторности.  Научные  концепции, 
основанные  на  анализе  факторов  риска  в  развитии,  прогрессировании  и 
исходе  патологических  состояний  и  заболеваний  имеют  практическое 
подтверждение  при  изучении  заболеваний  сердечно-сосудистой  системы. 

Одним  из  адекватных  методов  анализа  подобных  связей  является 
метод  логистической  регрессии,  который  в  отечественной  биомедицине 
начал  активно  применяться  в  последние  2-3  года.  Между  тем  в  зарубежной 
практике  это  довольно  популярный  метод.  Так,  например,  с  помощью 
логистических  моделей  изучали  прогноз  коронарной  болезни  сердца  и 
атеросклероза  у  взрослого  населения  [38]. 
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Целью  исследования  являлось  изучение  факторов  риска  сердечно¬ 
сосудистых  заболеваний,  выделение  среди  них  наиболее  значимых, 
оказывающих  наиболее  значительное  влияние  на  прогрессирование 
хронической  ИБС  в  виде  ишемической  делятационной  кардиомиопатии,  с 
использованием  бинарной  логистической  регрессии. 

Исследование  проводилось  в  два  этапа.  На  первом  этапе  при  помощи 
бинарной  логистической  функции  моделировалась  зависимость  наличия 
или  отсутствия  ишемической  дилятационной  кардиомиопатии  у 
пациентов,  находящихся  на  перитонеальном  диализе  от  наличия  ряда 
клинических  признаков.  На  втором  этапе  производилось  тестирование 
созданной  модели  с  позиций  ее  статистической  значимости. 

Все  пациенты  были  разделены  на  группы  следующим  образом:  1-я 
группа  -  пациенты  без  ишемии  миокарда;  2-я  группа  -  пациенты  с 
ишемической  дилятационной  кардиомиопатией. 

В  качестве  потенциальных  предикторов  использовали 

гемодинамические  факторы  (среднее  артериальное  давление,  гемоглобин, 
гематокрит  и  др.);  нутриционный  статус  (альбумин,  индекс  массы  тела, 
площадь  поверхности  тела  и  др);  фосфорно-кальциевый  обмен  (фосфор, 
кальций,  интактный  паратиреоидный  гормон  и  др.);  дислипидемия 
(липопротеиды  высокой  плотности,  липопротеиды  низкой  плотности  и  др); 
показатели  адекватности  перитонеального  диализа  (средние  молекулы, 
скорость  клубочковой  фильтрации  и  др);  системное  воспаление 
(интерлекины- 1  (3,  интерлейкин- 8,  С-реактивный  белок,  сывороточный 
амилоид  А,  фибриноген  ,  ферритин,  (32-  микроглобулин  и  др.). 

Все  показатели  были  закодированы  и  поставлены  соответственно  87- 
мерному  вектору,  который  учитывает  отсутствие,  наличие,  направленность 
и  величину  каждого  показателя. 

Математическая  обработка  результатов  проводилась  с 
использованием  и  пакета  прикладных  программ  8Р88  21  для  \Ѵі псіохѵз. 

Для  ввода  всех  закодированных  данных  в  компьютер,  после  запуска 
8Р88  для  \УіпсІо\У8  с  помощью  двойного  нажатия  левой  кнопкой  мыши  на 
значке  8Р88  открыли  редактор  данных  и  внесли  всю  необходимую 
информацию  (рис.  3.1). 
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Рисунок  3.1  -  Диалоговое  окно  8Р88.  Внесенная  информация 


Выбрали  в  меню  «Анализ»—»  «Регрессия»  — >  «Логистическая» 
(рис.  3.2).  Открылось  диалоговое  окно  «Логистическая  регрессия» 
(рис.  3.3).  Поместили  единственную  группирующую  переменную 
«группы»  (сравниваемые  группы),  содержащую  информацию  о 
принадлежности  к  одному  или  второму  классу  (нет  ИБС  или  ИДК),  в  поле 
для  зависимых  переменных,  а  все  необходимые  переменные  —  в  поле 
ковариат. 

В  качестве  метода  использования  переменных  в  вычислениях 
предварительно  бал  установлен  метод  «Принудительное  включение»,  при 
котором  в  расчёт  одновременно  вовлекаются  все  переменные  объявленные 
ковариатами.  Однако  количество  объявленных  ковариат  должно  быть 
меньше  чем  количество  наблюдений.  Альтернативой  здесь  являются 
прямая  и  обратная  селекции  (пошаговые  методы).  Метод  прямой  селекции 
начинается  с  использования  одних  лишь  констант  на  стартовом  этапе,  а 
затем  последовательно  подключаются  переменные,  которые 
демонстрируют  сильную  корреляцию  с  зависимыми  переменными.  Далее 
опять  следует  проверка  того,  какие  переменные  должны  быть  исключены, 
причём  в  качестве  критерия  проверки  выбирается  либо  статистика  \Уа1ё 
(«Включение:  Вальд»),  либо  функция  правдоподобия  («Включение:  ОП»), 
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либо  один  из  вариантов,  называемых  «условной  статистикой» 
(«Включение:  Условное»  -  которые,  однако,  не  рекомендуются)  [39]. 


Файл  Правка  Вид  Данные  Преобразовать 


Анализ  Графика  Сервис  Дополнения  Окно  Справка 


&0І  Е  Ъ  Г*  1.ІІ  #4  * 


1  :  Группы 


18 


1,0 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

Группы 


Поп 


Данные  Переменные 


Возраст 


Логистическая 


Отчеты  ► 

Описательные  статистики  ► 

Таблицы  ► 

Анализ  КРМ  > 

Сравнение  средних  > 

Общая  линейная  модель  > 

Обобщенные  линейные  модели  ► 
Смешанные  модели  ► 

Корреляции  > 

Реірессия  > 

Логлинейный  ► 

Нейронные  сети  ► 

Классификация  > 

Снижение  размерности  ► 

Шкалирование  ► 

Непараметрические  критерии  ► 

Прогнозирование  > 

Выживаемость  > 

Множест  венныеответы  ► 

И  Анализ  пропущенных  значений... 
Множественная  импутация  ► 

Сложные  выборки  ► 

Контроль  качества  > 

0  КОС  кривые... 


Показать  переменные:  1 1 6  из  1 1 6 


СКФ.Кач 

СКФ1 

ППТ 

,40 

1 

0,00 

2,03 

,80  1 

4,90 

1,75 

,40  1 

2,10 

1,74 

,50  1 

0,90 

1,84 

,10  1 

0,00 

1,99 

Л 


ппп 


Линейная... 

Подгонка  кривых... 

Метод  частичных  наименьших  квадратов .. 


Логистическая... 

Мультиномиальная  логистическая.. 
Порядковая... 

Пробит... 


Нелинейная... 

Взвешенная... 

Двух  этапный  МН  К... 


Категориальная.. 


ППТ1 


5Р55  Зіаііеіісз  Ргосессог  іс  геасіу 


Рисунок  3.2  -  Диалоговое  окно  8Р88.  Выбор  метода  анализа 


Метод  обратной  селекции  сначала  берёт  в  расчёт  все  переменные,  а 
затем  в  обратном  порядке  происходит  исключение  малозначимых 
переменных.  Мы  же  воспользовались  пошаговым  методом  «Включение 
ОП»  (рис.  3.3),  при  котором  в  качестве  критерия  проверки 
последовательно  включенных  переменных  выбирается  функция 
правдоподобия. 

Используя  кнопку  «Категориальные»  подготовили  для  расчета 
категориальные  переменные  (то  есть  переменные,  принадлежащие  к 
номинальной  шкале).  На  рис.  3.4  представлено  диалоговое  окно  8Р88, 
которое  открывается  после  нажатия  на  кнопку  «Категориальные»  для 
задания  категориальных  переменных. 
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Рисунок  3.3  -  Диалоговое  окно  8Р88.  Логистическая  регрессия 


После  установки  всех  необходимых  настроек  начали  расчёт 
нажатием  кнопки  «ОК»  (рис.  3.3).  Результаты  проведенного  расчета 
приведены  ниже. 


Рисунок  3.4  -  Диалоговое  окно  8Р88.  Логистическая  регрессия: 
Задать  категориальные  переменные 
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С  учетом  рассмотренных  показателей  составили  уравнение 
логистической  регрессии,  по  которому  определяется  вероятность  наличия 
ишемической  дилятационной  кардиомиопатии  у  пациентов,  находящихся 
на  перитонеальном  диализе: 


Р  =  [1  +  ехр(г)]_1 , 


где  Р  -  вероятность  того,  что  произойдет  интересующее  событие; 

2  -  значение  логистической  функции,  которое  определяется 
согласно  формуле: 


2=Ь1*Х1  +  ъ2*х2+...+  ъп*хп  +ь0, 

где  Хі  -  фактор  или  предиктор, 

Ь;  -  коэффициенты,  значения  которых  указывают  на  изменение 
отношения  шансов  при  изменении  Х;. 

Отношением  шансов  является  отношение  вероятности  того,  что 
событие  произойдет  к  вероятности  того,  что  оно  не  произойдет, 
высчитывается  как  Р/(1-Р). 

Параметры  модели  оценивались  с  помощью  метода  максимального 
правдоподобия  [40]. 

Метод  пошагового  включения  предикторов,  который  ранжирует 
признаки  в  соответствии  с  их  вкладом  в  модель,  применялся  при  оценке 
уравнений  регрессии.  Итогом  стало  построение  регрессионной  функции,  в 
которую  вошли  4  показателя:  -  площадь  поверхности  тела  (в  м2);  Х2  - 

сывороточный  амилоид  А  (в  мг/л);  Х3  -  интерлейкин  8  (в  пг/мл);  Х4  - 
липопротеиды  высокой  плотности  (в  ммоль/л). 

Подобранная  из  обучающей  выборки  бинарная  логистическая 
функция  имеет  вид: 


Р  =  [Г +  ехр(68,739  •  X,  +  77,496  •  Х2  -  0,038  •  Х3  - 10,257  •  Х4  - 134,818)] 1 , 


где  Р  -  вероятность  появления  ишемической  дилятационной 
кардиомиопатии  у  пациентов,  находящихся  на  перитонеальном  диализе. 
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Таким  образом,  сравнение  групп  1  и  2  показало,  что  из  всех 
оцениваемых  факторов  статистически  значимое  влияние  на  вероятность 
развития  ишемической  дилятационной  кардиомиопатии  у  пациентов, 
находящихся  на  перитонеальном  диализе  оказывали  липопротеиды 
высокой  плотности,  площадь  поверхности  тела,  интерлейкин-8, 
сывороточный  амилоид  А. 

Проверка  значимости  подобранных  коэффициентов  Ь1  проводилась 
при  помощи  статистики  Вальда,  представляющей  собой  квадрат 
отношения  соответствующего  коэффициента  к  его  стандартной 
ошибке  [41]: 


Все  переменные,  согласно  статистики  Вальда  (табл.  3.1),  значимы 
(р<0,05)  и  подобраны  правильно. 

Рассчитанные  коэффициенты  регрессионной  функции  и  результаты 
проверки  их  значимости  представлены  в  табл.  3.1. 


Таблица  3.1 

Коэффициенты  регрессионной  функции 


Показа¬ 

тели 

(Хі) 

Коэфи- 

циенты(Ьі) 

Стандартные 
ошибки  (8і) 

Критерии 

Вальда 

ОѴі) 

Значимость 

(Рі) 

Критерий 

0(М§  Кайо 

95%  Доверительный 

интервал  для 

0(М8  Кайо 

Нижняя 

Верхняя 

X, 

68,739 

28,110 

5,980 

0,014 

7Д27Е29 

841956,444 

6,032Е53 

х2 

77,496 

33,777 

5,264 

0,022 

4,531ЕЗЗ 

80338,851 

2,556Е62 

Х3 

-0,038 

0,016 

6,020 

0,014 

0,963 

,934 

,992 

Х4 

-10,257 

5,049 

4,127 

0,042 

0,000 

,000 

,697 

Константа 

-134,818 

55,353 

5,932 

0,015 

0,000 
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Полученные  результаты  свидетельствуют  о  вероятности  перехода  от 
состояния  «без  ишемии»  к  «ишемической  дилятационной 
кардиомиопатии» . 

Столбец  Ъі  показывает  изменение  вероятности  отнесения  пациента  к 
группе  с  ишемической  дилятационной  кардиомиопатией  (в  %)  при 
изменении  каждого  фактора  на  единицу. 

Кроме  того,  в  таблице  представлен  критерий  Ос1с1$  Кайо  -  отношение 
шансов  «без  ишемии  миокарда»  и  «ишемическая  дилятационная 
кардиомиопатия»,  который  рассчитывается  согласно  формуле: 

ОсісИ  Пайо  =  еЬ/ . 

Вывод  доверительных  интервалов  для  ОйсВ  Кайо  осуществляется 
путем  нажатием  кнопки  «Опции»  в  диалоговом  окне  «Логистическая 
регрессия»  (рис.  3.3)  и  установкой  флажка  в  новом  окне  (рис.  3.5)  напротив 
поля  «ДИ  для  ехр(В)». 


Рисунок  3.5  -  Диалоговое  окно  8Р88. 
Логистическая  регрессия:  Параметры 
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В  области  «Критерии  шагового  отбора»  имеются  поля  «Включение» 
и  «Исключение».  В  первом  указывается  значение  вероятности,  служащее 
критерием  включения  предиктора  в  уравнение  регрессии,  а  во  втором  — 
значение  вероятности  для  исключения  из  уравнения  регрессии.  По 
умолчанию  в  поле  «Включение»  указано  значение  0,05,  а  в  поле 
«Исключение»  —  значение  0,10.  При  желании  исследователь  может 
изменять  их. 

Созданная  регрессионная  модель  имеет  высокое  качество  и  для 
97,7%  объектов  результаты  прогноза  оказались  верными  (табл.  3.2  и  3.3). 

Таблица  3.2 

Характеристики  модели  бинарной  логистической  регрессии, 
созданной  для  определения  риска  ишемической  дилятационной 
кардиомиопатии  у  пациентов,  находящихся  на  перитонеальном  диализе 


Результаты 

заключительного 

шага  анализа 

-2  Ьо§ 

Правдоподобие 

(О) 

К 2  Кокса  и 

Шелла 

К2 

Нэйджелкерка 

х2 

Р 

10,940 

0,665 

0,944 

95,180 

0,001 

Классификационные  результаты 


Таблица  3.3 


Наблюдаемые 

группы 

Спрогнозированные  результаты  (и) 

%  верно  предсказанных 

показателей 

Без  ишемии 

миокарда 

Ишемическая 

дилятационная 

кардиомиопатия 

Без  ишемии 

25 

1 

96,2 

миокарда 

Ишемическая 

1 

60 

98,4 

дилятационная 

кардиомиопатия 

Суммарный  показатель 

97,7 

Качество  приближения  регрессионной  модели  оценивается  при 
помощи  функции  подобия  [42].  Мерой  правдоподобия  служит 
отрицательное  удвоенное  значение  логарифма  этой  функции  (-2ІХ).  В 
качестве  начального  значения  для  -2ІХ  применяется  значение,  которое 
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получается  для  регрессионной  модели,  содержащей  только  константы. 
Значение  -2  Ьо§  Правдоподобия  определяется  как: 


0  =  ~2  (ЬЬпи11-ЬЬк\ 


где  ^пиіі  _  логарифм  правдоподобия  модели,  состоящей  только  из 
констант  [43  -  44]. 

ііегсерг  -  логарифм  прав допо добия  модели,  состоящей  из  к 


предикторов.  При  нулевой  гипотезе,  Н0:Ъі=Ъ2=. .  .=0  статистика  О  имеет  хи- 
квадрат  распределение  и  оценивает,  насколько  хорошо  независимые 
факторы  влияют  на  переменную.  Полученное  в  исследовании  значение  О  = 
10,940  при  Р=0,001  показывает,  что  в  целом  независимые  переменные 
имеют  значительный  вклад,  чтобы  спрогнозировать  зависимую 
переменную  [45]. 

Показатели  Кокса  и  Шелла,  Нэйджелкерка,  являются  мерами 
определённости.  Они  указывают  на  ту  часть  дисперсии,  которую  можно 
объяснить  с  помощью  логистической  регрессии.  Мера  определённости  по 
Коксу  и  Шеллу  имеет  тот  недостаток,  что  значение  равное  1  является 
теоретически  не  достижимым;  этот  недостаток  устранен  благодаря 
модификации  данной  меры  по  методу  Наделькеркеса.  Часть  дисперсии, 
объяснимой  с  помощью  логистической  регрессии,  в  рассматриваемом 
нами  случае  составляет  94,4  %. 

Формула  Кокса  и  Шелла  заключается  в  следующем  [46]: 


Показатель  Нэйджелкерка  находится  как: 


1  -  (-2 Ыпи11)У" 
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На  рис.  3.6  представлена  диаграмма  классификации.  В  пакете  8Р88 
вывод  диаграмм  классификации  осуществляется  установкой  флажка  в 
диалоговом  окне  «Логистическая  регрессия:  Параметры»  (рис.  3.5) 
напротив  поля  «Графики  классификации». 


к| 

к+ 

К| 

К| 

К| 

к+ 

К| 

К| 

К| 

к+ 

К| 

К| 

К| 

О  0,1  0,2  0,3  0,4  0,5  0,6  0,7  0,8  0,9  1 

Предсказанная  вероятность 

Рисунок  3.6  -  Диаграмма  классификации. 

Символы:  и  -  без  ишемии  миокарда;  к  -  дилятационная 
кардиомиопатия 


Точкой  разделения  служило  значение  р=0,5. 

По  горизонтальной  оси  отложены  значения  предсказанной 
вероятности,  вычисляемые  по  уравнению  регрессии,  по  вертикали  - 
частоты.  Таким  образом,  каждый  столбик  на  диаграмме  соответствует 
определенной  предсказанной  вероятности,  а  его  высота  —  количеству 
объектов,  для  которых  предсказана  данная  вероятность.  Чем  ближе 
значение  прогнозируемой  вероятности  к  единице,  тем  вероятнее  развитие 
дилятационной  кардиомиопатии  у  пациента,  находящегося  на 
перитонеальном  диализе. 

О  количестве  правильных  и  неправильных  предсказаний  позволяет 
судить  классификационная  табл.  3.3. 

Из  таблицы  можно  сделать  вывод  о  том,  что  из  общего  числа 


80  + 


60  н 


(О  I 

I 


20  + 
|н 

|Н 

|н 


79 


больных  с  ишемической  дилятационной  кардиомиопатией,  равного  61, 
тестом  были  признаны  верно  60  и  1  ошибочно  отнесен  к  группе  «Без 
ишемии  миокарда».  Из  общего  числа  здоровых  («Без  ишемии  миокарда»), 
равного  26,  тестом  были  признаны  верно  25  и  1  ошибочно  отнесен  к 
группе  больных.  В  общем,  правильно  были  распознаны  85  случаев  из  87, 
это  составляет  97,7  %. 

Общая  оценка  согласия  между  влиянием  выявленных  в  модели 
факторов  риска  и  реально  зафиксированным  наступлением 
неблагоприятного  исхода  производилась  с  использованием  теста  согласия 
Хосмера-Лемешова  ( Нь ),  в  котором  значение  р  тем  выше,  чем  меньше 
различия  между  частотой  наблюдаемых  и  предсказанных  на  основании 
данных  регрессионной  модели  исходов. 

Критерий  Хосмера-Лемешова  [47  -  48],  который  показывает  степень 
различия  между  оцененными  и  фактическими  значениями  зависимой 
переменной  в  обеих  группах  (значения  бинарной  переменной  0  и  1) 
определяется  как: 


8 


7=1 


где  п]  -  количество  наблюдений  в  /-й  группе; 

°.і=Цу«  =  наблюдаемое  число  случаев  в  /-й  группе; 

і 

А 

=  ожидаемое  число  случаев  в  ]- й  группе. 

і 

Полученное  наблюдаемое  О і  и  ожидаемое  Еі  число  случаев  в  /-й 

группе  для  каждого  шага  представлено  в  табл.  3.2. 

Согласно  полученным  <9уи  Еі  (табл.  В.1,  Приложение  В)  для 

каждого  шага  были  рассчитаны  значения  критерия  Хосмера-Лемешова 
(табл.  3.4).  В  8Р88  расчет  теста  согласия  Хосмера-Лемешова 
осуществляется  путем  установки  опции  «Критерий  согласия  Хосмера- 
Лемешова»  в  диалоговом  окне  8Р88.  Логистическая  регрессия:  Параметры 
(рис.  3.5). 
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Критерий  Хосмера-Лемешова 


Таблица  3.4 


Шаг 

Степень 

свободы 

Значимость 

Р, 

1 

10,066 

7 

0,185 

2 

10,564 

8 

0,228 

3 

0,254 

5 

0,998 

4 

1,227 

4 

0,874 

Полученное  на  последнем  шаге  низкое  значение  Нь  =1,227,  при 
уровне  значимости  Р> 0,05  (р=0,874),  свидетельствует  о  приемлемом 
качестве  подобранной  модели. 

Вывод  остатков  (рис.  3.7)  осуществляется  путем  установки  в 
диалоговом  окне  «Логистическая  регрессия:  Параметры»  опции 
«Поточечный  вывод  остатков»  (рис.  3.5). 


Рисунок  3.7  -  Поточечный  вывод  остатков 
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В  первом  столбце  таблицы  (рис.  3.7)  указан  номер  исследуемого 
пациента,  во  втором  столбце  буквой  О  обозначены  отобранные 
наблюдения,  *  *  -  неверно  классифицированные  наблюдения.  Таким 
образом,  пациент  №16  относящийся  к  группе  «Без  ишемии  миокарда» 
ошибочно  отнесен  к  группе  «Ишемическая  дилятационная 
кардиомиопатия»,  а  пациент  под  номером  №43,  относящейся  к  группе 
«Ишемическая  дилятационная  кардиомиопатия»,  ошибочно  отнесен  к 
группе  «Без  ишемии  миокарда». 

В  пакете  8Р88  в  диалоговом  окне  «Логистическая  регрессия: 
Параметры»  (рис.  3.5)  флажок  «Корреляции  оценок»  указывает  на 
необходимость  построения  корреляционной  матрицы  для  всех 
переменных,  включенных  в  регрессионное  уравнение.  Полученная 
корреляция  оценок  представлена  на  рис.  3.8. 


Рисунок  3.8  -  Корреляция  оценок 


Также  в  пакете  8Р88  есть  возможность  просмотра  истории  итераций 
(рис.  В.1,  Приложение  В)  путем  установки  опции  «История  итераций»  в 
диалоговом  окне  «Логистическая  регрессия:  Параметры»  (рис.  3.5). 
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Для  оценки  эффективности  модели  использовался  КОС-анализ.  Для 
проведения  КОС-анализа  на  первом  этапе  получили  дополнительные 
переменные.  При  помощи  кнопки  «Сохранить»  в  диалоговом  окне  8Р88. 
«Логистическая  регрессия»  (рис.  3.3)  перешли  в  окно  «Логистическая 
регрессия:  Сохранить»  (рис.  3.9)  и  установили  флажки  «Вероятности»  и 
«Принадлежность  к  группе».  После  чего  в  файл  добавились  новые 
переменные  РКЕ_1  (предсказанные  вероятности)  и  РОК_1  (предсказанные 
принадлежности  в  группе). 


Рисунок  3.9  -  Диалоговое  окно  8Р88.  «Логистическая  регрессия: 

Сохранить» 

Далее  в  меню  выбрали  «Анализ»  — >  «РОС-кривые»  как  показано  на 
рис.  3.10.  В  открывшемся  диалоговом  окне  «КОС-кривые»  (рис.  3.11) 
выбрали  в  качестве  проверяемой  переменной  полученную  переменную 
«РОК1»  (предсказанные  принадлежности  к  группе),  а  в  качестве 
переменной  состояния  -  переменную  «группы»  (фактическая 
принадлежность  к  группе).  Значение  переменной  состояния  установили  2, 
т.к.  2  -  прогнозируемая  группа  (ишемическая  дилятационная 

кардиомиопатия).  После  установки  всех  необходимых  настроек  нажали 
кнопку  «ОК». 
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Файл  Правка  Вид  Данные  Преобразовать 

Анализ  Графика  Сервис  Дополнения 

Окно  Справка 

&  0  й 

К  ГГІ  §9 

ПШД  §§• 

#4  * 

Отчеты  ► 

22 :  Группы 

2 

Описательные  статистики  ► 

Показать  переменные:  1 1 8  из  1 

N 

Группы 

Таблицы  ► 

Мес  ПД 

Диагноз 

С 

і 

39 

1 

Бабаев 

Анализ  ВРМ  ► 

і 

60 

2 

2 

40 

1 

Виниче 

Сравнение  средних  ► 

4 

2 

3 

41 

1 

Ковале 

Об  шея  линейная  модель  ► 

18 

8 

4 

42 

1  Кобец^ 

Обобщенные  линейные  модели  ► 

28 

2 

5 

43 

1 

Когтенн 

Смешанные  модели  ► 

9 

38 

2 

6 

44 

1 

Новоха 

Корреляции  ► 

5 

60 

2 

7 

45 

1 

Рыбалн 

Регрессия  ► 

5 

43 

5 

8 

46 

1  Довжег 

Логлинейный  ► 

4 

14 

2 

9 

47 

1  Духани 

Нейронные  сети  ► 

1 

52 

2 

10 

48 

1  Денисе 

Классификация  ► 

1 

20 

8 

11 

49 

1  Гаврил 

Снижение  размерности  ► 

2 

24 

2 

12 

50 

1  Грицен 

Шкалирование  ► 

2 

42 

8 

13 

51 

1 

Камене 

Непараметрические  критерии  ► 

2 

30 

5 

14 

52 

1 

Плохот 

Протозирование  ► 

7 

41 

8 

15 

53 

1 

Панюш 

Выживаемость  ► 

4 

67 

2 

16 

54 

1  Троцен 

Множест  в  енные  от  в  еты  ► 

9 

21 

8 

17 

55 

2  Асрян 

Щэ)  Анализ  пропущенных  значений... 

4 

14 

9 

18 

56 

2  Бондас 

Множест  в  енна  я  импу  таци  я  ► 

3 

61 

1 

< 

Сложные  выборки  ► 

I* 

Данные  Переменные 

Контроль  качества  ► 

РОС  кривые... 

^  РОС  кривые... 

іісз  Ргосеезог  іе  геасіу 

Рисунок  3.10  -  Диалоговое  окно  8Р88.  Выбор  метода 


Рисунок  3.11  -  Диалоговое  окно  8Р88.  КОС-кривые 


Результаты  КОС-анализа  приведены  в  табл.  3.5  и  на  рис.  3.12. 
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Таблица  3.5 

Результаты  КОС-анализа 


Характеристики  КОС  Кривой 

95%  Доверительный  интервал 

Площадь 

Стандартная 

ошибка 

Значимость  Р 

Нижняя 

граница 

Верхняя  граница 

0,974 

0,023 

0,0001 

0,0001 

1,000 

і 

о 

Л 

к 

о 

О 

5° 

с 

ф 

ь~ 

X 

а 

н  о 
о 

т 

Э- 

Т 

О 


Значение  площади  под  кривой  АІТС  (Агеа  Тіпсіег  Сигѵе)  составило 

0,974. 

Кос-анализ  полученной  модели  (рис.  3.12,  табл.  3.5)  выявил  ее 
характеристики,  являющиеся  показателем  отличного  качества  исследуемой 
модели. 

Таким  образом,  у  пациентов,  находящихся  на  перитонеальном 
диализе,  большая  площадь  поверхности  тела,  низкое  содержание 
липопротеидов  высокой  плотности  может  служить  фактором  риска 
развития  и  прогрессирования  хронической  ишемической  болезни  сердца  в 
виде  ишемической  дилятационной  кардиомиопатии  со  стойкой 
эндотелиальной  дисфункцией  коронарных  артерий  и  диастолической 
дисфункцией  с  исходом  в  диффузный  кардиосклероз  и  сердечную 
и  ед  о  статочно  сть . 

Применение  синтезированной  логистической  модели  позволит 
проводить  своевременную  диагностику,  лечение  и  профилактику 
коронарных  событий  у  больных  с  почечной  недостаточностью. 


0.0  0.2  0.4  0.6  0.8  1, 

Специфичность 

Рисунок  3.12  -  КОС  Кривая 
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ГЛАВА  4.  ПРОГНОЗИРОВАНИЕ  ЗАБОЛЕВАЕМОСТИ 
РАССТРОЙСТВАМИ  ПСИХИКИ  И  ПОВЕДЕНИЯ 
НЕПСИХОТИЧЕСКОГО  ХАРАКТЕРА 

Из-за  социально-экономического  кризиса  претерпевают  быстрые, 
резкие  изменения  условия  существования  людей.  Растёт  безработица, 
меняется  общественный  менталитет,  образ  жизни,  характер  питания. 
Качество  жизни  и  социальное  функционирование  становятся  одним  из 
важнейших  показателей  психического  здоровья.  Термин  «психическое 
здоровье»,  введенный  Всемирной  организацией  здравоохранения, 
обозначает  успешное  выполнение  психических  функций,  имеющих 
результатом  продуктивную  деятельность,  установление  отношений  с 
другими  людьми  и  способность  адаптироваться  к  изменениям  и 
справляться  с  неприятностями.  Психическое  здоровье  -  это  в  первую 
очередь  отсутствие  психических  расстройств  и  форм  инвалидности.  От 
психических  расстройств  страдает  более  450  миллионов  человек.  Помимо 
этого  многие  люди  имеют  разнообразные  психические  проблемы  [49  -  50]. 

Исследования  последних  лет  показали,  что  у  сельских  жителей  по 
частоте  инвалидизации  психические  расстройства  стоят  среди  первых  трех 
категорий  наряду  с  сердечно-сосудистыми  заболеваниями  и  травмами. 

За  период  с  2000  по  2010  гг.  отмечалось  увеличение  показателя 
общей  заболеваемости  психическими  и  поведенческими  расстройствами  у 
сельского  населения  Украины.  При  этом  в  структуре  общей 
заболеваемости  первое  ранговое  место  занимают  непсихотические 
психические  расстройства  [51]. 

Весомая  доля  в  структуре  психических  расстройств  и  отмечаемый 
рост  заболеваемости  непсихотических  психических  расстройств 
указывают  на  актуальность  изучения  тенденций  психического  здоровья 
сельского  населения  Украины  и  составления  прогнозов  заболеваемости, 
необходимых  для  планирования  ресурсов  и  объемов  медицинской  помощи, 
последующей  разработки  научно-практического  комплекса  медико¬ 
социальных  мероприятий,  направленных  на  оптимизацию 
психиатрической  помощи  в  Украине,  в  частности,  сельскому  населению. 

Известно  множество  методов  прогнозирования,  но  все  они  имеют 
одну  общую  идею:  обнаружение  связей  между  прошлым  и  будущим, 
между  информацией  о  процессе  в  контролируемый  период  времени  и 
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характером  протекания  процесса  в  дальнейшем.  От  того,  насколько  точно 
описаны  исследуемые  связи,  будет  зависеть  точность  прогнозирования  [52]. 

В  качестве  прогнозной  модели  нами  была  выбрана  регрессионная 
модель  АКМА.  АКМА-процессы  -  это  класс  стохастических  процессов, 
используемых  для  анализа  временных  рядов.  Модели  позволяют  получать 
точные  прогнозы,  опираясь  только  на  информацию,  содержащуюся  в 
предыстории  прогнозируемых  рядов  медицинских  данных. 

В  практических  медицинских  задачах  можно  было  бы  обойтись  без 
использования  моделей  АКМА,  ограничиваясь  либо  АК  либо  МА 
моделями.  Однако,  при  таком  подходе,  количество  коэффициентов, 
подлежащих  оцениванию,  может  оказаться  слишком  большим  (что 
снижает  точность  оценивания)  и  даже  превосходить  количество 
имеющихся  наблюдений.  В  этом  смысле  модели  АКМА  являются  «более 
экономными». 

Анализ  научных  публикации  последнего  десятилетия  показал,  что  в 
большом  количестве  работ  по  прогнозированию  классическими  методами 
используется  именно  АКМА,  как  один  из  наиболее  обоснованных  и 
надежных  алгоритмов  (из  статистических). 

Прогнозирование  заболеваемости  непсихотическими  психическими 
расстройствами  сельского  населения  различных  регионов  Украины 
позволит  усовершенствовать  систему  оказания  психиатрической  помощи 
сельскому  населению  Украины. 

Метод  прогнозирования  заболеваемости  сельского  населения 
Украины  непсихотическими  психическими  расстройствами  включает  в 
себя  ряд  этапов. 

Этап  1.  Сбор  статистической  информации,  характеризующей 
заболеваемость  непсихотическими  психическими  расстройствами 
сельского  населения  Украины  за  несколько  лет  с  шагом  в  полгода. 

Для  построения  модели  прогнозирования  заболеваемости 
расстройствами  психики  и  поведения  непсихотического  характера  среди 
сельского  населения  Украины,  проживающего  в  24  областях  страны, 
использовали  статистические  данные  за  2008  -2012  гг.  Исходные  данные 
сформировали  согласно  регионального  принципа  разделения:  первый 
регион  -  западноукраинский  (Волынская,  Ровенская,  Львовская, 
Тернопольская,  Ивано  -  Франковская,  Закарпатская,  Черновицкая 
области);  второй  регион  -  центрально-южный  (Киевская,  Винницкая, 
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Кировоградская,  Одесская,  Николаевская,  Херсонская  области  и  АР 
Крым);  третий  регион  -  центрально-северный  (Житомирская, 
Хмельницкая,  Черкасская,  Полтавская,  Сумская,  Черниговская  области); 
четвертые  регион  -  восточный  (Харьковская,  Донецкая,  Луганская, 
Днепропетровская,  Запорожская  области). 

Этап  2.  Предварительный  анализ  данных. 

На  данном  этапе  проводится  структурирование  информации,  расчет 
показателей,  характеризующих  изменение  временного  ряда,  выявление 
аномальных  отклонений. 

Сформированные  данные  о  количестве  случаев  заболеваемости 
расстройствами  психики  и  поведения  представляют  собой  временные  ряды 
(табл.  4.1). 

Таблица  4.1 

Заболеваемость  расстройствами  психики  и  поведения 
непсихотического  характера  сельского  населения  в  4-х  регионах  Украины 


за  2008-2012  года 


Период 

Регионы  (количество  человек} 

і 

1-й  регион 

2-й  регион 

3-й  регион 

4-й  регион 

с  01.01.2008  по 
30.06.2008 

6585 

7222 

5305 

2588 

с  01.07.2008  по 
31.12.2008 

6362 

6891 

5302,5 

2625,5 

с  01. 01. 2009  по 
по  30.06.2009 

6139 

6560 

5300 

2663 

с  01. 07.2009  по 
31.12.2009 

6303,5 

6685,5 

5151 

2604,5 

с  01.01.2010  по 
30.06.2010 

6468 

6811 

5002 

2546 

с  01.07.2010  по 
31.12.2010 

7398,5 

7784,5 

6170,5 

3268 

с  01.01.2011  по 
по  30.06.2011 

8329 

8758 

7339 

3990 

с  01.07.2011  по 
31.12.2011 

8161 

8720,5 

7209 

3935 

с  01.01.2012  по 
по  30.06.2012 

7993 

8683 

7079 

3880 

с  01.07.2012  по 
31.12.2012 

7825 

8646 

6949 

3825 
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По  данным,  представленным  на  графиках  (рис.  4.1),  можно  сделать 
вывод,  что  заболеваемость  непсихотическим  психическими 
расстройствами  у  сельского  населения  на  протяжении  2011  -  2012  г.г. 
носила  отрицательный  характер,  то  есть  наметился  спад  уровня 
заболеваемости  по  всем  регионам  Украины  ,  что  коррелирует  с 
изменением  уровня  инфляции  в  этот  период. 


а)  б) 


4200 
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3600 

3400 

3200 

3000 

2800 
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2400 


Дата 


Дата 
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Рисунок  4.1-  Распределение  заболеваемости  непсихотическими 
психическими  расстройствами  сельского  населения  за  2008-2012  г.  по  4-м 
регионам  Украины:  а)  -  по  1-му  региону;  б)  -  по  2-му  региону; 
в)  -  по  3-му  региону;  г)  -  по  4-му  региону 

Рост  заболеваемости  в  2008  -  201 1  гг.  и  снижение  в  2012  г.  повторяет 
тенденции,  характерные  для  страны  в  целом,  а  именно:  низкий  уровень 
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жизни,  высокий  уровень  безработицы,  высокая  миграционная  активность 
населения,  низкий  уровень  рождаемости. 

Рассмотрим  этот  факт  подробнее. 

Большой  вред  Украине  нанес  мировой  финансовый  кризис. 

На  рис  4.2.  представлен  график  изменения  индекса  инфляции  в 
Украине,  начиная  с  2001  года  нарастающим  итогом. 

Согласно  данным  Государственного  комитета  статистики  Украины 
[53]  количество  официально  зарегистрированных  безработных  в  2011  году 
увеличилось  на  5,3%  по  сравнению  с  2010  годом. 

Однако  в  2012  году  безработица  в  Украине  постепенно  начала 
снижаться.  Уровень  зарегистрированной  безработицы  в  целом  по  Украине 
в  конце  2011  года  составил  1,5%  от  общего  количества  трудоспособного 
населения,  что  на  0,  7  %  меньше,  чем  1  марта  2011  года. 


Рисунок  4.2  -  График  изменения  индекса  инфляции  в  Украине 

При  этом,  следует  отметить,  что  2  миллиона  сельских  жителей  живут 
только  за  счет  приусадебного  хозяйства,  многие  находятся  в  условиях 
неполной  занятости. 

На  этапе  предварительного  анализа  данных  проводилось  выявление 
аномальных  отклонений  с  помощью  критерия  Ирвина: 
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У  г  -У  1-1 
8 

(  -У 

Уі  -  У 

ѵ  ) 

5 

где  8  -  несмещённое  среднеквадратическое  отклонение  данного  ряда, 

У  -  выборочное  среднее  ряда  или  средний  уровень  ряда. 

Далее  проверили  гипотезу  Н0:  «аномальные  данные  отсутствуют». 
По  таблицам  для  критерия  Ирвина  определили  Акритическое . 

Если  А,  <  ^критическое ,  то  у,  считается  нормальным; 
если  Аг  >  Акритическое,  то  у,  считается  аномальным. 

Результаты  расчетов  значений  критерия  Ирвина  для  всех  4-х 
регионов  Украины  представлены  в  табл.  4.2.,  критические  значения  этого 
критерия  -  в  табл.  4.3. 


Таблица  4.2 

Результаты  расчетов  значений  критерия  Ирвина  для  4-х  регионов 

Украины 


Наблюдение,  п 

Критерий  Ирвина,  А1 

1-й  регион 

2-й  регион 

3-й  регион 

4-й  регион 

1 

0,9201553854 

0,9408362265 

0,8724324502 

0,8106499608 

2 

0,2568567858 

0,3505195001 

0,0025785949 

0,0577362738 

3 

0,2568567858 

0,3505195001 

0,0025785949 

0,0577362738 

4 

0,1894750729 

0,1329008981 

0,1536842581 

0,0900685872 

5 

0,1894750729 

0,1329008981 

0,1536842581 

0,0900685872 

6 

1,0717723729 

1,0309085599 

1,2052352723 

1,1116157258 

7 

1,0717723729 

1,0309085599 

1,2052352723 

1,1116157258 

8 

0,1935064574 

0,0397114237 

0,1340869366 

0,0846798683 

9 

0,1935064574 

0,0397114237 

0,1340869366 

0,0846798683 

10 

0,1935064574 

0,0391819381 

0,1340869366 

0,0846798683 
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Таблица  4.3 

Критические  значения  критерия  Ирвина 


Число  наблюдений  и 

0 

Л критическое 

Р=0,95 

Р=0,99 

2 

2,8 

3,7 

3 

2,2 

2,9 

10 

1,5 

2,0 

20 

1,3 

1,8 

30 

1,2 

1,7 

50 

1,1 

1,6 

100 

1,0 

1,5 

400 

0,9 

1,3 

1000 

0,8 

1,2 

Согласно  табл.  4.3  получили  4<Ѵ,„е 


Таким  образом,  в 


рассматриваемых  нами  случаях  аномальные  данные  отсутствуют. 

Для  подтверждения  стационарности  исследуемых  временных  рядов 
провели  расширенный  тест  Дики-Фуллера  (АБР-тест),  где  нулевой 
гипотезой  (Н0)  являлся  тот  факт,  что  исследуемый  временной  ряд 
нестационарен  и  описывается  моделями  авторегрессии  первого  порядка. 
После  чего  вычисляли  значение  1-статистики  для  проверки  нулевой 


гипотезы:  { 


— — — —  •  (п  -  Ь) ,  где  и  -  число  наблюдений,  Ь  -  оператор 

А  -Ауі 

лагов  (задержки).  Полученное  значение  сравнивали  с  критическим  (Ікрих). 
Г ипотеза  о  нестационарности  временного  ряда  отвергали,  если  1;  <  Ікрит. 

В  табл.  4.4  приведены  значения  статистики  Дикки-Фуллера  (Іф),  а 
также  критические  значения  (Ікрит)  этой  статистики  на  1%-м,  5%-м  и  10%-м 
уровнях  значимости  для  1-го  региона  Украины. 

Как  видно  из  таблицы  результатов  АБР -теста  временной  ряд  для 
первого  региона  можно  признать  стационарными,  так  как  критическое 
значенея  при  уровне  значимости  а  =  0,01  превышает  расчетное  значение 
статистики  Дикки-Фуллера  (1ф  <  Ікрит). 
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Таблица  4.4 

Результаты  теста  Дикки-Фуллера  для  1  -го  региона 


Статистика  Дикки-Фуллера  (Ц) 

1-Статистика 

Значимость* 

1.558967 

0.4561 

Критические  значения 

(Крит)  • 

1 %  уровень 

4.582648 

5%  уровень 

3.320969 

1 0%  уровень 

2.801384 

Анализ  результатов  АБР-теста  для  Центрально-Южного, 
Центрально-Северного  и  Восточного  регионов  Украины  также  подтвердил 
предположение  о  стационарности  исследуемых  временных  рядов. 

Средние  и  дисперсию  временных  рядов  рассчитывали  по  формулам 

[54]: 


у  =  Ту(  з2  =-1,{уг -у)2 

(=1  и  п  (= 1 

Выборочную  автоковариацию  к- го  порядка  вычисляли  как: 


1  п~к  /  — \ 

С  к  =-Т{уг-у){уг+к  -х) 

п  і=і 

Количественно  автокорреляция  устанавливается  с  помощью 
линейного  коэффициента  корреляции  между  уровнями  изучаемого  ряда  и 
некоторым  лагом  т=к.  Величина  лага  определяет  порядок  коэффициента 
корреляции.  С  увеличением  лага  число  пар  уровней,  по  которым 
рассчитывается  коэффициент  автокорреляции,  уменьшается.  Считается 
целесообразным  использовать  коэффициенты  автокорреляции  с 
порядками:  к<п-4. 

Расчет  коэффициентов  АКФ  к-го  порядка  проводится  согласно 
формулы  [55]: 
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уг-  значения 


гк 


И(Уг 


і—к-\- 1 


У2к-і)(Уі-к  ~У2к) 


-Уік- 1)2  Т(у(-к 

і=к-\-  1 


гДе  ?2*-і  =— *-7*  Уг*  =— 

Н  Л  (=к+1  Н  к  {=к+ 1 

временного  ряда,  і,  к,  п-  порядки  взятых  значений. 

Коэффициент  автокорреляции  первого  порядка  можно  рассчитать 

как: 

=  соѵ(уг_! ;  )  =  ум  •  уг  -  уг_ і  •  у, 

СГѴ  •  СГѴ  /~ 2  2  /~ 2  2~ 

>г-'  г, 


]  п  |  /г-1 

где  уг  = - -Хл.Ум  = - -ІУг 

и  - 1  г-2  п-  1м 

Для  проверки  гипотезы  о  том,  что  наблюдаемые  данные  являются 
реализацией  «белого  шума»,  использовали  (^-статистику  Льюинга  -  Бокса, 
которую  рассчитывали  по  формуле  [56]: 

2 

~  /  \т  Г, 

С  =  и(л  +  2)Х^у,  (4.1) 

к=\п-к 

где  гк  -  автокорреляция  к-го  порядка, 

ш  -  число  проверяемых  лагов; 

и  -  число  наблюдений  (п=10). 

В  предположении,  что  ряд  представляет  собой  белый  шум, 
использовали  следующий  доверительный  интервал  для  отдельного 
коэффициента  автокорреляции: 


п 


п(п  +  2) 


к  2 

Х\ —а,т  ■ 


гк  + 


П 


п(п  +  2) 


к  2 

Х[—а,т 


где  Х\- а  т  —  квантиль  распределения  % 2  с  т-степенями  свободы  и 
уровнем  значимости  а. 
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Нулевая  гипотеза  в  (^-критерии  заключается  в  том,  что  ряд 
представляет  собой  белый  шум,  то  есть  является  чисто  случайным 
процессом.  Используется  стандартная  процедура  проверки:  если  расчетное 

значение  ^-статистики  больше  заданного  квантиля  распределения  %х__а  т  > 

то  нулевая  гипотеза  отвергается  и  признается  наличие  автокорреляции  до 
т-го  порядка  в  исследуемом  ряде. 

Этап  3.  Определение  параметров  моделей  (параметрическая 
идентификация). 

Так  как  исследуемые  ряды  являются  стационарными,  то  для 
прогноза  заболеваемости  психическими  расстройствами  непсихотического 
типа  сельского  населения  Украины  была  выбрана  модель  Бокса- 
Дженкинса  АКМА{р,  с/),  которая  обобщает  стационарные  линейные 
модели  авторегрессии  и  скользящего  среднего  (р  -  параметры 
авторегрессии,  ^  -  параметры  скользящего  среднего). 

Таким  образом,  временной  ряд  представили  в  виде  суммы  двух 
процессов:  авторегрессии  порядка  р,  и  скользящего  среднего  порядка  ц: 

Уі  =(Р\Уг-\  +-+<РрУі-р  +  */ -ѲХ8,_Х -...-Ѳд8{_д  +8{ 

9 

где  уг  -  зависимая  переменная  в  момент  времени  і\ 

(р{,  (р2,...(рк_  оцениваемые  параметры  (параметры  авторегрессии); 

Ѳі  -  коэффициенты  процесса  скользящего  среднего, 

8,  -  случайная  величина. 

Задача  заключалась  в  выборе  наилучшего  сочетания  порядка 
авторегрессии  р  и  порядка  скользящего  среднего  ц. 

При  попытке  построить  лучшую  модель  испытывались  разные 
порядки  авторегрессии  и  скользящего  среднего. 

Построение  моделей  прогнозирования  осуществили  с  помощью 
пакета  прикладных  программ  8Р88  21.  В  редакторе  данных  ввели 
информацию  о  количестве  случаев  заболеваемости  расстройствами 
психики  и  поведения  непсихотического  характера  сельского  населения  в 
4-х  регионах  Украины  (рис.  4.3). 
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Рисунок  4.3  -  Диалоговое  окно  8Р88 

В  меню  программы  выбрали  Апаіуге— »Рогеса8Ііп§— ►Сгеаіе 
МобеІ8...(рис.  4.4). 


Рисунок  4.4  -  Выбор  метода  прогнозирования 
В  появившемся  окне  «Тіте  8егіе8  Мобеіег»  (рис.  4.5)  в  поле 
«Оерепсіепі:  ѴагіаЫез»  внесли  зависимую  переменную  «регион  і».  Для 
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получения  моделей  по  остальным  3-м  регионам  соответственно  внесли  в 
качестве  зависимых  переменных  «регион_2»  или  «регион_3»  или 
«регион_4».  В  поле  «МеСІюсІ»  выбрали  метод  «АКІМА»  и  нажали  кнопку 
«Сгйегіа. . .» (рис.4.6). 


Рисунок  4.5  -  Диалоговое  окно  «Тіте  8егіев  Мосіеіег» 

В  открывшемся  окне  «Тіте  8егіе8  Мосіеіег:  АКІМА  Сгйегіа» 
(рис.  4.6)  задали  параметры  модели  АКМА  р  ид.  Для  этого  в  поле 
«Мои8еа80иа1»  напротив  «Аиіоге§ге88Іѵе  (р)»  двойным  нажатием 
графического  манипулятора  ввели  значение  параметра  р= 2  напротив 
«Моѵіи§  Аѵега§е  (ц)»  -  д=  \ .  Для  построения  моделей  Бокса- Дженкинса 
АКМА  (1,  2)  и  АКМА  (2,  2)  аналогичным  образом  задали  параметры  р  ид. 
Далее  нажали  кнопку  «СопПпие». 

Для  получения  расчетных  коэффициентов  моделей  АКМА  для 
первого  региона  в  диалоговом  окне  Тіте  8егіе8  Мосіеіег»  (рис.  4.5) 
перешли  на  вкладку  «8іаіІ8ііс8»  (рис.  4.7)  и  напротив  опции  «Рагатеіег 
езіітаіез»  установили  «галочку». 


97 


Рисунок  4.6  -  Диалоговое  окно  «Тіте  8егіев  Моёеіег:  АКІМА  Сгііегіа» 


Рисунок  4.7  -  Диалоговое  окно  «Тіте  8егіев  Моёеіег».  Вкладка  «Зіаіізіісв» 
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Полученные  расчетные  коэффициенты  моделей  АКМА  для  первого 
региона  представлены  в  табл.  4.5.  Для  остальных  регионов  расчетные 
значения  представлены  в  табл.  Г.1-Г.З  (Приложение  Г). 


Таблица  4.5 

Параметры  моделей  АКМА  для  первого  региона 


Модель 

Бокса- 

Дженкинса 

Значения 

Стандартная 
ошибка 
АКФ,  8Е 

1- 

критерий 

Значимость, 

Р. 

Константа 

7154,930 

155,058 

46,144 

0,0001 

АКМА  (2,  1) 

Коэффициенты 

АК 

Лаг  1 

1,656 

0,152 

10,929 

0,0001 

Лаг  2 

-0,943 

0,160 

-5,893 

0,0010 

Коэффициент 

МА 

Лаг  1 

0,997 

85,557 

0,120 

0,0090 

АКМА  (1,2) 

Константа 

7115,963 

652,394 

10,907 

0,0001 

Коэффициент 

АК 

Лаг  1 

0,672 

0,466 

1,443 

0,0199 

Коэффициенты 

МА 

Лаг  1 

-1,045 

8,817 

-0,119 

0,9100 

Лаг  2 

-0,060 

0,662 

-0,091 

0,0310 

АКМА  (2,  2) 

Константа 

7117,078 

377,626 

18,847 

0,0001 

Коэффициенты 

АК 

Лаг  1 

1,577 

0,583 

2,706 

0,0420 

Лаг  2 

-0,793 

0,432 

-1,837 

0,0126 

Коэффициенты 

МА 

Лаг  1 

0,025 

38,746 

0,001 

1,0000 

Лаг  2 

0,961 

39,330 

0,024 

0,0010 

Закономерность  исследуемого  процесса  находит  отражение  в 
наличии  определенных  статистических  свойств  остаточной  компоненты,  а 
именно:  независимости  уровней,  их  случайности,  соответствия 

нормальному  закону  распределения  и  равенства  нулю  средней  ошибки. 

Пусть  построена  некоторая  модель,  на  основании  которой  вычислено 

Л  Л 

значение  у( .  Тогда  е1=уі-уі  -  остаток  в  определении  уп  где  у{  - 
наблюдение. 

Для  того,  чтобы  модель  могла  считаться  хорошей,  необходимо 
доказать,  что  ряд  остатков  ег  является  случайным  и  подчиняется 
нормальному  закону  распределений.  Для  доказательства  используется: 

1)  проверка  равенства  нулю  математического  ожидания. 
Выдвигается  основная  гипотеза 


99 


Н0 :  /л  =  0,  то  есть  проверяется  е  =  О 

Для  проверки  этой  гипотезы  строится  критерий  или  случайная 
величина  Стьюдента: 


Ік=гк/8Е(гк), 


где  5Е(гк ) 


1  +  2- 


І=  1 
П 


-  стандартная  ошибка  для  гк. 


По  таблице  распределения  Стьюдента  при  заданном  а  =  0,05 
находим  ікр  =Нс(,п)'  Если  і  <  ікр  -  считают,  что  /л-0,  если  і  >  ікр  - 
считают,  что  /л  ^  0 . 

Оценим  полученные  результаты. 

Сравнив  полученные  значения  [-критерия  с  табличным  при  заданном 
уровне  значимости  а=0,05  для  модели  АКМА  (2,  1)  отвергаем  гипотезу  о 
незначимости  коэффициентов  модели  (Р<  0.05). 

2)  Проверка  наличия  или  отсутствия  автокорреляции  в  отклонении 
модели.  Независимость  остаточной  компоненты  означает  отсутствие 
автокорреляции  между  остатками.  Автокорреляция  означает,  что 
некоторое  значение  ег  зависит  от  одного  или  нескольких  своих 
предыдущих  значений. 

Тест  Дарбина-Уотсона  связан  с  проверкой  гипотезы  об  отсутствии 
автокорреляции  первого  порядка,  т.е.  автокорреляции  между  соседними 
остаточными  членами  ряда.  При  этом  статистика  определяется  по  формуле: 


-ег_х)2 


і=і 


Можно  показать,  что  величина  0~2-(1-г1),  где  /у  -  коэффициент 
автокорреляции  первого  порядка  (т.е.  парный  коэффициент  корреляции 
между  двумя  последовательностями  остатков  еь  е2,  ...  ,еп_і  и  е2,  е3,  ... ,  еп). 
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Из  (4.2)  видно,  что  близость  значения  статистики  И  к  нулю  означает 
наличие  высокой  положительной  автокорреляции  (коэффициент  Гі  близок  к 
единице);  близость  значения  статистики  И  к  четырем  означает  наличие 
высокой  отрицательной  автокорреляции  (коэффициент  ід  близок  к  минус 
единице).  Естественно,  в  случае  отсутствия  автокорреляции  значение 
статистики  Б  будет  близким  к  двум  (коэффициент  г,  не  сильно  отличается 
от  нуля). 

Применение  на  практике  критерия  Дарбина-Уотсона  основано  на 
сравнении  расчетного  значения  статистики  Б  с  пороговыми,  граничными 
значениями  Бь  и  Бу.  Граничные  значения  Бь  и  Бу,  зависящие  от  числа 
наблюдений  п,  количества  объясняющих  переменных  в  модели,  уровня 
значимости  а,  находятся  по  таблицам  с  критериальными  значениями. 

Алгоритм  выявления  автокорреляции  остатков  на  основе  критерия 
Дарбина-Уотсона  следующий. 

Выдвигается  гипотеза  Н0  об  отсутствии  автокорреляции  остатков. 
Пусть  альтернативная  гипотеза  состоит  в  наличии  в  остатках 
положительной  автокорреляции  первого  порядка.  Тогда  при  сравнении 
расчетного  значения  статистики  Б  с  Бь  и  Би  возможны  следующие 
варианты: 

-  если  Б  <  Яд,  то  существует  положительная  корреляция,  гипотеза 
Н0  отвергается; 

-  если  Б  >  4 -Оь,  то  наблюдается  отрицательная  автокорреляция, 
гипотеза  Н0  не  отвергается; 

-  если  <  Б  <  4  —  Б^ ,  то  корреляция  отсутствует,  гипотеза  Н0  не 
отвергается; 

-  имеются  области  неопределенности:  Б^  <  Б  <А/> 

4-Эи  <  Б  <  4  —  И ь 

• 

Данный  критерий  нельзя  использовать,  если  среди  объясняющих 
переменных  содержатся  легированные  значения  результативного 
показателя  (например,  он  не  применим  к  моделям  авторегрессии).  Таким 
образом,  можно  считать,  что  в  случае  отсутствия  автокорреляции  в 
остатках  расчетное  значение  статистики  «не  слишком  отличается»  от  2. 

Мы  рассчитали  интеграционную  статистики  Дарбина-Уотсона 
(ІБ\У).  Значение  статистики  вычисляется  по  формуле: 
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ю\ѵ 


±(Ѵ,  -г,_о2 

1=2 


І<Т, -У)2 

1=1 


где  У  -  среднее  значение  ряда. 

Рассчитанные  значения  Госстатистики  сравниваются  с 
критическими  и  в  случае  ГО\У  <  ГО\\У  гипотеза  о  стационарности  ряда 
принимается. 

Получили,  например,  для  ряда  первого  региона  значение 
интеграционной  статистики  Дарбина-Уотсона  ГО\У  ~  1,47.  Сравнивая  это 
значение  с  верхним  критическим  значением  статистики  ГО\Уу~1,6, 
принимаем  гипотезу  о  стационарности  ряда  остатков.  Следовательно, 
можно  сказать,  что  ряд  остатков  действительно  является  «белым  шумом», 
а  модель  точной.  Такое  же  поведение  остатков  наблюдается  и  для 
остальных  регионов. 

При  графическом  анализе  учитывали,  что  ЧАКФ  (частная 
автокорреляционная  функция)  и  АКФ  белого  шума  не  должны  выходить  за 
пределы  красной  штриховой  линии,  та  как  остатки  должны  быть  малыми  и 
не  должно  быть  закономерных  корреляций. 

Этап  4.  Выбор  лучшей  модели. 

При  выборе  параметров  модели  учитывали  показатели  ее  точности. 
Они  описывают  величины  случайных  ошибок,  полученных  при 
использовании  модели. 

Ошибка  прогноза  -  величина,  характеризующая  расхождение  между 
фактическим  и  прогнозным  значением  показателя  [57]. 

Абсолютная  ошибка  прогноза  определяется  по  формуле: 

\  =  У,~У, 

где  У,  -  прогнозное  значение  показателя; 

у,  -  фактическое  значение. 

Эта  характеристика  имеет  ту  же  размерность,  что  и  прогнозируемый 
показатель,  и  зависит  от  масштаба  измерения  уровней  временного  ряда. 
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Ошибки  прогнозов,  полученные  ретроспективно  характеризуют 
точность  применяемой  модели. 

При  сравнении  прогнозных  и  действительных  значений  ряда 
рассчитали  ошибки.  Среди  множества  типов  ошибок,  самой  подходящей 
является  абсолютная  относительная  ошибка  -АРЕ  [58]: 


АРЕк  = 


У  к  ~  У  к 
Ук 


где  ук  -  данные,  ук  -  прогноз. 

Среднее  значение  ошибки  АРЕ  обозначается  МАРЕ  [59]: 


1  п 

МАРЕ  =  — 

п  к=і 


У  к  ~  У  к 
Ук 


•100% 


Если  МАРЕ  <  10%,  то  прогноз  сделан  с  высокой  точностью,  10%  < 
20%  -  прогноз  хороший,  20%  <  МАРЕ  <  50%  -  прогноз 

удовлетворительный,  МАРЕ  >  50%  -  прогноз  плохой. 

МАРЕ  подчеркивает,  насколько  велики  ошибки  прогноза  в  сравнении 
с  действительными  значениями  ряда. 

Естественно  ожидать,  что  правильно  подобранный  метод  даст  малые 
ошибки  прогноза. 

Кроме  ясности,  МАРЕ  обладает  большей  устойчивостью  к  выбросам 
по  сравнению  со  средней  квадратичной  ошибкой. 

Среднеквадратическая  ошибка  прогноза  (Меап  Зциагесі  Еггог,  М8Е) 
используется  в  качестве  оценки  метода  прогнозирования.  Поскольку 
каждое  значение  отклонения  возводится  в  квадрат,  то  подчеркиваются 
большие  ошибки  прогноза: 


М5Е  =  -Х(Ук-Ук)2 

п  *=, 


Следовательно,  прогнозы  на  основе  МАРЕ  будут  более  надежны. 
Модель  ряда  считается  хорошей,  если  ошибка  прогноза  МАРЕ  мала  и 
остатки  ряда:  остатки  =  данные  -  прогноз  имеют  вид  белого  шума,  т.е.  их 
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дисперсия  постоянна  и  они  некоррелированы  [60  -  63]. 

Качество  моделей  сравнивалось  по  информационному  критерию 
Шварца  ( 8ВСс  )■ 

Байесовский  критерий  Шварца  (Могтаіігеё  ВІС)  для  моделей 
регрессии  можно  представить  в  виде: 


8ВСс  =  1п(<52)  +  (—)■/? 

п 

где  п  -  объём  выборочной  совокупности; 

р  -  число  параметров  в  модели  . 

^  (М8Е”)  -  оценка  максимального  правдоподобия  дисперсии 

остатков  е,  модели  регрессии. 

Выбирается  та  модель,  для  которой  значение  8ВСс  является 
наименьшим. 

В  пакете  8Р88  для  получения  показателей  качества  моделей  АКМА 
(1,  2),  (2,  1)  и  (2,  2)  в  диалоговом  окне  «Тіте  8егіев  Моёеіег».  Вкладка 
«Зіайвйсв»  (рис.  4.7)  выбрали  основные  критерии  их  оценки,  а  именно  «К 
8 с]  иге»,  «Кооі  теап  8циаге  еггог»,  «Меап  аЬзоІиІе  регсеп1;а§е  еггог», 
«Махітит  аЬзоІиІе  регсепіа^е  еггог»,  «Ыогтаіігесі  ВІС».  Результаты  для 
первого  региона  приведены  в  табл.  4.6. 


Таблица  4.6 

Основные  критерии  оценки  качества  моделей  АКМА 


Модель 

Критерии  оценки  качества 

к- 

М8Е 

МАРЕ 

МахАРЕ 

8ВСс 

АКМА  (1,2) 

0,804 

470,577 

4,033 

11,343 

13,229 

АКМА  (2,  2) 

0,835 

472,470 

3,841 

8,655 

13,467 

АКМА  (2,  1) 

0,836 

430,458 

4,297 

8,646 

13,051 
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Коэффициент  детерминации  К2  во  всех  трех  моделях  выше  0,8,  т.е. 
модели  работают  хорошо  (имеют  высокую  значимость).  Анализируя 
полученные  значения  М8Е,  МАРЕ,  МахАРЕ  и  8ВСс  для  четырех 
регионов  Украины  можно  сделать  вывод,  что  модель  АРМА  (2,1)  является 
лучшей. 

Таким  образом,  для  прогнозирования  заболеваемости  сельского 
населения  НПР  для  1-4  регионов  предлагаем  использовать  следующие  4 
модели  временного  ряда: 

-  для  первого  региона 

уг  =  7154.930  +  1.656- ум  -0.943 -уг_2  +  0.997 -*м; 

-  для  второго  региона 

уг  =  7713.684  + 1.677 -ум  -  0.952  •  уг_2  +  0.993  •  ег_х ; 

-  для  третьего  региона 

уг  =  6040. 129  +  0.745  •  у м  -  0. 1 16  •  уг_2  -  0.993  •  е м ; 

-  для  четвертого  региона 

уг  =  3154.232  +  0.780- ум  -0.084- уг_2  -0.999 -*м. 

Полученные  модели  можно  признать  удачными,  все  коэффициенты  в 
модели  являются  статистически  значимыми  (Р  <  0.05),  в  модели 
присутствует  автокорреляция,  остатки  моделей  стационарны  и  имеют 
нормальное  распределение. 

Оценку  АКФ  и  ЧАКФ  остатков,  получили  в  8Р88  (рис.  4.8)  путем 
установки  в  диалоговом  окне  «Тіте  8егіе§  Мобеіег».  Вкладка  «Р1о1»опций 
«Рейсіиаі  аиШсоггеІаНоп  ІипсНоп  (АСЕ)»  и  «Кезіёиаі  рагйаі  аиіосоггеіайоп 
Етсйоп  (РАСЕ)». 
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АКФ  остатков  ЧАКФ  остатков 


Рисунок  4.8  -  АКФ  и  ЧАКФ  остатков  для  1-го  региона 

Из  рисунка  следует,  что  ряд  остатков  является  «белым  шумом»,  а 
построенная  модель  прогноза  для  первого  региона  является  точной. 

Модели  прогноза  для  остальных  регионов  также  являются  точными 
рис.  Г1  -  ГЗ  (Приложение  Г). 

Анализ  и  верификация  прогностических  моделей  показали,  что 
точность  прогнозирования  достаточно  высокая,  полученные  результаты 
приемлемы  для  принятия  управленческих  решений  в  процедуре  выбора  и 
формирования  лечебно-  профилактических  мероприятий 

Этап  6.  Прогнозирование  с  помощью  построенной  модели. 

Было  проведено  прогнозирование  заболеваемости  непсихотическими 
психическими  расстройствами  седьского  населения  Украины  с 
использованием  синтезированных  моделей. 

В  пакете  8Р88  для  получения  показателей  прогнозирования  по 
заболеваемости  непсихотическими  психическими  расстройствами  по 
каждому  региону  в  диалоговом  окне  «Тіше  8егіе§  Мобеіег».  Вкладка 
«Зіаіівіісв»  (рис.  4.7)  установили  флажок  напротив  опции  «Бізріау 
Ібгесазіз». 

Результаты  прогнозирования  по  заболеваемости  непсихотическими 
психическими  расстройствами  представлены  в  табл.  4.7. 
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Для  построения  графиков  прогноза  заболеваемости  сельского 
населения  непсихотическими  психическими  расстройствами  по  всем  4-м 
регионам  Украины  в  диалоговом  окне  «Тіте  8егіев  Моёеіег»  (рис.  4.9) 
перешли  на  вкладку  «Р1оЬ>  и  установили  отметку  в  соответствующем  поле 
«Ріоіз  іог  іпсіі ѵісіиаі  Мосіеіз».  После  этого  нажали  кнопку  «ОК». 


Рисунок  4.9  -  Диалоговое  окно  «Тіте  8егіе§  Мобеіег».  Вкладка  «Ріоі» 

Из  полученных  графиков  (рис.  4.10  -  4.13)  видно,  что  прогнозные 
кривые  практически  повторяют  фрагменты  кривой  исходного  ряда,  причём 
все  наблюдаемые  значения  ряда  попадают  в  доверительный  интервал. 


—  Наблюдаемые  значения 

- Значения  модели 

Нижняя  граница  д.и. 
Врхняя  граница  д.и. 
Прогноз 


Рисунок  4.10  -  Прогноз  заболеваемости  сельского  населения  Украины 
непсихотическими  психическими  расстройствами  в  первом  регионе 
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Таблица  4.7 

Результаты  прогнозирования  по  заболеваемости  непсихотическими 
психическими  расстройствами  для  четырех  регионов  Украины 


Регио 

н 

Год 

2013 

2014 

2015 

1 

Верхний 

доверительный 

интервал 

8115,71 

7782,70 

7411,10 

7267,90 

7507,82 

8044,50 

Прогноз 

7303,62 

6769,51 

6376,50 

6229,16 

6355,62 

6703,95 

Нижний 

доверительный 

интервал 

6491,53 

5756,33 

5341,91 

5190,42 

5203,43 

5363,41 

2 

Верхний 

доверительный 

интервал 

8973,37 

8599,03 

8108,83 

7793,88 

7864,12 

8301,01 

Прогноз 

8159,07 

7572,57 

7052,92 

6740,17 

6710,67 

6959,06 

Нижний 

доверительный 

интервал 

7344,78 

6546,12 

5997,01 

5686,45 

5557,22 

5617,12 

Верхний 

доверительный 

интервал 

7739,20 

8566,31 

8769,49 

8796,40 

8785,53 

8772,56 

Прогноз 

6586,51 

6341,42 

6201,01 

6124,92 

6084,58 

6063,38 

Нижний 

доверительный 

интервал 

6433,83 

4916,52 

4632,53 

4653,45 

4683,63 

4854,20 

і 

Верхний 

доверительный 

интервал 

3615,77 

3457,89 

3352,32 

3283,24 

3238,22 

3208,91 

Прогноз 

4333,66 

4865,15 

5013,87 

5045,54 

5041,86 

5029,81 

Нижний 

доверительный 

интервал 

2897,89 

2050,63 

2690,77 

1920,94 

2434,58 

2388,01 

Таким  образом,  на  основании  полученных  результатов  следует 
отметить,  что  заболеваемость  непсихотическим  психическими 
расстройствами  у  сельского  населения  на  протяжении  2011  -  2012  гг. 
носила  отрицательный  характер,  то  есть  наметился  спад  уровня 
заболеваемости  по  всем  регионам  Украины. 
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—  Наблюдаемые  значения 

—  Значения  модели 
Нижняя  граница  д.и. 
Врхняя  граница  д.и. 
Прогноз 


Рисунок  4.11  -  Прогноз  заболеваемости  сельского  населения  Украины 
непсихотическими  психическими  расстройствами  во  втором  регионе  на  2015  год 


—  Наблюдаемые  значения 

—  Значения  модели 
Нижняя  граница  д.и. 
Врхняя  граница  д.и. 
Прогноз 


Рисунок  4.12  -  Прогноз  заболеваемости  сельского  населения  Украины 
непсихотическими  психическими  расстройствами  в  третьем  регионе  на  2015  год 


—  Наблюдаемые  значения 

- Значения  модели 

Нижняя  граница  д.и. 
Врхняя  граница  д.и. 
Прогноз 


Рисунок  4.13.  -  Прогноз  заболеваемости  сельского  населения  Украины 
непсихотическими  психическими  расстройствами  в  четвертом  регионе  на  20 1 5  год 


109 


Для  оптимизации  медицинского  обеспечения  лиц,  страдающих 
непсихотическими  психическими  расстройствами,  необходимо 
пользоваться  прогностическими  моделями  заболеваемости,  на  основе 
которых  должно  проводиться  планирование  ресурсов  и  объемов 
медицинской  помощи,  а  также  дальнейшее  усовершенствование 
психиатрической  службы,  направленное  на  повышение  качества, 
территориальной  и  экономической  доступности  психиатрической  помощи 
сельскому  населению  Украины. 


ПО 


ПРИЛОЖЕНИЕ  А 


Таблица  А.  1 

Кодировочная  таблица  к  задаче  «Способ  диагностики  хронической 
сердечной  недостаточности  у  подростков  на  ранней  стадии  развития» 


№ 

Показатель 

Описание 

Код  показателя 

1 

ёіа§по2 

Диагноз 

1-  миокардиофиброз 

2-  дизпластическая  кардимиоопатия 

3-  гипертрофическая  кардимиопатия 

4-  дилятационная  кардимиопатия 

5-  аритмия 

6-  синдром  перенапряжения  миокарда 

7-  коарктация  аорты 

8-  стеноз  аорты 

9-  стеноз  легочной  артерии 

10-  дефект  межжелудочковой  перегородки 

11-  дефект  межпредсердной  перегородки 

1 2-  сочетанные  пороки 

13-  тетрада  Фалло 

14-  здоров 

2 

Оіа§по2І 

Диагноз 

1-  миокардиофиброз 

2-  дизпластическая  кардиомиопатия 

3-  пороки  с  перегрузкой  левых  отделов 
сердца 

4-  пороки  с  перегрузкой  правых  отделов 
сердца 

5-  порки  по  смешанному  типу 

6-  синдром  перенапряжения  миокарда 

7-  условно  здоровые 

3 

8И 

Сердечная 

недостаточность 

1-  без  сердечной  недостаточности 

2-  1-я  стадия  сердечной  недостаточности 

3-  контрольная  группа 

4 

роі 

Пол 

1-  М 

2-  Ж 

5 

іт_8гЬ 

ИММУНОЛОГИЯ 
С-реактивный  белок 

0  -  отрицательный 

1-  Положительный  (определяется) 

6 

іт_а8І_о 

ИММУНОЛОГИЯ 

антистрептолизин-О; 

абс.ед. 

0  -  норма  (  0  -  200) 

1  -  более  200 
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Рисунок  А.  1  -  Общая  ковариационная  матрица  Н 
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Рисунок  А.2  -  Внутригрупповая  ковариационная  матрица  А 
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Рисунок  А.З  -  Объединенная  внутригрупповая  матрица  <2 


Семенная  недостаточность 

ЦИТОКИНЫ 
ІИ.  пг/мл 

ЦИТОКИНЫ 
И6,  пг/мл 

ЦИТОКИНЫ 
ФНО  а,  пг/мл 

СД-95,  % 
клеток 
апоптотичны 

X 

ЦАМФ, 

циклический 

уденозинмон 

офосфат, 

нмоль/мл 

!{ 

КАРБОНИЛИ 
РОВАННЫЕ 
БЕЛКИ,  у.е. 

НОРАДРЕНА 

ЛИН, 

нмоль/сут 

ЭХОКАРДИО 

ГРАФИЯ 

ПОСЛЕ 

НАГРУЗКИ 

ФВ.  % 

ИКДО, 

индекс  конечно 
го  диастоличес 
кого  объема 
лев.жел„у.е. 

Без  сердечной 
недостаточности 

ЦИТОКИНЫ  И 1 ,  пг/мл 

872,800 

15,798 

242,098 

-14,008 

281,766 

■6,358 

,454 

80,305 

186,777 

63,549 

ЦИТОКИНЫ  И  6,  пг/мл 

15,798 

17,066 

15,316 

3,915 

5,830 

-,524 

-.073 

7,975 

13,875 

-10,408 

ЦИТОКИНЫ  ФНО  а, 
пг/мл 

242,098 

15,316 

364,255 

10,996 

130,568 

-3,203 

-2,467 

-4,229 

98,819 

30,286 

СД-95,  %  клеток 
апоптотичных 

-14,008 

3,915 

10,996 

19,169 

10,250 

-1,139 

-,039 

46,894 

-4,865 

-11,062 

цАМФ,  циклический 

уденозинмонофосфат, 

нмоль/мл 

281,766 

5,830 

130,568 

10,250 

198,114 

-3,910 

-.324 

26,904 

54,290 

50,200 

РЕНИН,  нг/млЛчас 

-6,358 

-.524 

-3,203 

-1,139 

-3,910 

,469 

-,040 

-2,311 

-,260 

,878 

КАРБОНИЛИРОВАННЫЕ 
БЕЛКИ,  у.е. 

,454 

-.073 

-2,467 

-.039 

-.324 

-,040 

,301 

-1,351 

,975 

,457 

НОРАДРЕНАЛИН, 

нмоль/сут 

80,305 

7,975 

-4,229 

46,894 

26,904 

-2,311 

-1,351 

562,787 

37,902 

-65,347 

ЭХОКАРДИОГРАФИЯ 
ПОСЛЕ  НАГРУЗКИ  ФВ, 

% 

186,777 

13,875 

98,019 

■4,865 

54,290 

■,260 

,975 

37,902 

157,019 

■22,583 

ИКДО,  индекс  конечного 
диастолического  объема 

63,549 

■10,408 

30,286 

■11,062 

50,200 

,878 

,457 

■65,347 

■22,583 

149,372 

лев.  жел..  ѵ.е. 

1  стадия  сердечной 
недостаточности 

ЦИТОКИНЫ  ІИ,  пг/мл 

35,004 

■2,504 

25,026 

,832 

24,794 

■,184 

-,017 

4,525 

11,004 

■21,981 

ЦИТОКИНЫ  И  6,  пг/мл 

-2,504 

39,515 

■4,956 

4,218 

■7,048 

,649 

,051 

■11,015 

7,628 

-,256 

ЦИТОКИНЫ  ФНО  а, 
пг/мл 

25,026 

■4,956 

181,164 

4,505 

38,202 

1,910 

,048 

62,922 

5,215 

■32,262 

СД-95,  %  клеток 
апоптотичных 

,832 

4,218 

4,505 

10,476 

7,840 

■,464 

■,074 

34,513 

1,615 

■6,067 

цАМФ,  циклический 

уденозинмонофосфат, 

нмоль/мл 

24,794 

-7,048 

38,202 

7,840 

138,709 

-2,457 

,783 

70,012 

4,512 

■24,597 

РЕНИН,  нг/млАчас 

-,184 

,649 

1,910 

-,464 

■2,457 

,576 

,071 

4,930 

■,858 

1,814 

КАРБОНИЛИРОВАННЫЕ 
БЕЛКИ,  у.е. 

-,017 

,051 

,048 

-,074 

,783 

,071 

,213 

-3,696 

,781 

-,067 

НОРАДРЕНАЛИН, 

нмоль/сут 

4,525 

-11,015 

62,922 

34,513 

70,012 

4,930 

-3,696 

2688,478 

-128,801 

158,038 

ЭХОКАРДИОГРАФИЯ 
ПОСЛЕ  НАГРУЗКИ  ФВ, 

% 

11,004 

7,628 

5,215 

1,615 

4,512 

-,858 

,781 

-128,801 

78,973 

-26,215 

ИКДО,  индекс  конечного 
диастолического  объема 

-21,981 

-.256 

-32,262 

-6,067 

-24,597 

1,814 

-,067 

158,038 

-26,215 

240,147 

лев.  жел..  ѵ.е. 

Контрольная  группа 

ЦИТОКИНЫ  ІИ,  пг/мл 

4,818 

,981 

-,046 

,668 

2,437 

-,034 

,024 

-1,027 

1,174 

4,349 

ЦИТОКИНЫ  И  6,  пг/мл 

,981 

2,529 

-,062 

,075 

-.165 

,011 

,052 

-1,322 

-1,193 

2,748 

ЦИТОКИНЫ  ФНО  а, 
пг/мл 

-,046 

-.062 

2,655 

,339 

-,359 

-,081 

,069 

-3,611 

,365 

-,238 

СД-95,  %  клеток 
апоптотичных 

,668 

,075 

,339 

10,683 

-4,597 

-.111 

,068 

7,986 

1,651 

,665 

цАМФ,  циклический 

уденозинмонофосфат, 

нмоль/мл 

2,437 

-.165 

-,359 

■4,597 

61,440 

-,270 

,009 

,824 

1,975 

13,008 

РЕНИН,  нг/млАчас 

-,034 

,011 

-,081 

-.111 

-.270 

,114 

,006 

,835 

•,018 

-,105 

КАРБОНИЛИРОВАННЫЕ 
БЕЛКИ,  у.е. 

,024 

,052 

,069 

,068 

,009 

,006 

,167 

,830 

-,277 

,630 

НОРАДРЕНАЛИН, 

нмоль/сут 

■1,027 

■1,322 

■3,611 

7,986 

,824 

,835 

,830 

730,897 

61,595 

12,039 

ЭХОКАРДИОГРАФИЯ 
ПОСЛЕ  НАГРУЗКИ  ФВ, 

1,174 

■1,193 

,365 

1,651 

1,975 

■,018 

-,277 

61,595 

47,093 

■14,794 

ИКДО,  индекс  конечного 
диастолического  объема 

4,349 

2,748 

■,238 

,665 

13,008 

-,105 

,630 

12,039 

-14,794 

90,408 

лев.  жел..  ѵ.е. 

Рисунок  А.4  -  Групповая  ковариационная  матрица  К 
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Таблица  А.  2 

Фактическая  и  прогнозируемые  группы,  значения  дискриминантных 

функций 


Наивероятнейшая  группа 

Вторая 

вероятнейшая 

группа 

Значения 
дискриминантн 
ых  функций 

Номер 

наблюдения 

Фак¬ 

тичес¬ 

кая 

группа 

Предо 

казан 

ная 

груп¬ 

па 

Р 

ст.с 

в 

Р(0=§ 

Ю=ф 

Квадрат 

расстоя¬ 

ния 

Махалон 
обиса  до 
центра 

Г 

РУ 

п 

па 

Р(0=§ 

т=б) 

Квадрат 

расстоя¬ 

ния 

Махалон 
обиса  до 
центра 

Функ¬ 
ция  1 

Функ¬ 
ция  2 

1 

2 

3 

4 

5 

6 

7 

Т 

9 

10 

11 

12 

Исход  1 

1 

1 

,634 

2 

,987 

,911 

2 

,013 

11,351 

3,040 

2,534 

ные  2 

1 

1 

,015 

2 

,851 

8,440 

2 

,115 

14,206 

-,123 

,502 

3 

1 

1 

,437 

2 

,842 

1,657 

2 

,158 

6,764 

1,704 

1,149 

4 

1 

1 

,390 

2 

,969 

1,883 

2 

,031 

10,500 

1,088 

1,446 

5 

1 

1 

,290 

2 

,883 

2,475 

2 

Д17 

8,280 

1,229 

1,063 

6 

1 

2** 

,260 

2 

,794 

2,692 

1 

,206 

3,631 

3,188 

,811 

7 

1 

1 

,144 

2 

,680 

3,875 

2 

,320 

7,143 

1,167 

,640 

8 

1 

1 

,223 

2 

,612 

3,000 

2 

,388 

5,666 

3,625 

1,555 

9 

1 

1 

,004 

2 

1,000 

11,270 

2 

,000 

43,840 

2,355 

5,721 

10 

1 

1 

,688 

2 

,996 

,748 

2 

,004 

13,676 

2,829 

2,837 

11 

1 

1 

,435 

2 

,706 

1,664 

2 

,294 

5,177 

2,962 

1,416 

12 

1 

1 

,014 

2 

1,000 

8,604 

2 

,000 

34,716 

3,244 

5,074 

13 

1 

1 

,483 

2 

,861 

1,455 

2 

,139 

6,863 

3,155 

1,792 

14 

1 

2** 

,271 

2 

,975 

2,612 

1 

,025 

8,162 

1,830 

-,470 

15 

1 

1 

,074 

2 

,876 

5,202 

2 

,124 

10,867 

,504 

,744 

16 

1 

1 

,452 

2 

,959 

1,589 

2 

,041 

9,643 

3,356 

2,289 

17 

1 

1 

,000 

2 

1,000 

27,464 

3 

,000 

50,053 

-1,342 

6,327 

18 

1 

2** 

,260 

2 

,885 

2,692 

1 

,115 

5,019 

2,070 

,134 

19 

1 

1 

,070 

2 

1,000 

5,327 

2 

,000 

31,201 

,309 

3,831 

20 

1 

1 

,255 

2 

,947 

2,730 

2 

,053 

10,245 

3,751 

2,367 

21 

1 

1 

,061 

2 

1,000 

5,584 

2 

,000 

26,829 

3,373 

4,364 

22 

1 

1 

,176 

2 

1,000 

3,475 

2 

,000 

25,732 

2,547 

4,184 

23 

1 

1 

,924 

2 

,998 

,157 

2 

,002 

14,148 

1,748 

2,560 
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1 

2 

3 

1  4 

5 

6 

7 

I  8 

1  9 

10 

11 

12 

24 

1 

1 

Д12 

2 

1,000 

4,371 

2 

,000 

27,027 

2,792 

4,347 

25 

1 

1 

,217 

2 

,570 

3,054 

2 

,430 

5,379 

1,642 

,687 

26 

1 

1 

,584 

2 

,990 

1,075 

2 

,010 

11,957 

3,104 

2,629 

27 

1 

1 

,414 

2 

1,000 

1,766 

2 

,000 

22,802 

1,517 

3,569 

28 

2 

2 

,697 

2 

,986 

,721 

1 

,014 

7,453 

2,753 

-,276 

29 

2 

2 

,236 

2 

,990 

2,887 

1 

,010 

10,362 

4,621 

,375 

30 

2 

2 

,669 

2 

,977 

,803 

1 

,023 

6,513 

3,054 

,008 

31 

2 

2 

,883 

2 

,994 

,248 

1 

,006 

8,860 

3,366 

-,319 

32 

2 

2 

,622 

2 

1,000 

,950 

1 

,000 

17,205 

2,964 

-1,683 

33 

2 

2 

,005 

2 

,956 

10,423 

1 

,034 

15,317 

,188 

-1,042 

34 

2 

2 

,499 

2 

1,000 

1,389 

1 

,000 

19,511 

3,053 

-1,939 

35 

2 

2 

,644 

2 

,975 

,880 

1 

,025 

6,441 

3,241 

,108 

36 

2 

2 

,057 

2 

,999 

5,723 

1 

,001 

17,879 

1,211 

-1,760 

37 

2 

2 

,868 

2 

1,000 

,284 

1 

,000 

15,766 

3,504 

-1,339 

38 

2 

2 

,264 

2 

1,000 

2,661 

1 

,000 

23,121 

4,858 

-1,564 

39 

2 

2 

,818 

2 

1,000 

,401 

1 

,000 

15,684 

3,227 

-1,422 

40 

2 

2 

,801 

2 

,990 

,444 

1 

,010 

7,839 

3,031 

-,266 

41 

2 

2 

,142 

2 

,919 

3,901 

1 

,081 

7,013 

1,528 

-,212 

42 

2 

2 

,439 

2 

,999 

1,649 

1 

,001 

14,250 

2,264 

-1,397 

43 

2 

2 

,104 

2 

,908 

4,524 

1 

,092 

7,342 

4,455 

1,036 

44 

2 

2 

,746 

2 

,997 

,587 

1 

,003 

10,417 

3,948 

-,271 

45 

2 

2 

,599 

2 

,999 

1,025 

1 

,001 

13,708 

4,383 

-,540 

46 

2 

2 

,687 

2 

,993 

,751 

1 

,007 

8,799 

2,583 

-,552 

47 

2 

2 

,878 

2 

1,000 

,260 

1 

,000 

15,412 

3,764 

-1,181 

48 

2 

2 

,748 

2 

1,000 

,581 

1 

,000 

17,504 

3,890 

-1,407 

49 

2 

2 

,881 

2 

1,000 

,253 

1 

,000 

15,559 

3,679 

-1,240 

50 

2 

2 

,746 

2 

,992 

,586 

1 

,008 

8,562 

3,668 

-,095 

51 

2 

2 

,108 

2 

,989 

4,452 

1 

,011 

11,607 

4,993 

,578 

52 

2 

2 

,894 

2 

,994 

,225 

1 

,006 

8,848 

3,207 

-,386 

53 

2 

2 

,510 

2 

,999 

1,346 

1 

,001 

14,511 

4,539 

-,551 

54 

2 

2 

,344 

2 

1,000 

2,134 

1 

,000 

23,308 

3,396 

-2,276 

55 

2 

2 

,656 

2 

1,000 

,843 

1 

,000 

17,298 

4,249 

-1,186 

56 

2 

2 

,962 

2 

,999 

,077 

1 

,001 

13,313 

3,666 

-,921 

57 

2 

2 

,546 

2 

1,000 

1,210 

1 

,000 

16,866 

4,499 

-,958 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

58 

2 

2 

,618 

2 

1,000 

,963 

1 

,000 

14,660 

4,383 

-,699 

59 

2 

2 

,821 

2 

,991 

,395 

1 

,009 

7,967 

3,119 

-,258 

60 

2 

2 

,813 

2 

,991 

,415 

1 

,009 

8,052 

3,347 

-,174 

61 

2 

2 

,722 

2 

,983 

,651 

1 

,017 

7,035 

3,215 

-,032 

62 

2 

2 

,583 

2 

1,000 

1,079 

1 

,000 

16,441 

4,441 

-,936 

63 

2 

2 

,464 

2 

1,000 

1,535 

1 

,000 

16,761 

4,648 

-,830 

64 

2 

2 

,096 

2 

,962 

4,678 

1 

,038 

9,407 

4,769 

,866 

65 

2 

2 

,330 

2 

,999 

2,219 

1 

,001 

14,957 

2,082 

-1,493 

66 

2 

2 

,815 

2 

,998 

,409 

1 

,002 

11,478 

2,782 

-,944 

67 

2 

2 

,561 

2 

1,000 

1,156 

1 

,000 

20,436 

3,895 

-1,774 

68 

2 

2 

,894 

2 

,994 

,224 

1 

,006 

8,853 

3,211 

-,385 

69 

2 

2 

,521 

2 

1,000 

1,303 

1 

,000 

20,216 

4,272 

-1,562 

70 

2 

2 

,073 

2 

,996 

5,223 

1 

,004 

14,662 

1,185 

-1,344 

71 

2 

2 

,667 

2 

,986 

,811 

1 

,014 

7,563 

3,593 

,066 

72 

2 

2 

,438 

2 

,998 

1,653 

1 

,002 

12,372 

2,165 

-1,143 

73 

2 

2 

,473 

2 

,997 

1,499 

1 

,003 

11,077 

2,192 

-,953 

74 

2 

2 

,852 

2 

,993 

,321 

1 

,007 

8,402 

3,081 

-,353 

75 

2 

2 

,337 

2 

,934 

2,174 

1 

,066 

5,702 

3,724 

,625 

76 

2 

2 

,890 

2 

1,000 

,232 

1 

,000 

15,135 

3,762 

-1,142 

77 

2 

2 

,286 

2 

1,000 

2,505 

1 

,000 

25,115 

4,274 

-2,141 

78 

2 

2 

,903 

2 

,999 

,203 

1 

,001 

12,558 

3,022 

-1,047 

79 

2 

2 

,557 

2 

1,000 

1,171 

1 

,000 

16,926 

2,766 

-1,686 

80 

2 

2 

,683 

2 

1,000 

,763 

1 

,000 

14,846 

4,282 

-,800 

81 

2 

2 

,929 

2 

,999 

,148 

1 

,001 

13,310 

3,790 

-,858 

82 

2 

2 

,356 

2 

1,000 

2,065 

1 

,000 

19,743 

4,807 

-1,149 

83 

2 

2 

,832 

2 

,991 

,369 

1 

,009 

8,114 

3,098 

-,293 

84 

2 

2 

,820 

2 

1,000 

,397 

1 

,000 

13,990 

2,966 

-1,264 

85 

2 

2 

,181 

2 

,870 

3,416 

1 

,130 

5,468 

1,780 

,058 

86 

2 

2 

,495 

2 

1,000 

1,408 

1 

,000 

20,408 

3,258 

-1,992 

87 

2 

2 

,332 

2 

,879 

2,206 

1 

Д21 

4,411 

2,428 

,300 

88 

2 

2 

,441 

2 

,969 

1,636 

1 

,031 

6,784 

2,274 

-,224 

89 

2 

2 

,511 

2 

,997 

1,343 

1 

,003 

10,906 

4,306 

-,082 

90 

2 

2 

,422 

2 

1,000 

1,723 

1 

,000 

18,311 

4,706 

-1,019 

91 

2 

2 

,122 

2 

1,000 

4,207 

1 

,000 

27,214 

2,849 

-2,788 
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1 

2 

1  3 

1  4 

5 

6 

|  7 

8 

9 

10 

И 

12 

92 

2 

2 

,020 

2 

1,000 

7,818 

1 

,000 

37,771 

5,309 

-2,866 

93 

3 

3 

,767 

2 

1,000 

,530 

1 

,000 

53,056 

-4,825 

Д12 

94 

3 

3 

,705 

2 

1,000 

,698 

1 

,000 

52,702 

-4,855 

,287 

95 

3 

3 

,574 

2 

1,000 

1,111 

1 

,000 

41,028 

-3,403 

-,906 

96 

3 

3 

,669 

2 

1,000 

,803 

1 

,000 

55,658 

-5,008 

,102 

97 

3 

3 

,922 

2 

1,000 

,163 

1 

,000 

45,035 

-4,258 

,224 

98 

3 

3 

,791 

2 

1,000 

,469 

1 

,000 

36,772 

-3,503 

,052 

99 

3 

3 

,787 

2 

1,000 

,480 

1 

,000 

37,548 

-3,461 

-,202 

100 

3 

3 

,764 

2 

1,000 

,539 

1 

,000 

47,650 

-4,534 

,462 

101 

3 

3 

,057 

2 

1,000 

5,717 

1 

,000 

24,635 

-1,817 

-,676 

102 

3 

3 

,804 

2 

1,000 

,437 

1 

,000 

53,787 

-4,809 

-,090 

103 

3 

3 

,789 

2 

1,000 

,473 

1 

,000 

40,018 

-3,544 

-,486 

104 

3 

3 

,608 

2 

1,000 

,996 

1 

,000 

58,511 

-5,147 

-,078 

105 

3 

3 

,689 

2 

1,000 

,745 

1 

,000 

35,243 

-3,291 

-Д14 

106 

3 

3 

,839 

2 

1,000 

,352 

1 

,000 

46,454 

-4,415 

,367 

107 

3 

3 

,619 

2 

1,000 

,958 

1 

,000 

59,014 

-5,130 

-,227 

108 

3 

3 

,490 

2 

1,000 

1,428 

1 

,000 

34,188 

-3,001 

-,485 

109 

3 

3 

,982 

2 

1,000 

,036 

1 

,000 

44,151 

-3,997 

-,271 

110 

3 

3 

,630 

2 

1,000 

,923 

1 

,000 

34,417 

-3,461 

,500 

111 

3 

3 

Д71 

2 

1,000 

3,534 

1 

,000 

41,242 

-4,289 

1,709 

112 

3 

3 

,124 

2 

1,000 

4,168 

1 

,000 

22,313 

-2,189 

,391 

ИЗ 

3 

3 

,450 

2 

1,000 

1,599 

1 

,000 

54,672 

-4,244 

-1,427 

114 

3 

3 

,836 

2 

1,000 

,359 

1 

,000 

44,137 

-4,252 

,425 

115 

3 

3 

,608 

2 

1,000 

,995 

1 

,000 

55,519 

-5,049 

,271 

116 

3 

3 

,783 

2 

1,000 

,488 

1 

,000 

53,678 

-4,832 

-,002 

117 

3 

3 

,825 

2 

1,000 

,384 

1 

,000 

40,472 

-3,956 

,422 

118 

3 

3 

,169 

2 

1,000 

3,557 

1 

,000 

29,447 

-2,353 

-,729 

119 

3 

3 

,581 

2 

1,000 

1,087 

1 

,000 

59,965 

-5,193 

-,232 

120 

3 

3 

,862 

2 

1,000 

,297 

1 

,000 

50,921 

-4,650 

,056 

121 

3 

3 

,512 

2 

1,000 

1,340 

1 

,000 

48,687 

-4,710 

,849 

122 

3 

3 

,928 

2 

1,000 

,149 

1 

,000 

50,213 

-4,536 

-,116 

123 

3 

3 

,899 

2 

1,000 

,212 

1 

,000 

51,536 

-4,613 

-,173 

124 

3 

3 

,570 

2 

1,000 

1,125 

1 

,000 

55,419 

-5,070 

,367 

125 

3 

3 

,725 

2 

1,000 

,643 

1 

,000 

36,532 

-3,628 

,441 

126 

3 

,844 

2 

1,000 

,340 

1 

,000 

53,320 

-4,593 

-,547 
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ш 

2 

3 

I  4 

5 

1  6 

7 

8 

|  9 

|  10 

11 

12 

3 

3 

,962 

2 

1,000 

,077 

1 

,000 

42,260 

-3,970 

,043 

в 

3 

3 

,741 

2 

1,000 

,601 

1 

,000 

55,155 

-4,920 

-,055 

129 

3 

3 

,712 

2 

1,000 

,679 

1 

,000 

48,741 

-4,628 

,507 

130 

3 

3 

,699 

2 

1,000 

,715 

1 

,000 

56,787 

-4,998 

-,159 

131 

3 

3 

,996 

2 

1,000 

,009 

1 

,000 

44,918 

-4,079 

-,225 

132 

3 

3 

,365 

2 

1,000 

2,017 

1 

,000 

49,605 

-4,825 

1,085 

133 

3 

3 

,455 

2 

1,000 

1,573 

1 

,000 

64,029 

-5,318 

-,630 

134 

3 

3 

,027 

2 

1,000 

7,256 

2 

,000 

26,068 

-1,680 

-1,234 

135 

3 

3 

,939 

2 

1,000 

,126 

1 

,000 

47,650 

-4,168 

-,521 

136 

3 

3 

,533 

2 

1,000 

1,257 

1 

,000 

32,664 

-3,335 

,601 

137 

3 

3 

,828 

2 

1,000 

,378 

1 

,000 

45,907 

-3,921 

-,735 

138 

3 

3 

,735 

2 

1,000 

,616 

1 

,000 

56,711 

-4,899 

-,409 

139 

3 

3 

,708 

2 

1,000 

,692 

1 

,000 

35,085 

-3,349 

,049 

140 

3 

3 

,366 

2 

1,000 

2,008 

1 

,000 

33,455 

-2,825 

-,661 

141 

3 

3 

,922 

2 

1,000 

,162 

1 

,000 

45,017 

-4,257 

,224 

142 

3 

3 

,703 

2 

1,000 

,705 

1 

,000 

55,990 

-4,983 

-,042 

143 

3 

3 

,955 

2 

1,000 

,093 

1 

,000 

41,534 

-3,883 

-,024 

144 

3 

3 

,940 

2 

1,000 

,124 

1 

,000 

42,025 

-3,989 

,146 

145 

3 

3 

,486 

2 

1,000 

1,443 

1 

,000 

56,755 

-4,463 

-1,326 

146 

3 

3 

,815 

2 

1,000 

,409 

1 

,000 

38,584 

-3,774 

,350 

147 

3 

3 

,949 

2 

1,000 

,105 

1 

,000 

42,691 

-3,860 

-,306 

148 

3 

3 

,134 

2 

1,000 

4,020 

1 

,000 

69,352 

-4,998 

-1,984 

149 

3 

3 

,513 

2 

1,000 

1,334 

1 

,000 

59,310 

-5,269 

,131 

150 

3 

3 

,365 

2 

1,000 

2,014 

1 

,000 

56,616 

-4,305 

-1,577 

151 

3 

3 

,505 

2 

1,000 

1,365 

1 

,000 

48,385 

-3,817 

-1,285 

152 

3 

3 

,822 

2 

1,000 

,393 

1 

,000 

41,988 

-4,091 

,458 

153 

3 

3 

,227 

2 

1,000 

2,967 

1 

,000 

63,987 

-4,733 

-1,788 

154 

3 

3 

,568 

2 

1,000 

1,132 

1 

,000 

55,683 

-4,460 

-1,184 

155 

3 

3 

,278 

2 

1,000 

2,560 

1 

,000 

26,745 

-2,596 

,206 

156 

3 

3 

,592 

2 

1,000 

1,050 

1 

,000 

35,587 

-3,164 

-,435 

157 

3 

3 

,731 

2 

1,000 

,626 

1 

,000 

56,716 

-4,833 

-,571 

158 

3 

3 

,574 

2 

1,000 

1,109 

1 

,000 

53,694 

-4,981 

,484 

159 

3 

3 

,937 

2 

1,000 

,130 

1 

,000 

48,584 

-4,245 

-,515 

**.  Неправильно  классифицированное  наблюдение 
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ПРИЛОЖЕНИЕ  Б 


Таблица  Б.1 

Кодировочная  таблица  для  определения  систолической  дисфункции 

миокарда  у  подростков 


№ 

Показатель 

Описание 

Код  показателя 

1 

ёіа§пог 

Диагноз 

1-  миокардиофиброз 

2-  дизпластическая  кардимиопатия 

3-  гипертрофическая  кардиомипатия 

4-  дилятационная  кардимиопатия 

5-  аритмия 

6-  синдром  перенапряжения  миокарда 

7-  коарктация  аорты 

8-  стеноз  аорты 

9-  стеноз  легочной  артерии 

10-  дефект  межжелудочковой  перегородки 

1 1-  дефект  межпредсердной  перегродки 

12-  сочетанные  пороки 

13-  тетрада  Фалло 

14-  здоров 

2 

8П 

Сердечная 

недостаточность 

1  -  без  сердечной  недостаточности 

2-  1  стадия  сердечной  недостаточности 

3-  2  стадия  сердечной  недостаточности 

4-  контрольная  группа 

3 

гагпіса 

Физическая  нагрузка 

1-  адекватная  реакция  на  физическую 
нагрузку 

2-  неадекватная  реакция  а  физическую 
нагрузку 

3-  норма 

4 

роі 

Пол 

1-  м 

2-  ж 

5 

іт_8гЬ 

ИММУНОЛОГИЯ 
С-реактивный  белок 

0  -  отрицательный 

1  -  положительный  (  определяется) 

6 

іт_а8І_о 

ИММУНОЛОГИЯ 

антистрептолизин-О; 

абс.ед. 

0  -  норма  (  0  -  200) 

1  -  более  200 

121 


ПРИЛОЖЕНИЕ  В 


Таблица  В.1 

Таблица  сопряженности  для  критерия  Хосмера- Л емешова 


Группы  =  без  ишемии 
миокарда 

Группы  =  ишемическая 
дилятационная 
кардиомиопатия 

Наблюденные 

(°у) 

Ожидаемые 

(^) 

Наблюденные 

(°у) 

Ожидаемые 

(^) 

Всего 

(?т) 

Шаг  1  1 

5 

4,219 

0 

,781 

5 

2 

12 

15,226 

9 

5,774 

21 

3 

9 

5,162 

3 

6,838 

12 

4 

0 

1,095 

6 

4,905 

6 

5 

0 

,276 

9 

8,724 

9 

6 

0 

,020 

8 

7,980 

8 

7 

0 

,002 

9 

8,998 

9 

8 

0 

,000 

9 

9,000 

9 

9 

0 

,000 

8 

8,000 

8 

Шаг  2  1 

9 

8,500 

0 

,500 

9 

2 

8 

7,387 

1 

1,613 

9 

3 

7 

6,271 

2 

2,729 

9 

4 

0 

3,289 

9 

5,711 

9 

5 

2 

,535 

7 

8,465 

9 

6 

0 

,018 

10 

9,982 

10 

7 

0 

,000 

9 

9,000 

9 

8 

0 

,000 

9 

9,000 

9 

9 

0 

,000 

9 

9,000 

9 

10 

0 

,000 

5 

5,000 

5 

Шаг  3  1 

9 

8,997 

0 

,003 

9 

2 

10 

9,798 

0 

,202 

10 

3 

6 

6,290 

4 

3,710 

10 

4 

1 

,915 

8 

8,085 

9 

5 

,000 

9 

9,000 

9 

6 

,000 

6 

6,000 

6 

7 

,000 

34 

34,000 

34 

Шаг  4  1 

9 

9,000 

0 

,000 

9 

2 

10 

9,918 

0 

,082 

10 

3 

6 

6,672 

3 

2,328 

9 

4 

1 

,411 

9 

9,589 

10 

5 

0 

,000 

7 

7,000 

7 

6 

0 

,000 

42 

42,000 

42 
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Итеоаиия 

-2  1_од 

Правдоподоб 

ие 

Коэффициенты  1 

Сопзіапі 

САА 

ППТ 

ИЛ  8 

ЛПВП 

Шаг  1  1 

99,840 

,501 

,320 

2 

80,335 

-.150 

1,132 

3 

66,404 

-.795 

2,330 

4 

58,004 

-1,432 

3,923 

5 

54,502 

-1,961 

5,577 

6 

53,657 

-2,308 

6,809 

7 

53,594 

-2,426 

7,255 

8 

53,593 

-2,437 

7,295 

9 

53,593 

-2,437 

7,295 

Шаг  2  1 

91,945 

-4,213 

,278 

2,51  4 

2 

71,381 

-8,221 

,996 

4,409 

3 

55,973 

-11,916 

2,198 

6,057 

4 

46,433 

-14,951 

3,868 

7,321 

5 

41,133 

-16,989 

5,968 

8,054 

6 

38,668 

-19,586 

8,376 

9,088 

7 

38,070 

-21,827 

10,200 

10,037 

8 

38,029 

-22,631 

10,828 

10,385 

9 

38,029 

-22,700 

10,880 

10,415 

10 

38,029 

-22,700 

10,880 

10,415 

Шаг  3  1 

89,849 

-4,471 

,298 

2,726 

-.002 

2 

67,861 

-8,904 

1,040 

4,878 

-.003 

3 

50,105 

-13,668 

2,350 

7,129 

-.004 

4 

38,692 

-18,655 

4,217 

9,482 

-.006 

5 

31,796 

-22,986 

6,697 

11,469 

-.008 

6 

26,226 

-28,760 

10,777 

1  4,025 

-.010 

7 

21,883 

-38,646 

17,460 

18,503 

-.013 

8 

19,066 

-54,641 

27,470 

25,888 

-.017 

9 

17,983 

-72,916 

37,827 

34,428 

10 

17,804 

-84,136 

44,008 

39,676 

11 

17,798 

-86,668 

45,402 

40,859 

-.026 

12 

17,798 

-86,764 

45,456 

40,904 

-.026 

13 

17,798 

-86,765 

45,456 

40,904 

-,026 

Шаг  4  1 

87,894 

-2,516 

,315 

2,609 

-.002 

-1,634 

2 

66,268 

-6,905 

1,022 

4,980 

-.002 

-1,993 

3 

48,269 

-11,794 

2,332 

7,516 

-.004 

-2,356 

4 

36,946 

-16,668 

4,181 

10,165 

-.005 

-2,948 

5 

30,402 

-21,219 

6,554 

12,399 

-.008 

-3,123 

6 

24,474 

-27,586 

10,778 

1  4,972 

-.010 

-2,624 

7 

18,896 

-38,958 

18,668 

20,393 

-.012 

-3,115 

8 

14,276 

-60,809 

32,726 

31,253 

-.019 

-4,509 

9 

11,840 

-90,575 

50,581 

46,292 

-.027 

-6,672 

10 

11,048 

-118,135 

67,132 

60,278 

-.034 

-8,843 

11 

10,943 

-132,178 

75,793 

67,405 

-.037 

-10,023 

12 

10,940 

-134,738 

77,442 

68,699 

-.038 

-10,250 

13 

10,940 

-134,818 

77,496 

68,739 

-.038 

-10,257 

14 

10,940 

-134,818 

77,496 

68,739 

-.038 

-10,257 

Рисунок  В.  1  —  История  итераций 
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ПРИЛОЖЕНИЕ  Г 


Лаг  Корр  Станд 
откл. 

1  +,795  ,2739  ■ 


2  +,426  ,2532  - 


3  +,024  ,2415  ■ 


4  -,336  ,2236  - 


5  -,486  ,2041  - 


6  -,423  ,1826  • 


0 

-1 


<3  Р 


8,42  ,0037 


-  11,  15  ,0038 


11,16  ,0109 


13,42  ,0094 


■  19,09  ,0019 


24,46  ,0004 


Доверительный 


-0.5 


0.0 


0.5 


1.0 


о  интервал 


Лаг  Корр  Станд. 

откч. 

1  +,795  ,3162 

2  -,556  ,3162 

3  -,242  ,3162 

4  -,254  ,3162 

5  +,256  ,3162 

6  +,087  ,3162 

О 

-1,0  -0,5  0,0  0,5  1,0 


Дов  ері  п  епьнын 
ннтервач 


б) 


Рисунок  Г.1  -Автокорреляционные  функции  остатков  для  2-го  региона: 


а  -  АКФ; б  - ЧАКФ 
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Лаг  Корр.  Станд. 

откт. 

1  +,753  ,2739  - 


2  +,363  ,2582  • 


3  +,034  ,2415  - 


4  -,268  ,2236  - 


5  -,449  ,2041  - 


6  -,407  ,1826  - 


<3  Р 


7,56  ,0060 


9,53  ,0085 


9,55  ,0228 


10,  98  ,0268 


15,81  ,0074 


20,77  ,0020 


■1,0 


-0,5 


0,0 


0,5 


Доверитеть  ный 
о  интерват 


1,0 


Лаг  Корр.  Станд. 

откл. 

1  +,753  ,3162 

2  -,471  ,3162 

3  -,039  ,3162 

4  -,356  ,3162 

5  +,001  ,3162 

6  +,127  ,3162 


0  I— 

-1,0 


а) 


Доверительный 

интервал 


б) 


Рисунок  Г. 2  -  Автокорреляционные  функции  остатков  для  3-го  региона: 

а  -  АКФ; б  - ЧАКФ 
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Лаг  Корр 


Станд. 

откл. 


1  +,759  ,3162  ■ 


2  -,423  ,3162  ■ 


3  -,093  ,3162  ■ 


4  -,347  ,3162  ■ 


5  -,058  ,3162  ■ 


6  +,174  ,3162  • 


0 
-1,0 


а) 


Доверительный 


-0,5 


0,0 


0,5 


1  д  іштервал 


б) 


Рисунок  Г.З  -  Автокорреляционные  функции  остатков  для  4-го  региона: 

а  -  АКФ; б  - ЧАКФ 
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Таблица  Г.1 


Параметры  моделей  АКМА  для  второго  региона 


Модель  Бокса- 
Дженкинса 

Значения 

Стандартная 

ошибка 

АКФ,  8Е 

і- 

критерий 

Значимость, 

Р 

Константа 

7713,684 

176,831 

43,622 

0,0001 

АКМА  (2,  1) 

Коэффици¬ 
енты  АК 

га 

1,677 

0,122 

13,702 

0,0001 

-0,952 

0,133 

-7,164 

0,0001 

Коэффици¬ 
ент  МА 

■ 

0,993 

37,845 

0,026 

0,0098 

АКМА  (1,2) 

Константа 

7732,531 

756,785 

10,218 

0,0001 

Коэффици¬ 
ент  АК 

га 

0,702 

0,452 

1,553 

0,171 

-1,020 

2,784 

-0,366 

0,727 

Коэффици¬ 
енты  МА 

а 

-0,056 

0,693 

-0,081 

0,938 

АКМА  (2,  2) 

Константа 

7722,477 

304,076 

25,397 

0,0001 

Коэффици¬ 
енты  АК 

га 

1,653 

0,248 

6,673 

0,001 

-0,909 

0,231 

-3,926 

0,011 

Коэффици¬ 
енты  МА 

га 

0,721 

29,883 

0,024 

0,982 

га 

0,263 

8,599 

0,031 

0,977 
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Таблица  Г.2 


Параметры  моделей  АКМА  для  третьего  региона 


Модель 

Бокса- 

Дженкинса 

Значения 

Стандарт¬ 
ная  ошибка 

АКФ,  8Е 

1-критерий 

Значимость, 

Р 

6040,129 

675,520 

8,941 

0,0001 

АКМА  (2,  1) 

мЯнЯЛгЯ 

ЭД  ■  ■  И  в)  иряі  зі 

га 

0,745 

0,522 

1,428 

0,0203 

-0,116 

0,549 

-0,212 

0,0083 

■ 

-0,993 

20,456 

-0,049 

0,0430 

АКМА  (1,2) 

6045,998 

700,926 

8,626 

0,0001 

'  'ЯЯ'.'.  «зд 

ши 

0,650 

0,554 

1,173 

0,  0285 

-0,891 

1,718 

-0,519 

0,623 

ж 

0,031 

0,940 

0,033 

0,975 

АКМА  (2,  2) 

6001,893 

556,726 

10,781 

0,0001 

ши 

1,558 

0,937 

1,662 

0,0157 

ш 

-0,725 

0,618 

-1,173 

0,0294 

0,00032 

57,444 

0,0000005 

0,0100 

ши 

0,991 

66,784 

0,015 

0,0890 
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Таблица  Г.З 

Параметры  моделей  АКМА  для  четвертого  региона 


Модель 

Бокса- 

Дженкинса 

Значения 

Стандартная 

ошибка 

АКФ,  8Е 

1- 

критерий 

Значимость, 

Р. 

Константа 

3154,232 

499,656 

6,313 

0,0010 

АКМА  (2,  1) 

га 

0,780 

0,541 

1,441 

0,0200 

-0,084 

0,575 

-0,146 

0,0090 

Коэффици¬ 
ент  МА 

■ 

-0,999 

182,897 

-0,005 

0,0060 

АКМА  (1,2) 

Константа 

3153,977 

553,845 

5,695 

0,001 

Коэффици¬ 
ент  АК 

га 

0,711 

0,443 

1,606 

0,159 

-1,046 

2,334 

-0,448 

0,670 

Коэффици¬ 
енты  МА 

■ 

-0,099 

0,659 

-0,150 

0,886 

АКМА  (2,  2) 

Константа 

3139,145 

386,215 

8,128 

0,0001 

Коэффици¬ 
енты  АК 

шл 

1,628 

0,825 

1,974 

0,0105 

ш 

-0,779 

0,594 

-1,312 

0,0247 

Коэффици¬ 
енты  МА 

0,005 

52,787 

0,0009 

1,000 

ш 

0,985 

57,475 

0,017 

0,0070 
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